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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979269 1979278 10–776 36 [0] [0] 36 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)32 (0.660) 27/260 NT 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

ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCT  >  NC_000913/1979135‑1979419
                                                                                                                                                |                                                                                                                                            
atgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgtgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                          <  2:757281‑M2/5‑1 (MQ=255)
              ttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgctcatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacctaaatgTATAAGTCATACttttgtt                                                                                                                            >  2:103996‑M2/131‑149 (MQ=255)
                atctaagatttttcctaaatcgacgcaactgtagtcgtcactacacgcacatacaacggaagggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattcattaacatgaatgTATAAGTCATACtttattttt                                                                                                                          >  2:314143‑M2/129‑149 (MQ=255)
                  ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgcgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                        <  2:798132‑M2/23‑1 (MQ=255)
                    aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                      >  1:215190‑M2/125‑149 (MQ=255)
                     agatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                                                                     >  2:89435‑M2/124‑149 (MQ=255)
                            tcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggtgggcggcgtcgagtttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                              <  1:623873‑M2/33‑1 (MQ=255)
                               taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTtcaagc                                                                                                           >  2:188544‑M2/114‑147 (MQ=255)
                                aaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatccagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTTGGTGTATTTTAATCt                                                                                                          >  2:476390‑M2/113‑149 (MQ=255)
                                 aatcgacgcaaccgtacccgtcaccacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                         <  1:184636‑M2/38‑1 (MQ=255)
                                              gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                            >  1:289921‑M2/99‑149 (MQ=255)
                                                  tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTtggca                                                                                        >  2:637062‑M2/95‑148 (MQ=255)
                                                        ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                  >  2:523905‑M2/89‑149 (MQ=255)
                                                          acacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                >  1:534145‑M2/87‑149 (MQ=255)
                                                                        cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                  <  1:530921‑M2/77‑1 (MQ=255)
                                                                              ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                                                            >  2:315735‑M2/67‑149 (MQ=255)
                                                                              ggggctgcgattttcaataatacgtgatgcagagcacacaaaacactcaattacttaacataaatgTATAAGTGATACTTTAGTTTTGGGTGTATTTCAATCTGTTGAAAAGTTTTATGCTACGCTAGCAAGAaaaaaggaaaaaggat                                                            >  2:512013‑M2/67‑148 (MQ=255)
                                                                                ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                          <  2:115512‑M2/85‑1 (MQ=255)
                                                                                ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                          >  2:118187‑M2/65‑149 (MQ=255)
                                                                                     cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa                                                     >  1:620482‑M2/60‑149 (MQ=255)
                                                                                       attttcaataatgcctgatgcagatcacacaaaacactcaattacttaagataaatgTATAAGTCATACTTTTGTATTGGGTGTATTTCAATCCGTTAAAACGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                   >  1:673297‑M2/58‑149 (MQ=255)
                                                                                       attttcaataatgcatgatgcagatcacacaaaacactcaattacttaaaataaatgTATAAGTCATACTTTTGTGTTGGGTGTATTTCAATCTGTTAAAACGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                   >  1:673264‑M2/58‑149 (MQ=255)
                                                                                          ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTCTAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCACTACGCGAGCAATCAAAAATGGAACAGGAATGATcgatatag                                                >  2:129697‑M2/55‑147 (MQ=255)
                                                                                             aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATGGGACAGACTGCTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                             <  1:664260‑M2/98‑1 (MQ=255)
                                                                                                     gtgatgcagatcacacaaaacacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                     <  2:253920‑M2/106‑1 (MQ=255)
                                                                                                     gtgatgcagatcacaaaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                     <  1:597272‑M2/106‑1 (MQ=255)
                                                                                                      tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTCTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATTTTGCGcac                                    >  2:592260‑M2/43‑149 (MQ=255)
                                                                                                         tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                 >  1:767552‑M2/40‑149 (MQ=255)
                                                                                                         tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                 >  1:720264‑M2/40‑149 (MQ=255)
                                                                                                           cagatcacacaaaacactccattccttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                               <  1:197187‑M2/112‑1 (MQ=255)
                                                                                                              atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                            <  2:363352‑M2/115‑1 (MQ=255)
                                                                                                              atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                            >  2:488366‑M2/35‑149 (MQ=255)
                                                                                                              atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                            >  1:753176‑M2/35‑149 (MQ=255)
                                                                                                                                   cttaacataaatgTATAAGTCATACTTTTCTTTTGGGTGTATTTCAATCTGTTTAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAATCAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAATTCACCGTAAATAatcat       <  2:673264‑M2/136‑1 (MQ=255)
                                                                                                                                   cttaacataaatgTATAAGTCATACTTTTCTTTTGGGTGTATTTCAATCTGTTCAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAGCAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat       <  2:673297‑M2/136‑1 (MQ=255)
                                                                                                                                        aataaatgTATAGGTCCTACTTTTGTTTTGGGTGTATTTCAATGTGTTAAAAAGTTTGTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCt  <  2:656817‑M2/141‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCT  >  NC_000913/1979135‑1979419

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