Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 1,978,503 Δ776 bp 100% insB5insA5 insB5, insA5

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979270 1979278 9–776 20 [0] [1] 21 insB5–insA5 insB5,insA5

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19785020 (0.000)19 (1.430) 17/384 NT 100% intergenic (‑305/+16) flhD/insB5 DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB
?NC_000913 1979279 = 0 (0.000)intergenic (‑64/‑474) insA5/uspC IS1 family protein InsA/universal stress protein C

CATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTT  >  NC_000913/1979029‑1979536
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                 
tgggccccgggaagccccggtaaaaaataattagcattagaagaattgcgataagctgcaataagcagaaccacctttttggttgaatatgtccttacaaagagaaatgggtctttacacttatctaagatttttcctaattcgacgcaactgtactcgtcactacacgcacatacaacggaggtgggctgcgattttcaataatgcgtgatgcagatcacataaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                                                                                                                                           <  2:168248‑M2/25‑1 (MQ=255)
                                                       ctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                                                                                                                                                                                 >  2:141124‑M2/196‑214 (MQ=255)
                                                           aataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                                                                                                                             <  1:141124‑M2/23‑1 (MQ=255)
                                                                                                acaaatagaaatgggtctttacacttatctcagatttttcctaaatcgacgcaactgtactagtccctacacgcaaatacaccggcggggggctgcgatcttcaataatgcgtgctgcagatcacacaaaacacccaattacttaccctacatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAACAAGTTTTTcgc                                                                                                                                                                                                              >  1:326‑M2/155‑208 (MQ=255)
                                                                                                                                                                cactacacgcacataccacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATActtttgtttg                                                                                                                                                                                                                                               >  1:130100‑M2/91‑110 (MQ=255)
                                                                                                                                                                        gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAATTCATACTTTTTTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGAAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACttt                                                                                                        >  1:31378‑M2/83‑238 (MQ=255)
                                                                                                                                                                        gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGCAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAAAAATCATCTGCTATAAATAATCACttt                                                                                                        >  1:19330‑M2/83‑238 (MQ=255)
                                                                                                                                                                          acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTCGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATa                                                                                                                                               >  1:46258‑M2/81‑197 (MQ=255)
                                                                                                                                                                             tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACttt                                                                                                        >  1:98730‑M2/78‑233 (MQ=255)
                                                                                                                                                                                 acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATg                                                                                                    <  2:98730‑M2/160‑1 (MQ=255)
                                                                                                                                                                                     ggggggtttttatttttaaatattgcgtaatacagatcacaaaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGGTTTTCGCTACGCTAGCAAGCAAAAATGACACAGGAATAATCGAAATGGGATGTTGCGAACAGTCAAAAGAACTCACCGTAAATAATCATCTGCTATAAACAATCACTTTCATGCACGACCAGCTAAGCTATTTTTNCACAGAAACTTACCGCAcaanaa                                                      <  2:261880‑M2/206‑1 (MQ=255)
                                                                                                                                                                                      ggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                                                                                                                                                     >  1:224066‑M2/69‑163 (MQ=255)
                                                                                                                                                                                          ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGCTGTTGCGCACAGTCaaaa                                                                                                                                                 >  1:49176‑M2/65‑179 (MQ=255)
                                                                                                                                                                                          ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGAGCACAGTCAAAATAACTCACCTTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGcaca                                                          >  1:261880‑M2/65‑266 (MQ=255)
                                                                                                                                                                                          ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATCTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                                                                                                                                                 <  2:224066‑M2/99‑1 (MQ=255)
                                                                                                                                                                                              gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAc                                                                                                                                             <  2:49176‑M2/119‑1 (MQ=255)
                                                                                                                                                                                                                             acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTAtttt                                                                               >  1:201827‑M2/30‑210 (MQ=255)
                                                                                                                                                                                                                                 aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAACAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTaaa                                                                           <  2:201827‑M2/185‑1 (MQ=255)
                                                                                                                                                                                                                                    actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACAGTAAATAATCATCTGCTATAAATAATCACTTTCATGCACTACCAGATAAGCTATTTTTAAACAGACACTTACCg                                                              >  1:75554‑M2/23‑220 (MQ=255)
                                                                                                                                                                                                                                         atttcataaaataaatttagaATTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGGTACGCGAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATACTCAACTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgccatctgg  <  1:124821‑M1/265‑48 (MQ=255)
                                                                                                                                                                                                                                                       atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACttt                                                                                                        >  1:169316‑M2/4‑159 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                 
CATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTT  >  NC_000913/1979029‑1979536

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: