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

CGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAAT  >  NC_000913_3_pae_tpiA/1979130‑1979428
                                                                                                                                                     |                                                                                                                                                     
tagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                                                                                                       >  1:23932‑M2/150‑150 (MQ=255)
       gggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                                                                                <  2:325981‑M2/8‑1 (MQ=255)
                 acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                                      >  1:355450‑M2/133‑150 (MQ=255)
                                   ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                                  >  2:585910‑M2/115‑136 (MQ=255)
                                   ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                                  <  1:585910‑M2/22‑1 (MQ=255)
                                               aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                                                                        >  1:39531‑M2/103‑150 (MQ=255)
                                                actgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                                       >  1:455793‑M2/102‑150 (MQ=255)
                                                       tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                                >  1:56196‑M2/95‑150 (MQ=255)
                                                         gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGGGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                                                                              >  1:278688‑M2/93‑150 (MQ=255)
                                                          tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCCATCTGTTAAAAAGTTTTTcgctacgc                                                                                             <  2:355450‑M2/59‑1 (MQ=255)
                                                          tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                             <  1:585703‑M2/59‑1 (MQ=255)
                                                                 acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                                      >  2:89537‑M2/85‑150 (MQ=255)
                                                                       atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                                >  2:93178‑M2/79‑150 (MQ=255)
                                                                                       ttgctattttaaataatgagtgatggagatcacacaaatcactcaattacttaccataaatgTATAAGTCACACTTTTGTTTTGGGTGTATTTCTATCTGTTAAAAAGTTTTTCGCTACGCTAGCACGCAAGAATGAAACAGGAATAATc                                                                <  1:547734‑M2/88‑1 (MQ=255)
                                                                                        tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                               >  1:363408‑M2/62‑150 (MQ=255)
                                                                                           gattttcaataatgcgtgatggagatcacacaaaacactaaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                            <  2:270272‑M2/92‑1 (MQ=255)
                                                                                                 caatcatgcgtgatgaagatcacacaaaatactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGTTGTATTTCAATCTGTTAAAAAGTTTTTAGCTACGGTAGAGAGAAACAATAAAACGAGAATAATCGAAATGGTAt                                                      >  2:207082‑M2/53‑150 (MQ=255)
                                                                                                            gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGcgcacgg                                           >  1:245743‑M2/42‑148 (MQ=255)
                                                                                                                  gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                         >  1:421226‑M2/36‑114 (MQ=255)
                                                                                                                  gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                         <  2:421226‑M2/79‑1 (MQ=255)
                                                                                                                  aatcacacaaaacacgaaattaaataacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGAAAGAAAAAAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                     <  2:460758‑M2/115‑1 (MQ=255)
                                                                                                                   atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                        <  2:118493‑M2/80‑1 (MQ=255)
                                                                                                                   atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                        >  1:118493‑M2/35‑114 (MQ=255)
                                                                                                                          aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGTTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGCTGTTGCGCACAGTCAAAATAACTCAc                             >  2:397855‑M2/28‑150 (MQ=255)
                                                                                                                                 tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGGCAAAATAACTCACCGtaaaa                      >  2:297110‑M2/21‑149 (MQ=255)
                                                                                                                                       acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat                <  2:363408‑M2/136‑1 (MQ=255)
                                                                                                                                         ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCAACGCTAGCAAGCAAAACTCAACCATCCAGAATCGAAATTGTATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCt              >  2:278821‑M2/13‑150 (MQ=255)
                                                                                                                                         ttaacataaatgTATAAGTCATACTTTTGTTTTGGCGGAACTACAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATAGAAAGGGGATGTTGAGCACAGGCAAAATAACTGAACGTTATTAATCATCt              >  1:144198‑M2/13‑150 (MQ=255)
                                                                                                                                              ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCACACAGTCAAAATAACTCACCGTGAATAATCATCTGCtat         >  2:219717‑M2/8‑150 (MQ=255)
                                                                                                                                                     tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAaataat  >  1:228110/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAAT  >  NC_000913_3_pae_tpiA/1979130‑1979428

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.