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–1979270 1979278 9–776 51 [0] [0] 51 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)47 (0.880) 37/268 0.5 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

CGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979140‑1979421
                                                                                                                                           |                                                                                                                                              
tctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                             <  1:65066‑M2/4‑1 (MQ=255)
tctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                             >  2:616198‑M2/140‑143 (MQ=255)
 ctgtacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                            <  2:442384‑M2/5‑1 (MQ=255)
       gctgatccaagacttctgctaaagcgaggcaactgtactcgtcactacacgcacatacaacggaggggggctgcgcgttgcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                      <  2:323417‑M2/11‑1 (MQ=255)
       acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaccactcaattacttaacataaatgTATAAGTCATa                                                                                                                                      >  1:551330‑M2/133‑143 (MQ=255)
       acctatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaactgaggtgggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                      <  2:620204‑M2/11‑1 (MQ=255)
                  agttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggagggagggagggatgtgcaagaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                           <  1:467853‑M2/22‑1 (MQ=255)
                   tttttcctacatcgacgccactgtactcgtcactacacgcacatacagcggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  1:120884‑M2/23‑1 (MQ=255)
                   tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgccgatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  1:331549‑M2/23‑1 (MQ=255)
                         ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                                                                    >  1:397627‑M2/115‑143 (MQ=255)
                         ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTGTGTTATGGGTGTAt                                                                                                                    >  2:731673‑M2/115‑143 (MQ=255)
                          taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaacacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                   >  2:103875‑M2/114‑143 (MQ=255)
                                     aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                                                                        >  1:334658‑M2/103‑143 (MQ=255)
                                         gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                    <  1:564918‑M2/45‑1 (MQ=255)
                                             tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAGACTATTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                                >  2:522489‑M2/95‑143 (MQ=255)
                                               gtcactgcacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagaacacacattgccatcaatccctcaccataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                              <  1:469653‑M2/51‑1 (MQ=255)
                                                  actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTcgc                                                                                           >  1:14973‑M2/90‑143 (MQ=255)
                                                    tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgata                                                                                         >  1:206319‑M2/88‑143 (MQ=255)
                                                      cacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                                                                       <  1:451070‑M2/58‑1 (MQ=255)
                                                          cacatacaacgtaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                   <  2:369113‑M2/62‑1 (MQ=255)
                                                           acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                  <  1:766755‑M2/63‑1 (MQ=255)
                                                             atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTAGAAGTCAAACTTTTGTATTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                >  2:761077‑M2/79‑143 (MQ=255)
                                                                  acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                                                           >  1:901309‑M2/74‑143 (MQ=255)
                                                                    ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                         <  1:174461‑M2/72‑1 (MQ=255)
                                                                      aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                       >  1:334445‑M2/70‑143 (MQ=255)
                                                                           ggctgcgtgtgtcaataatgcgtgatgccgatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                  <  2:632798‑M2/79‑1 (MQ=255)
                                                                           ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAATTTTTCGCTACGCTAGCAAGCAAAAatcaacaa                                                                  >  2:594833‑M2/65‑140 (MQ=255)
                                                                           ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTAGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTAGCATGCAAAAATGAAACa                                                                  >  2:595247‑M2/65‑143 (MQ=255)
                                                                            gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                 <  2:222203‑M2/80‑1 (MQ=255)
                                                                            gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGAGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                 <  2:276697‑M2/80‑1 (MQ=255)
                                                                              tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                               <  2:274690‑M2/82‑1 (MQ=255)
                                                                                cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                             >  1:102544‑M2/60‑143 (MQ=255)
                                                                                        aatactccgtgatgcagctcacacaaaacacacaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                     <  1:350949‑M2/92‑1 (MQ=255)
                                                                                        aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                     <  1:62163‑M2/92‑1 (MQ=255)
                                                                                             tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                                <  1:616198‑M2/97‑1 (MQ=255)
                                                                                              gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATCCGTACTCTCCTAGTTTTGGGTGTATTTCTAACAGTTAAAAAGTTTTTCGCTACGTCAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                               >  1:470160‑M2/46‑143 (MQ=255)
                                                                                                gtgatgcagatcacacaaaaaactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                             >  2:811335‑M2/44‑143 (MQ=255)
                                                                                                        gatcacagaaatcgcgcaattagtgaacataaatgTATAAGTCATACTGTTGTTTTGGGTGTATTTCAATCTGTGAAAAAGGTTTTCGCTACGCGAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTGGCGcaca                                     <  2:644213‑M2/108‑1 (MQ=255)
                                                                                                            acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                 >  1:321708‑M2/32‑143 (MQ=255)
                                                                                                            acacaaaacactcaattacttaacataaatgTATAAGTCATACTCTTGTTTTGGGTATATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                 >  2:230139‑M2/32‑143 (MQ=255)
                                                                                                              acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                               >  2:339343‑M2/30‑143 (MQ=255)
                                                                                                              acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                               >  1:642268‑M2/30‑143 (MQ=255)
                                                                                                              acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                               >  1:30599‑M2/30‑143 (MQ=255)
                                                                                                              acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATAAAACAGGAATAATCGAAATAGTATGTTACGCACAGTcaac                               >  2:333671‑M2/30‑142 (MQ=255)
                                                                                                                aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATATCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAAACGAAATGGGAAGTTGACCACAATCAAAAt                             >  2:421362‑M2/28‑143 (MQ=255)
                                                                                                                  aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATTTTGCGCACAGTAAAAATaa                           >  2:352458‑M2/26‑143 (MQ=255)
                                                                                                                                  acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc           >  2:543162‑M2/10‑143 (MQ=255)
                                                                                                                                    ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTGTTTACTACTCTCGAATGCCAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat         >  2:214276‑M2/8‑143 (MQ=255)
                                                                                                                                          gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTa   >  1:736413‑M2/2‑143 (MQ=255)
                                                                                                                                          gTATAAGTAATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTa   <  2:743177‑M2/142‑1 (MQ=255)
                                                                                                                                           tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCtat  <  1:118448/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979140‑1979421

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: