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 40 [0] [0] 26 pyk Pyruvate kinase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 48921130 (0.000)26 (0.660)
+GCGGCCGC
26/534 0.8 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

AGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGC  >  AE015451/4893238‑4893773
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
gctcagcctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcaccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGc                                                                                                                                                                                                                                                    <  2:179208‑M2/6‑1 (MQ=255)
      cctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTg                                                                                                                                                                                                                                       <  1:157632‑M2/19‑1 (MQ=255)
           ccgaaagcctgaaaggcctgccgggcgtttacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGtt                                                                                                                                                                                                                                        <  2:305690‑M2/18‑1 (MQ=255)
              aaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacggcttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCt                                                                                                                                                                                                                                >  2:105728‑M2/275‑300 (MQ=255)
               aagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgccccaacgtcaacgacttccgccccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGgc                                                                                                                                                                                                                                             >  2:110810‑M2/274‑286 (MQ=255)
                     gaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGcccc                                                                                                                                                                                                                                 <  1:110814‑M2/25‑1 (MQ=255)
                     gaaaggcctgccgggcgggtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCggg                                                                                                                                                                                                                                                <  1:387311‑M2/10‑1 (MQ=255)
                                cgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGc                                                                                                                                                                                                              >  2:171956‑M2/257‑300 (MQ=255)
                                  ggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTcc                                                                                                                                                                                                           <  1:334334‑M2/47‑1 (MQ=255)
                                         acgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccgttgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCcac                                                                                                                                                                                                        <  1:105732‑M2/50‑1 (MQ=255)
                                           gccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCCCCCCCACAGGGc                                                                                                                                                                                                   >  2:138873‑M2/246‑300 (MQ=255)
                                                     gtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCCCCCACAGGGCCACGCTGGc                                                                                                                                                                                          >  1:86828‑M2/236‑299 (MQ=255)
                                                           ccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAg                                                                                                                                                                                     <  1:171962‑M2/69‑1 (MQ=255)
                                                                ggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAgg                                                                                                                                                                                                                  <  2:74258‑M2/40‑1 (MQ=255)
                                                                ggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAgg                                                                                                                                                                                                                  >  1:74262‑M2/225‑264 (MQ=255)
                                                                        tggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacaacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACCCAAGGCCGCTccc                                                                                                                                                                                                          >  2:194331‑M2/217‑264 (MQ=255)
                                                                        tggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacaacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTccc                                                                                                                                                                                                          <  1:194338‑M2/48‑1 (MQ=255)
                                                                         ggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCtgtg                                                                                                                                                                    <  1:138878‑M2/86‑1 (MQ=255)
                                                                               gaataaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAggg                                                                                                                                                                                   <  1:420677‑M2/71‑1 (MQ=255)
                                                                               gaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACCAGGCCGCTCCCACAGGGCCACGCt                                                                                                                                                                                             >  2:420665‑M2/210‑270 (MQ=255)
                                                                                                  tcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCTCGTTGCGCCCCTTTCGCGACACAAGGCCGCCCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCt                                                                                                                                                          >  1:308900‑M2/191‑286 (MQ=255)
                                                                                                                    acgccagcgccgaagccttgggcctgtcggccagcgacgggctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGc                                                                                                                                                                              <  2:86824‑M2/76‑1 (MQ=255)
                                                                                                                                   ccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGc                                                                                                                                                                                          >  2:30777‑M2/158‑221 (MQ=255)
                                                                                                                                   ccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGc                                                                                                                                                                                          <  1:30780‑M2/64‑1 (MQ=255)
                                                                                                                                                                                              gcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgccaaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCt                                                                                                                                                          <  2:308892‑M2/96‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    atgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGCGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCACAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGc  >  2:321525‑M2/13‑260 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
AGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGC  >  AE015451/4893238‑4893773

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