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

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

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

TAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCT  >  NC_000913_3_pae_tpiA/1978354‑1978645
                                                                                                                                                    |                                                                                                                                               
tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                                                                >  1:23932‑M1/1‑149 (MQ=255)
       gggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                                                                         <  2:325981‑M1/150‑9 (MQ=255)
                 acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                               >  1:355450‑M1/1‑132 (MQ=255)
                                   cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                           >  2:585910‑M1/1‑114 (MQ=255)
                                   cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                           <  1:585910‑M1/136‑23 (MQ=255)
                                               aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                                 >  1:39531‑M1/1‑102 (MQ=255)
                                                aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                >  1:455793‑M1/1‑101 (MQ=255)
                                                       tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                         >  1:56196‑M1/1‑94 (MQ=255)
                                                         gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggggtatttcaatctgttaaaaagtttttcgctacg                                                                                       >  1:278688‑M1/1‑92 (MQ=255)
                                                          tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttccatctgttaaaaagtttttcgctacgc                                                                                      <  2:355450‑M1/150‑60 (MQ=255)
                                                          tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                      <  1:585703‑M1/150‑60 (MQ=255)
                                                                 acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                                                               >  2:89537‑M1/1‑84 (MQ=255)
                                                                       aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                         >  2:93178‑M1/1‑78 (MQ=255)
                                                                                       ttgCTATTTTAAATAATGAGTGATGGAGATCACACAAATCACTCAATTACTTACCATAAATGtataagtcacacttttgttttgggtgtatttctatctgttaaaaagtttttcgctacgctagcacgcaagaatgaaacaggaataatc                                                         <  1:547734‑M1/149‑89 (MQ=255)
                                                                                        tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                        >  1:363408‑M1/1‑61 (MQ=255)
                                                                                           gATTTTCAATAATGCGTGATGGAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                     <  2:270272‑M1/150‑93 (MQ=255)
                                                                                                 cAATCATGCGTGATGAAGATCACACAAAATACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtatttcaatctgttaaaaagtttttagctacggtagagagaaacaataaaacgagaataatcgaaatggtat                                               >  2:207082‑M1/1‑52 (MQ=255)
                                                                                                            gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacgg                                    >  1:245743‑M1/1‑41 (MQ=255)
                                                                                                                  gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                                  >  1:421226‑M1/1‑35 (MQ=255)
                                                                                                                  gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                                  <  2:421226‑M1/114‑80 (MQ=255)
                                                                                                                  aatCACACAAAACACGAAATTAAATAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagaaagaaaaaaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                              <  2:460758‑M1/149‑116 (MQ=255)
                                                                                                                   aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                 <  2:118493‑M1/114‑81 (MQ=255)
                                                                                                                   aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                 >  1:118493‑M1/1‑34 (MQ=255)
                                                                                                                          aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggctgttgcgcacagtcaaaataactcac                      >  2:397855‑M1/1‑27 (MQ=255)
                                                                                                                                 tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacaggcaaaataactcaccgtaaaa               >  2:297110‑M1/1‑20 (MQ=255)
                                                                                                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat         <  2:363408‑M1/150‑137 (MQ=255)
                                                                                                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcaacgctagcaagcaaaactcaaccatccagaatcgaaattgtatgttgcgcacagtcaaaataactcaccgtaaataatcatct       >  2:278821‑M1/1‑12 (MQ=255)
                                                                                                                                         ttAACATAAATGtataagtcatacttttgttttggcggaactacaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatagaaaggggatgttgagcacaggcaaaataactgaacgttattaatcatct       >  1:144198‑M1/1‑12 (MQ=255)
                                                                                                                                              aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcacacagtcaaaataactcaccgtgaataatcatctgctat  >  2:219717‑M1/1‑7 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCT  >  NC_000913_3_pae_tpiA/1978354‑1978645

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: 
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.