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

CGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGCGC  >  NC_000913/3951877‑3952170
                                                                                                                                                 |                                                                                                                                                    
ttatattgaaattattaaacgcatcataaaaatcggccaaaaaatatcttgtactatttacaaaacctatggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGcc                                                                                                                                                 <  2:74786‑M2/6‑1 (MQ=255)
                                             atcttgtactatttacaaaacctatggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGaa                                                                                                    >  1:156389‑M2/101‑151 (MQ=255)
                                                    actatttacaaaacctatggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCa                                                                                              >  1:198677‑M2/94‑150 (MQ=255)
                                                       atttacaaaacctatggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCg                                                                                          >  1:21615‑M2/91‑151 (MQ=255)
                                                       atttacaaaacctatggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCg                                                                                          >  2:268658‑M2/91‑151 (MQ=255)
                                                         ttacaaaacctatggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGa                                                                                         >  1:139867‑M2/89‑150 (MQ=255)
                                                                     tggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAAc                                                                            <  1:68992‑M2/75‑1 (MQ=255)
                                                                     tggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAAc                                                                            <  2:117781‑M2/75‑1 (MQ=255)
                                                                     tggtaactctttaggcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAAc                                                                            <  2:175871‑M2/75‑1 (MQ=255)
                                                                                   gcattccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagctAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCa                                                                <  1:51954‑M2/87‑1 (MQ=255)
                                                                                     attccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATAcc                                                            <  2:170913‑M2/91‑1 (MQ=255)
                                                                                        ccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTg                                                          >  1:165634‑M2/58‑150 (MQ=255)
                                                                                        ccttcgaacaagatgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTg                                                          >  2:261456‑M2/58‑150 (MQ=255)
                                                                                                     tgcaagaaaagacaaaatgacagcccttctacgagtgattagccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAAt                                            <  2:83641‑M2/107‑1 (MQ=255)
                                                                                                                                               ccAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGgcgc  <  1:185403‑M2/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGCGC  >  NC_000913/3951877‑3952170

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.