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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_bme_pgi 1978503 1979278 776 40 [0] [0] 41 insB1–insA insB1,insA

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

GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_bme_pgi/1979137‑1979427
                                                                                                                                              |                                                                                                                                                    
gggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                                                                                <  1:395222‑M2/7‑1 (MQ=255)
          acttatctaagatgtttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattactcaacataaatgTATAAGTCATACttttg                                                                                                                                      >  1:363325‑M2/133‑149 (MQ=255)
                ctaagatttttcctaaatcgacgcagctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatacgtgatgcagatcacacaaagcactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  1:754221‑M2/23‑1 (MQ=255)
                ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  1:289135‑M2/23‑1 (MQ=255)
                 taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTtggt                                                                                                                               >  1:225341‑M2/126‑148 (MQ=255)
                 taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                               >  1:590161‑M2/126‑149 (MQ=255)
                        tttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGCAttt                                                                                                                        <  2:637032‑M2/31‑1 (MQ=255)
                        tttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggatgggattttgaataatgagcgaggcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGGTTTGGGTGTAttt                                                                                                                        >  1:823392‑M2/119‑149 (MQ=255)
                            ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    >  2:203858‑M2/115‑149 (MQ=255)
                            ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    >  2:406046‑M2/115‑149 (MQ=255)
                                 tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                               >  1:500192‑M2/110‑149 (MQ=255)
                                                tcgtcactacacgcacatacaacggaggagggatgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                                >  2:825967‑M2/95‑149 (MQ=255)
                                                      ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                          <  1:751264‑M2/61‑1 (MQ=255)
                                                          acgcacatacaacggaggggggctgcgattttcaacaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                      <  1:365409‑M2/65‑1 (MQ=255)
                                                            gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                                                                    >  2:113150‑M2/83‑149 (MQ=255)
                                                             gacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTGGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                                                                   <  2:424952‑M2/68‑1 (MQ=255)
                                                               catacaacggaggggggctgcgattttaaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                          <  1:549606‑M2/61‑1 (MQ=255)
                                                               catacaacggaggggggctgcgattttaaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                          >  2:549606‑M2/80‑140 (MQ=255)
                                                                    aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                            <  2:459057‑M2/59‑1 (MQ=255)
                                                                    aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                            >  1:459057‑M2/75‑133 (MQ=255)
                                                                                ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                                <  1:203858‑M2/87‑1 (MQ=255)
                                                                                    gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                >  1:566185‑M2/59‑81 (MQ=255)
                                                                                    gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  2:566185‑M2/23‑1 (MQ=255)
                                                                                       tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAACTg                                                         >  1:173178‑M2/56‑149 (MQ=255)
                                                                                       tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg                                                         >  1:576423‑M2/56‑149 (MQ=255)
                                                                                            ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                                    >  1:432978‑M2/51‑149 (MQ=255)
                                                                                              aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAAACGAAATGGGATGtt                                                  >  1:494622‑M2/49‑149 (MQ=255)
                                                                                                 gcgtgatgcagatcacacaaaacactcaaatacttaacgtaaatgTATAACTCATACGTTTGTTTTGGGTGTATTTCAATCTGTTAGAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGCTgcg                                               >  1:707683‑M2/46‑149 (MQ=255)
                                                                                                   gtgatgcagatcacacaaaacactgaattacttaacataaatgTATAAGTCATACTTTTGTTTTAGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTACGca                                             >  1:725816‑M2/44‑149 (MQ=255)
                                                                                                        gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                         >  2:35553‑M2/39‑68 (MQ=255)
                                                                                                        gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                         <  1:35553‑M2/30‑1 (MQ=255)
                                                                                                        gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTc                                        >  1:514993‑M2/39‑149 (MQ=255)
                                                                                                            atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGCAGCATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                    >  1:539106‑M2/35‑149 (MQ=255)
                                                                                                                     aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGCTGTTGCTCACAGTCAAAATCACTCAcc                           >  2:337224‑M2/26‑149 (MQ=255)
                                                                                                                      acactcaatcacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCg                          <  1:561527‑M2/125‑1 (MQ=255)
                                                                                                                             attacttaacataaatgCATAACTCATACTTTTGTTTTGGGTGTATTTCAATCTGTCAAAAAGTTTTTCGATACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCACAATAACTCACAGTAaataa                   <  2:707683‑M2/132‑1 (MQ=255)
                                                                                                                                    aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg            <  2:514993‑M2/139‑1 (MQ=255)
                                                                                                                                         aaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                                                                        >  1:557840‑M2/6‑52 (MQ=255)
                                                                                                                                         aaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                                                                        <  2:557840‑M2/47‑1 (MQ=255)
                                                                                                                                           atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa     <  2:363325‑M2/146‑1 (MQ=255)
                                                                                                                                              tataAGTCATACTTTTGGTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  <  2:82405/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_bme_pgi/1979137‑1979427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: