Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,950,350 Δ1,672 bp [ilvL][ilvG] [ilvL], ilvX, ilvG, [ilvG]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3950350 3952021 1672 15 [0] [0] 15 [ilvL]–[ilvG] [ilvL], ilvX, ilvG, [ilvG]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 39503490 (0.000)14 (0.720) 10/286 0.7 100% coding (28/99 nt) ilvL ilvG operon leader peptide
?NC_000913 3952022 = 0 (0.000)pseudogene (481/663 nt) ilvG pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

TTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGG  >  NC_000913/3950205‑3950498
                                                                                                                                                |                                                                                                                                                     
ttATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgcc                                                                                                                                                 <  2:74786‑M1/151‑7 (MQ=255)
                                             atCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaa                                                                                                    >  1:156389‑M1/1‑100 (MQ=255)
                                                    aCTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggca                                                                                              >  1:198677‑M1/1‑93 (MQ=255)
                                                       aTTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcg                                                                                          >  1:21615‑M1/1‑90 (MQ=255)
                                                       aTTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcg                                                                                          >  2:268658‑M1/1‑90 (MQ=255)
                                                         ttACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcga                                                                                         >  1:139867‑M1/1‑88 (MQ=255)
                                                                     tGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaac                                                                            <  1:68992‑M1/151‑76 (MQ=255)
                                                                     tGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaac                                                                            <  2:117781‑M1/151‑76 (MQ=255)
                                                                     tGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaac                                                                            <  2:175871‑M1/151‑76 (MQ=255)
                                                                                   gCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCTagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaacagtgatgggcca                                                                <  1:51954‑M1/149‑88 (MQ=255)
                                                                                     aTTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaacagtgatgggccatacc                                                            <  2:170913‑M1/151‑92 (MQ=255)
                                                                                        ccTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaacagtgatgggccatacctg                                                          >  1:165634‑M1/1‑57 (MQ=255)
                                                                                        ccTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaacagtgatgggccatacctg                                                          >  2:261456‑M1/1‑57 (MQ=255)
                                                                                                     tGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaacagtgatgggccatacctgcttcatgtctcaat                                            <  2:83641‑M1/151‑108 (MQ=255)
                                                                                                                                               ccagcgccttcggcatccatggccaacacatcacccggaaagaccaggttgaagcggcactcgacaccatgctgaacagtgatgggccatacctgcttcatgtctcaatcgacgaacttgagaacgtctggccgctggtgccgcctggcgc  <  1:185403‑M1/151‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGG  >  NC_000913/3950205‑3950498

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.