Predicted mutation
evidence seq id position mutation annotation gene description
MC JC AE015451 4,892,114 1412 bp→GCGGCCGC coding (5‑1416/1416 nt) pyk → Pyruvate kinase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ AE015451 4892114 4893525 1412 69 [0] [0] 68 pyk Pyruvate kinase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 48921130 (0.000)66 (0.640)
+GCGGCCGC
60/584 0.5 100% coding (1275/1275 nt)
coding (4/1416 nt)
ttuD
pyk
putative hydroxypyruvate reductase
Pyruvate kinase
?AE015451 4893526 = 0 (0.000)intergenic (+1/‑207) pyk/PP_4302 Pyruvate kinase/Urea transporter

GTCAACTACAGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTGCGCAA  >  AE015451/4893229‑4893818
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                    
ctccccacacgaggcccttccgatctgcctgaaaggccggccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGAtt                                                                                                                                                                                                                                                                                                   <  2:406606‑M2/4‑1 (MQ=255)
  agttcctgctcagcctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGc                                                                                                                                                                                                                                                                                                 >  3:528732‑M2/296‑301 (MQ=255)
          ctcggcctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGgcg                                                                                                                                                                                                                                                                                         <  1:209813‑M2/14‑1 (MQ=255)
           tcagcctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgacggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGgcgc                                                                                                                                                                                                                                                                                        <  3:348054‑M2/15‑1 (MQ=255)
                     cgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGcccc                                                                                                                                                                                                                                                                              >  4:312860‑M2/277‑301 (MQ=255)
                          gcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcgattagatcggaagagcacacgtctgaact                                                                                                                                                                                                                                                                         >  4:406719‑M2/272‑275 (MQ=255)
                                      tgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGcc                                                                                                                                                                                                                                                             >  4:259014‑M2/260‑301 (MQ=255)
                                              ggtgacgccctggccggtgacacccacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCcac                                                                                                                                                                                                                                                     <  2:96132‑M2/50‑1 (MQ=255)
                                                ctccgccctggcagttgcaaccgccggcagcgatgtcccggaagaaaacgccagtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaaccacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAg                                                                                                                                                                                                                                                   <  2:480610‑M2/52‑1 (MQ=255)
                                                 tacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAgg                                                                                                                                                                                                                                                  <  3:120632‑M2/53‑1 (MQ=255)
                                                   cgccctggccggtgacaccgacggcatcgatggctcgggagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGc                                                                                                                                                                                                                                                <  2:258945‑M2/55‑1 (MQ=255)
                                                    gccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGcc                                                                                                                                                                                                                                               <  3:174280‑M2/56‑1 (MQ=255)
                                                      cctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCAc                                                                                                                                                                                                                                             >  4:746447‑M2/244‑301 (MQ=255)
                                                      ccgggccgctgaccccgcccgcatcgatggctcgcaagcaaacgccggttccttcatgacccccgccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCAc                                                                                                                                                                                                                                             <  4:24745‑M2/58‑1 (MQ=255)
                                                           ccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTgg                                                                                                                                                                                                                                        <  1:164037‑M2/63‑1 (MQ=255)
                                                            gggtgaacccgagggcaccgatggctcggaagaaaaggccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGc                                                                                                                                                                                                                                       <  2:312779‑M2/64‑1 (MQ=255)
                                                            acgtgtgctcttccgatctgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGc                                                                                                                                                                                                                                       <  3:8294‑M2/64‑1 (MQ=255)
                                                              gtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGcac                                                                                                                                                                                                                                     <  1:329817‑M2/66‑1 (MQ=255)
                                                                 acaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACgctggccccag                                                                                                                                                                                                                                  >  2:656581‑M2/233‑296 (MQ=255)
                                                                    ccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGa                                                                                                                                                                                                                               <  4:837419‑M2/72‑1 (MQ=255)
                                                                    ccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatcagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGCCCCCCCCCCCAGGGCCCCGCTGGCACAagggg                                                                                                                                                                                                                               >  2:288167‑M2/230‑300 (MQ=255)
                                                                      tggagttcagacgtgtgctcttccgatctcggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgggactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAAt                                                                                                                                                                                                                             <  3:412737‑M2/74‑1 (MQ=255)
                                                                        cggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACCGGGCCACGCTGGCCCAAGGGAATGc                                                                                                                                                                                                                           >  4:403422‑M2/226‑301 (MQ=255)
                                                                         ggcatcgatggctcggaaaaaaacgccggggccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCa                                                                                                                                                                                                                          <  2:746190‑M2/77‑1 (MQ=255)
                                                                          gcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAg                                                                                                                                                                                                                         <  4:656802‑M2/78‑1 (MQ=255)
                                                                               gatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTcccacagggccacgctggcagatcggaagagcggcgtg                                                                                                                                                                                                                    >  1:8293‑M2/219‑283 (MQ=255)
                                                                                             aaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCCCGACACAAGGCCCCTCCCACAGGGCCACGCTCGCCCAAGggcctgcagtccctggtggcgcgccctc                                                                                                                                                                                                      >  2:464922‑M2/205‑288 (MQ=255)
                                                                                                   cggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGTCCGCtcccacagggccacgctggcacaagggaatagatcggaagagcgtcgtgtagggaaaga                                                                                                                                                                                                >  1:412726‑M2/199‑272 (MQ=255)
                                                                                                        ccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGAtt                                                                                                                                                                                           >  1:73192‑M2/194‑301 (MQ=255)
                                                                                                            catgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCCGGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCAACGGGCCCCGCTGGCACAAGGGAATGCAGTCCCTTTGGTAGCGGCcgtgtgtcgcggatggag                                                                                                                                                                                       >  2:776907‑M2/190‑294 (MQ=255)
                                                                                                                     ggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCGCGCTTGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTCGTCTCGCCATTggagggcgacgcc                                                                                                                                                                              >  2:273111‑M2/181‑295 (MQ=255)
                                                                                                                       ccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGccct                                                                                                                                                                            <  4:465060‑M2/123‑1 (MQ=255)
                                                                                                                              cgccggcgccgaggcctggggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGaa                                                                                                                                                                     <  2:403309‑M2/130‑1 (MQ=255)
                                                                                                                              cgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGaa                                                                                                                                                                     >  3:676031‑M2/172‑301 (MQ=255)
                                                                                                                                ccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATa                                                                                                                                                                   <  4:511911‑M2/132‑1 (MQ=255)
                                                                                                                                          agccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGCATGCAGTCCCTGTTGGACCGGCCTTGTGTCGCGATTGGCGGGCAAAGCCCTCCCTTAataccgcccctcc                                                                                                                                                         >  2:703999‑M2/160‑295 (MQ=255)
                                                                                                                                           gccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgccttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAACGCCCTCCCTGAATACCGCTCTAccc                                                                                                                                                        <  1:528719‑M2/143‑1 (MQ=255)
                                                                                                                                             cttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGa                                                                                                                                                      >  4:68165‑M2/157‑301 (MQ=255)
                                                                                                                                               tgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGAcc                                                                                                                                                    <  4:777170‑M2/147‑1 (MQ=255)
                                                                                                                                                      gtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGaaa                                                                                                                                             <  4:273180‑M2/154‑1 (MQ=255)
                                                                                                                                                        ctgcccgcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAAc                                                                                                                                           <  4:716870‑M2/156‑1 (MQ=255)
                                                                                                                                                        cggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAAc                                                                                                                                           >  4:27348‑M2/146‑301 (MQ=255)
                                                                                                                                                          gccagcgacgagctggacaacacaaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCg                                                                                                                                         <  2:828864‑M2/158‑1 (MQ=255)
                                                                                                                                                          gccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGGCCCTGTGGGCACGGCCTTGGTTCGCGATTGGAGGGCAAAGCCTTCCCTTAATACCGCCCTACCCGCCCGgggaaaagca                                                                                                                                         >  2:716627‑M2/144‑298 (MQ=255)
                                                                                                                                                                             acaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCAcc                                                                                                                      <  4:288243‑M2/177‑1 (MQ=255)
                                                                                                                                                                              aaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCg                                                                                                                     <  1:676017‑M2/178‑1 (MQ=255)
                                                                                                                                                                                           ggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                                                                                        <  4:704241‑M2/191‑1 (MQ=255)
                                                                                                                                                                                            gctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAACCCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAcc                                                                                                       >  4:162286‑M2/110‑301 (MQ=255)
                                                                                                                                                                                                  cccccgcgctggatggactgaccggcaccaaaccgacccgaaccaccgtcaacgccttccgcgccttcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCact                                                                                                 <  2:68146‑M2/198‑1 (MQ=255)
                                                                                                                                                                                                         gctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCCACCTCCACTGACTGcc                                                                                          >  3:346217‑M2/97‑301 (MQ=255)
                                                                                                                                                                                                          ctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGTAAACCGATGCCCAAAGAGGCCCACCGCTCCCCCAtcaacccccaccgccccccc                                                                                         >  2:291977‑M2/96‑292 (MQ=255)
                                                                                                                                                                                                                         ctccccgaaccgccccgcaccaacgtcaacgacttccgcgccatcctgatccttgggactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaaaa                                                                          <  3:73193‑M2/221‑1 (MQ=255)
                                                                                                                                                                                                                         ctcaccgaaccgacccgcacccacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGCTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaaaa                                                                          <  1:346207‑M2/221‑1 (MQ=255)
                                                                                                                                                                                                                           caccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAACGAGGCCCACCGCTCCCCCATccacctccactgactgccccctgcaccccacaactt                                                                        >  2:567820‑M2/79‑285 (MQ=255)
                                                                                                                                                                                                                              cgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTtcc                                                                     >  1:197978‑M2/76‑300 (MQ=255)
                                                                                                                                                                                                                                aaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTc                                                                   >  1:13238‑M2/74‑301 (MQ=255)
                                                                                                                                                                                                                                            ccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGcccc                                                       >  1:572719‑M2/62‑301 (MQ=255)
                                                                                                                                                                                                                                              aacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCTCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCcgc                                                     >  1:253768‑M2/60‑300 (MQ=255)
                                                                                                                                                                                                                                                 ttcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTg                                                  <  3:197985‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                                                        cttaccgcgccgtccatgtcctttggactacccaagcatgagcggccggGATTGCCGGGGGCGCGTTGCGCCCCTTGCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACg                                           <  2:27334‑M2/252‑1 (MQ=255)
                                                                                                                                                                                                                                                                gccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTAGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTacc                                   >  1:20916‑M2/42‑300 (MQ=255)
                                                                                                                                                                                                                                                                    ccctgatccttgggactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCTCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAAc                               <  3:253773‑M2/264‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCCCCCCCATCCACCTCCACTGACGGCCCCATGTACACCAAACATTTCGTCAACCCGTGCCCCGCCTGGgccccggcgtcaccccacg                              >  2:748270‑M2/37‑292 (MQ=255)
                                                                                                                                                                                                                                                                            cttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGgcgttactcaagatcggaag                       >  4:156612‑M2/30‑292 (MQ=255)
                                                                                                                                                                                                                                                                                agactgcccaatcatgagcgaccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAg                   >  3:670931‑M2/26‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           tcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCCCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGAATGGGCACCGGCGTTACCCAACGGCTTCAGCCCGgtgctgctgcc        >  4:642917‑M2/15‑300 (MQ=255)
                                                                                                                                                                                                                                                                                              ggcgcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTgcg     <  4:748527‑M2/290‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 gcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTGCGCaa  <  4:292053‑M2/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                    
GTCAACTACAGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTGCGCAA  >  AE015451/4893229‑4893818

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: