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 1979278 776 31 [0] [0] 32 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)29 (1.400) 26/274 0.0 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

GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_vch_pgi/1978356‑1978650
                                                                                                                                                  |                                                                                                                                                    
gAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCCCACAAAACACTCAATTACTTACCATAAATGta                                                                                                                                                    >  1:142714‑M1/1‑147 (MQ=255)
 aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACGCGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAAGACTCAATTACTTAACATAAATGtat                                                                                                                                                   >  1:382283‑M1/1‑146 (MQ=255)
  aacgggtctTTACACTTATCTAAGATTTTTCCTAAATCGACTTCACTGCACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                                                                  <  1:158628‑M1/146‑5 (MQ=255)
        tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                            <  2:151341‑M1/149‑11 (MQ=255)
           ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                                         >  1:239751‑M1/1‑136 (MQ=255)
                     cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatccttttgttttgg                                                                                                                               >  1:270669‑M1/1‑126 (MQ=255)
                     cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                               <  1:122031‑M1/149‑24 (MQ=255)
                       aaGATTTTTCCTAAATCGACCCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                             <  1:243343‑M1/149‑26 (MQ=255)
                                 cTAAATCGACGCAACTGTACTCGTCACTACACGCACATAAAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaat                                                                                                                   <  2:5762‑M1/149‑36 (MQ=255)
                                              aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                                      <  1:270626‑M1/149‑49 (MQ=255)
                                                  tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                  <  1:369928‑M1/149‑53 (MQ=255)
                                                     tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcaaacttttgttttgggtgtatttcaatctgttcaacagtttttcgct                                                                                               >  2:25355‑M1/1‑94 (MQ=255)
                                                       gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                             <  1:237318‑M1/149‑58 (MQ=255)
                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                                   >  1:159627‑M1/1‑82 (MQ=255)
                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatacttttgttttgggtgtttttcaatctgttaaaaattttttcgctactctttcacgc                                                                                   >  2:169568‑M1/1‑82 (MQ=255)
                                                                              aGGGGGGCTGCCATTTTCAATAATGCGTCATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgatacgctagcaagcaaaaatgaaacag                                                                      <  2:11840‑M1/149‑81 (MQ=255)
                                                                                  gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacatgaat                                                                  >  2:186393‑M1/1‑65 (MQ=255)
                                                                                   ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                                 <  2:208109‑M1/149‑86 (MQ=255)
                                                                                            tttCAAAAATGCGGTATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgtgacaggaataatcgaaatg                                                        <  2:327809‑M1/149‑95 (MQ=255)
                                                                                              tCAATAATGCGAGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                      >  1:243301‑M1/1‑53 (MQ=255)
                                                                                               cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                                     <  1:193127‑M1/149‑98 (MQ=255)
                                                                                                aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                                    <  2:312979‑M1/149‑99 (MQ=255)
                                                                                                          gatatagataaCATAAAACGCTCAATTACTTAACAAAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaaaaaaaaaaaaaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                          <  2:374496‑M1/144‑109 (MQ=255)
                                                                                                          gATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatcttttaaaaagttattcgctacgctagcaagcaaaaatgaaacaggaataagcgaaatgggatgttgcgcaca                                          <  1:331048‑M1/149‑109 (MQ=255)
                                                                                                                        aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                            <  1:165627‑M1/149‑123 (MQ=255)
                                                                                                                               tCAATTACTTAACATAAATGtataagttatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                     <  1:34161‑M1/149‑130 (MQ=255)
                                                                                                                               tCAATTACTTAACACAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttacgcacagtcaaaataactcaccgtaaa                     <  2:382283‑M1/149‑130 (MQ=255)
                                                                                                                                      cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat              <  1:274532‑M1/149‑137 (MQ=255)
                                                                                                                                       ttAACAGAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc             <  1:320623‑M1/149‑138 (MQ=255)
                                                                                                                                               aaTGtataagtcatactttggttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataa     <  2:270669‑M1/149‑146 (MQ=255)
                                                                                                                                                  gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaata  >  1:63024‑M1/1‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_vch_pgi/1978356‑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: 
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.