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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979270 1979278 9–776 41 [0] [3] 44 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)39 (1.070) 32/426 0.2 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

GCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCAT  >  NC_000913/1978243‑1978727
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 
gCACCCCGTGGTGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAAACGACGCAACTGTACTCGTCACTACACGCACATACATCGGAGGGAGGCTGGGATTTTCAATAATGCGGGATGCAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                                          <  1:401061‑M1/285‑26 (MQ=255)
        tGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCCACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAGTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                                          <  1:95923‑M1/277‑26 (MQ=255)
              cGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCCATAAGCAGAACAACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTCTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGTGTTTCAATAATGCGTGATGCAGACCACACAAAAAACACAATTCCTTAACATAAATGtttaagtaaaactttattttttggggggcgcccacacc                                                                                                                                                                                             >  2:95923‑M1/1‑246 (MQ=255)
                  gggAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACATTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGCCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGACGGCGATTTTCATTAATGTGTGATGCAGATCAAGACAAACACTCAATTACTTACAATAAATGtatatatcataaattttgtttggggggtgttttaattttttaa                                                                                                                                                                                        >  2:558242‑M1/1‑242 (MQ=255)
                                                       ttGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTCAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTCCACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtatatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                                                                                                   >  1:262523‑M1/1‑205 (MQ=255)
                                                                    cAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                                                                                                        <  1:562142‑M1/283‑92 (MQ=255)
                                                                     ncnccccccaaccacgttttTATTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGAATTTTCCTAGATGCACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAAGGCGTGATGGAGATCAAACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                                                                                                     <  2:171706‑M1/276‑95 (MQ=255)
                                                                                   ccTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTGATCTAAGATTTTTCTTAAATCTACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATATTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtatttcaatctgtgaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataataaaaatagtaggttccgcc                                                                                                                         >  1:171706‑M1/1‑177 (MQ=255)
                                                                                   ccTTTTTGGTTTAATAAGTACTTACAAATAGAAATTGGTCTTTACACTTATCTAAGATTTTTCCTAAAACGACGAAACTGTACGCGTCACTACACGCACATACAACGGAGGGGGGCTGCAATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttgtgttttgggtgtatttcaaactgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                                                                                                          <  1:558242‑M1/282‑106 (MQ=255)
                                                                                                          aCAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactcttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcacacagtcaaaataactcacagtaaata                                                                                                >  1:157825‑M1/1‑154 (MQ=255)
                                                                                                                                   tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                                                                                                              >  2:355413‑M1/1‑129 (MQ=255)
                                                                                                                                   tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                                                                                                              <  1:355413‑M1/182‑54 (MQ=255)
                                                                                                                                                             aacgaatccacaacactacaccgaaatACACCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCAAAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagc                                                    <  2:42964‑M1/265‑176 (MQ=255)
                                                                                                                                                                          cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaactggg                                                                                                                                   >  1:420085‑M1/1‑90 (MQ=255)
                                                                                                                                                                          cACTACACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                                                                                                   <  2:420085‑M1/186‑97 (MQ=255)
                                                                                                                                                                             tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttt                                              >  1:42964‑M1/1‑87 (MQ=255)
                                                                                                                                                                                    acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                                                                                                          <  1:409341‑M1/169‑90 (MQ=255)
                                                                                                                                                                                    acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                                                                                                          >  2:409341‑M1/1‑80 (MQ=255)
                                                                                                                                                                                      aTACAACGGAGGGGGTCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGgataagtcatacttttgttttgggtgtatttcaagctgttaaaaagttttgcgaaacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttacgcacaatcaaaataactcaccgtaaataatcat                                                                                           <  1:268780‑M1/214‑137 (MQ=255)
                                                                                                                                                                                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgaacaacaaac                   >  1:380861‑M1/1‑77 (MQ=255)
                                                                                                                                                                                            cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatg                                                                   >  2:453466‑M1/1‑72 (MQ=255)
                                                                                                                                                                                            cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatg                                                                   <  1:453466‑M1/232‑161 (MQ=255)
                                                                                                                                                                                            cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                                                                                                       <  2:128520‑M1/180‑109 (MQ=255)
                                                                                                                                                                                            cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                                                                                                       >  1:128520‑M1/1‑72 (MQ=255)
                                                                                                                                                                                            cGGAGGGGGGCTGCGATTTTCAAGAATGCATGATGAAGCTCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttccatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaagg              <  1:251861‑M1/285‑214 (MQ=255)
                                                                                                                                                                                              gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatatgctatcaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcaccacaaactaagga             >  1:436277‑M1/1‑70 (MQ=255)
                                                                                                                                                                                                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATAACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaattttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgaacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttacagcacaacaaactaagaaagt          >  1:222738‑M1/1‑68 (MQ=255)
                                                                                                                                                                                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                                                                                                            >  1:237784‑M1/1‑64 (MQ=255)
                                                                                                                                                                                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                                                                                                            <  2:237784‑M1/151‑88 (MQ=255)
                                                                                                                                                                                                         cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaaggaaggtgcgaaca  >  1:13812‑M1/1‑59 (MQ=255)
                                                                                                                                                                                                         cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaaggaaggtgcgaaca  >  1:582090‑M1/1‑59 (MQ=255)
                                                                                                                                                                                                                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                                                                                                   <  2:166986‑M1/125‑97 (MQ=255)
                                                                                                                                                                                                                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                                                                                                   >  1:166986‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                                                                                                >  2:428074‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                                                                                                <  1:428074‑M1/80‑52 (MQ=255)
                                                                                                                                                                                                                                           aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcacttt                                                                       <  1:65581‑M1/181‑157 (MQ=255)
                                                                                                                                                                                                                                           aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcacttt                                                                       >  2:65581‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                                                        ttAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                                          >  1:506780‑M1/1‑12 (MQ=255)
                                                                                                                                                                                                                                                        ttAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                                          <  2:506780‑M1/37‑26 (MQ=255)
                                                                                                                                                                                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgtacagt                                                                                                                     >  1:353510‑M1/1‑3 (MQ=255)
                                                                                                                                                                                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                                                                                                     <  2:353510‑M1/113‑111 (MQ=255)
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 
GCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCAT  >  NC_000913/1978243‑1978727

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