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] 40 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)34 (0.910) 31/270 0.3 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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTC  >  NC_000913_3_bme_pgi/1978357‑1978646
                                                                                                                                                 |                                                                                                                                                
aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                               >  2:188185‑M1/1‑146 (MQ=255)
          ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGAGGCTGCGATTTTCAATAATGCGTGATGCAGATGACACAAAACACTCAATTACTTAACATAAAAGtataagtcatact                                                                                                                                     >  1:221837‑M1/1‑136 (MQ=255)
              acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                                                                 >  2:288451‑M1/1‑132 (MQ=255)
                         aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttcttttgggtgta                                                                                                                      <  1:425400‑M1/149‑29 (MQ=255)
                              tCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttggtttgggtgtatttca                                                                                                                 >  1:234131‑M1/1‑116 (MQ=255)
                                cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaaa                                                                                                               >  2:579377‑M1/1‑114 (MQ=255)
                                 tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                              >  1:420200‑M1/1‑113 (MQ=255)
                                  aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                             <  1:288033‑M1/149‑38 (MQ=255)
                                        aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                       >  2:90435‑M1/1‑106 (MQ=255)
                                                gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                               <  1:688189‑M1/149‑52 (MQ=255)
                                                  aCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                             <  2:356118‑M1/149‑54 (MQ=255)
                                                     cgccaccacacgCACATACAACGGAGGGGGGCTGCGATTTTCAATAACGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAGTGtataagtcacacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                          <  1:98820‑M1/146‑57 (MQ=255)
                                                         aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttctgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                                                                      <  1:546041‑M1/149‑61 (MQ=255)
                                                          cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAAGAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                     <  1:150762‑M1/149‑62 (MQ=255)
                                                             cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                                                                  <  1:449070‑M1/149‑65 (MQ=255)
                                                                        aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                       >  1:146274‑M1/1‑74 (MQ=255)
                                                                        aaCGGAGGGCGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                       <  2:551050‑M1/149‑76 (MQ=255)
                                                                               gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                                <  1:90435‑M1/149‑83 (MQ=255)
                                                                                    cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                           >  1:193606‑M1/1‑62 (MQ=255)
                                                                                     tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttttttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                          >  2:48647‑M1/1‑61 (MQ=255)
                                                                                            ttCAATAAGGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                   <  1:403978‑M1/149‑96 (MQ=255)
                                                                                                    tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtgtcttgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                           <  2:422984‑M1/149‑104 (MQ=255)
                                                                                                         gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                      <  2:681971‑M1/149‑109 (MQ=255)
                                                                                                               gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttttttttttttttcaaccaaaaatggaacaggaataatcgaaacgggatgtggcgcacagtctag                                >  2:160105‑M1/1‑35 (MQ=255)
                                                                                                               gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                                >  2:665336‑M1/1‑35 (MQ=255)
                                                                                                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                              <  1:47329‑M1/149‑117 (MQ=255)
                                                                                                                   acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                            >  2:95167‑M1/1‑31 (MQ=255)
                                                                                                                   acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                  >  1:254205‑M1/1‑31 (MQ=255)
                                                                                                                   acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                  <  2:254205‑M1/111‑81 (MQ=255)
                                                                                                                     acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttaatttgggtgtattacaatctgttaaaaaatttttcgctacgctagcaagcaaaaatgaaacagaaataatcgaaatggaatgttacgcacagtcaaaataact                          >  2:123400‑M1/1‑29 (MQ=255)
                                                                                                                                 aTTACTTAACATAAATGtatcagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa              <  1:441685‑M1/149‑133 (MQ=255)
                                                                                                                                 aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa              <  1:441160‑M1/149‑133 (MQ=255)
                                                                                                                                   tACTTAACATAAATGtatacgtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaattaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc            <  2:386961‑M1/149‑135 (MQ=255)
                                                                                                                                        aaCAGAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaagtaactcaccgtaaataatcatctg       <  1:584535‑M1/149‑140 (MQ=255)
                                                                                                                                          cttaaATGtataagttatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct     <  2:649724‑M1/147‑142 (MQ=255)
                                                                                                                                          cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct     <  1:280653‑M1/149‑142 (MQ=255)
                                                                                                                                           aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta    >  1:586933‑M1/1‑7 (MQ=255)
                                                                                                                                            tacatgtataagtcagactgttgttttgggtgtatttcaacctgttaaaaagtttttcgctacgcgagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctat   <  2:367647‑M1/146‑144 (MQ=255)
                                                                                                                                            tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctat   <  1:188185‑M1/149‑144 (MQ=255)
                                                                                                                                             aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctata  <  1:67853‑M1/149‑145 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTC  >  NC_000913_3_bme_pgi/1978357‑1978646

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