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

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

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

CCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAG  >  NC_000913/1979212‑1979344
                                                                   |                                                                 
gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                   <  1:309583‑M2/1‑1 (MQ=255)
gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                   <  2:315824‑M2/1‑1 (MQ=255)
gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                   <  1:702273‑M2/1‑1 (MQ=255)
      tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                             <  1:172841‑M2/7‑1 (MQ=255)
            tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATATGTCATACt                                                       <  2:229063‑M2/13‑1 (MQ=255)
                 ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                  <  2:301561‑M2/18‑1 (MQ=255)
                   aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                >  2:642948‑M2/49‑68 (MQ=255)
                   aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                >  2:165304‑M2/49‑68 (MQ=255)
                        gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                           <  1:272471‑M2/25‑1 (MQ=255)
                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                       >  1:337871‑M2/40‑68 (MQ=255)
                            cgcagatcacacaaatcactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                       <  1:140903‑M2/29‑1 (MQ=255)
                             gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                      >  2:85439‑M2/39‑68 (MQ=255)
                                  tcgcacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:176832‑M2/34‑68 (MQ=255)
                                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                       >  2:128210‑M2/34‑62 (MQ=255)
                                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                       <  1:128210‑M2/29‑1 (MQ=255)
                                                  attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                 >  2:22981‑M2/18‑68 (MQ=255)
                                                            ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta       >  1:523922‑M2/8‑68 (MQ=255)
                                                            ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta       >  2:543314‑M2/8‑68 (MQ=255)
                                                                 tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  2:15278‑M2/3‑68 (MQ=255)
                                                                   |                                                                 
CCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAG  >  NC_000913/1979212‑1979344

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