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

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

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

ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913/1979135‑1979427
                                                                                                                                                |                                                                                                                                                    
atgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                                  <  1:230419‑M2/5‑1 (MQ=255)
 tgggtctttactcttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatccaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                                 <  2:254161‑M2/6‑1 (MQ=255)
                  ctaagattttgcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacacttaagtacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  1:294631‑M2/23‑1 (MQ=255)
                                 aatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                                 <  2:141761‑M2/38‑1 (MQ=255)
                                  atcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAtctat                                                                                                                >  1:162410‑M2/111‑147 (MQ=255)
                                             tgtactcgtcacttcacgtacatacaactgaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATTCTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                                     <  2:373489‑M2/50‑1 (MQ=255)
                                              gtactcgtcactccacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                    <  1:656765‑M2/51‑1 (MQ=255)
                                                       actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                                                                           >  1:497178‑M2/90‑149 (MQ=255)
                                                         tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                         >  1:433619‑M2/88‑149 (MQ=255)
                                                                       acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                           >  1:620534‑M2/74‑149 (MQ=255)
                                                                          gatggtggctgcgattttcaattatgcgtgatgcagatcacacaaaacactcaattacgtaacataaatgTATAAGTCATACTTCTGTTTTGGGTGTATTTCAACCTGTCAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                        <  2:475809‑M2/79‑1 (MQ=255)
                                                                          gagggggtctgcgattttccagaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                        <  2:127306‑M2/79‑1 (MQ=255)
                                                                          gaggggggcggcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                        <  2:285709‑M2/79‑1 (MQ=255)
                                                                                  cagcgattttaaataatgcgtgatacagatcaaacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTCGGTGTATTTCAATCTGTTAAAAAGTTTTTCACTAAAATAAAAAACAAAAATGAAACAGGaataat                                                                <  1:6446‑M2/87‑1 (MQ=255)
                                                                                                 atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGACATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGGAACAGGAATAATCAAAATGGGATGTTg                                                 >  2:48002‑M2/48‑149 (MQ=255)
                                                                                                            agatcacaaaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaa                                      <  1:310481‑M2/113‑1 (MQ=255)
                                                                                                                              aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                    >  1:141162‑M2/19‑149 (MQ=255)
                                                                                                                                                tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  >  2:183032/1‑149 (MQ=255)
                                                                                                                                                tataAGTCATACTTTTGCTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  <  1:406252/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913/1979135‑1979427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ 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.