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

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

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

TCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_vch_pgi/1978364‑1978650
                                                                                                                                          |                                                                                                                                                    
tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATAACACAAAACACGCAATTACTTAACATAAATGtataagtcat                                                                                                                                            >  1:102150‑M1/1‑139 (MQ=255)
           aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACCATCAAGGACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                                 <  1:411077‑M1/149‑22 (MQ=255)
           aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATATATGtataagtcatacttttgtttt                                                                                                                                 <  1:411057‑M1/149‑22 (MQ=255)
             cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                               <  1:165334‑M1/149‑24 (MQ=255)
              tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                              >  1:419745‑M1/1‑125 (MQ=255)
              tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtgttggg                                                                                                                              >  1:419934‑M1/1‑125 (MQ=255)
                   tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtat                                                                                                                         >  1:419180‑M1/1‑120 (MQ=255)
                          tAAATCGACGAAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                  >  1:540607‑M1/1‑113 (MQ=255)
                           aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                                 >  1:429108‑M1/1‑112 (MQ=255)
                           aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTAACAATAGTGCGTGATGCAGATCACACAACACACTCAATTACTTAACATAATTGtataagtcatacttttgttgtgagtgtatttcaatcg                                                                                                                 >  1:78764‑M1/1‑112 (MQ=255)
                             aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                               <  2:24346‑M1/149‑40 (MQ=255)
                                   gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                         >  1:283181‑M1/1‑104 (MQ=255)
                                          tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                  <  1:426525‑M1/149‑53 (MQ=255)
                                          gacgtgtcactACACGCTCATACAACGGAGGGGGGTTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                  <  2:78764‑M1/144‑53 (MQ=255)
                                               gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                             >  1:320646‑M1/1‑92 (MQ=255)
                                               gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                             <  2:320646‑M1/101‑10 (MQ=255)
                                                                 aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                           >  2:79652‑M1/1‑74 (MQ=255)
                                                                   cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                         <  1:120928‑M1/149‑78 (MQ=255)
                                                                   cGGAGGAGTGCTGCGATTTTCAATAATGCGTGATGCAGATCATACAAAACACTCAAGTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                         <  2:436872‑M1/149‑78 (MQ=255)
                                                                           ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                                 <  2:280741‑M1/149‑86 (MQ=255)
                                                                                   ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                         <  1:331030‑M1/149‑94 (MQ=255)
                                                                                        aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                                    >  1:206680‑M1/1‑51 (MQ=255)
                                                                                        aaaaatgcgTGATGCAGATCACACAAAACACTCAACTACTTAACATAAATGtataagtcatccttttggtttgggtgtatttcaatctattaaaaagtttttcgctacgctagcaaacaaacattaaacaaggacaatacaaatgggat                                                    >  2:484770‑M1/4‑51 (MQ=255)
                                                                                         agaatgcGTGATGCAGATGGCACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataattgaaattggatg                                                   >  2:467885‑M1/3‑50 (MQ=255)
                                                                                            atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                                <  2:419745‑M1/149‑103 (MQ=255)
                                                                                               cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctccgctagcaagcaaaaatgaaacaggaataagcggaatgtgatgttgcgc                                             >  1:172778‑M1/1‑44 (MQ=255)
                                                                                                  gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctatcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                          <  2:232242‑M1/149‑109 (MQ=255)
                                                                                                         aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                       <  1:525919‑M1/65‑32 (MQ=255)
                                                                                                         aTCACACAAAACACTCAATTACTTAACATAAACGtataagtcatacttttgttttgggtgtattt                                                                                                                       >  2:525919‑M1/1‑34 (MQ=255)
                                                                                                             cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                               >  2:294342‑M1/1‑30 (MQ=255)
                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                              <  2:416242‑M1/149‑121 (MQ=255)
                                                                                                                         aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                   >  2:527277‑M1/1‑18 (MQ=255)
                                                                                                                               ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc             >  2:462830‑M1/1‑12 (MQ=255)
                                                                                                                                    aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta        >  1:384560‑M1/1‑7 (MQ=255)
                                                                                                                                          gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaata  >  1:448583‑M1/1‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
TCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_vch_pgi/1978364‑1978650

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: