Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913_3_vch_pgi 1,978,503 Δ776 bp insB1insA insB1, insA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_vch_pgi 1978503–1979269 1979278 10–776 9 [0] [0] 9 insB1–insA insB1,insA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi = 19785020 (0.000)7 (0.490) 7/274 0.8 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913_3_vch_pgi 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

GTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_vch_pgi/1978405‑1978650
                                                                                                 |                                                                                                                                                    
gTACTCGTCTCTACACGCAGATACAACGGAGTGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                   <  1:65745‑M1/149‑52 (MQ=255)
              acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                     <  2:239596‑M1/149‑66 (MQ=255)
                    aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                               >  2:199288‑M1/1‑78 (MQ=255)
                        agcggAGGGGGGCTGCGATTTTCAATAATGGGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtggatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                           <  2:252683‑M1/147‑76 (MQ=255)
                                          ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtaaaaaaatttttcgctacgctagcaagcaaaaatgaaaacggaaaacacgaaat                                                         >  2:13122‑M1/1‑56 (MQ=255)
                                                                                aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                   >  1:174649‑M1/1‑18 (MQ=255)
                                                                                         aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgc          <  2:180263‑M1/149‑141 (MQ=255)
                                                                                                tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgccacgctgccaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaat   <  1:262131‑M1/149‑148 (MQ=255)
                                                                                                 gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaata  >  1:200773‑M1/1‑1 (MQ=255)
                                                                                                 |                                                                                                                                                    
GTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_vch_pgi/1978405‑1978650

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ 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.