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 1979278 776 26 [0] [4] 30 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)26 (0.540) 24/396 1.0 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

GTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCA  >  NC_000913/1978297‑1978726
                                                                                                                                                                                                             |                                                                                                                                                                                                                                
gatgtgataagCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGTGTCTTTACACTTATCTAAGGTATTGCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACACCGGAGGGCGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                                                                                                             <  1:565299‑M1/222‑22 (MQ=255)
      aTAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                                                                                               >  2:745562‑M1/1‑200 (MQ=255)
      aTAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                                                                                               <  1:745562‑M1/203‑4 (MQ=255)
             gCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTGTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGTGCTGCGATTTTCAATAATGCGTGATGCAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                                                                                                    <  1:817493‑M1/287‑95 (MQ=255)
              cAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAAAGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                                                                                                  <  1:302305‑M1/288‑97 (MQ=255)
                        aaCCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAATATTCTTCCTAAATCGACGCAACTGTACTCGTCCCTACACGCCCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATTCAGAACACACCAAACACTCAATTACTTACCATAAATGtataagtcatacttttgttttggttgaatttcaatctggtaaaaagttttttgctactctatcaagcaaaaatacaaaagacataantaaaacggganntttttnt                                                                                                                        >  2:562058‑M1/1‑182 (MQ=255)
                                               tCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                                                                                                 >  1:837318‑M1/1‑159 (MQ=255)
                                                    aCAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATCAATGtataagtcatacttttgttttgcgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggaggttgcgcacagtcagaataactcaccgtaaataatc                                                                                            >  1:773218‑M1/1‑154 (MQ=255)
                                                                    cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaaaaatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataat                                                                            <  1:505698‑M1/288‑151 (MQ=255)
                                                                                                       aaacGGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCAATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatttgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaa                                         <  1:784203‑M1/287‑186 (MQ=255)
                                                                                                         aCTGTACTCGTCCCTACACGCACATACAACGGGGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATGACTTAACATAAATGgataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcacgcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg                                                                                       <  1:746916‑M1/240‑140 (MQ=255)
                                                                                                                  gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataa                                                                                 >  1:416896‑M1/1‑92 (MQ=255)
                                                                                                                     aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                                                                                                           >  1:736931‑M1/1‑89 (MQ=255)
                                                                                                                       tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                                         <  1:671886‑M1/112‑26 (MQ=255)
                                                                                                                       tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                                         >  2:671886‑M1/1‑87 (MQ=255)
                                                                                                                              acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                                                                                                  <  1:50804‑M1/176‑97 (MQ=255)
                                                                                                                              acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                                                                                                  >  2:50804‑M1/1‑80 (MQ=255)
                                                                                                                                      cGGAGGGGGGCTGCGATTTTCAATAAAGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatccttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagccagcaaaaatgaaacaggaatactcgaaatggcatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataattaatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaaggaag          >  1:25725‑M1/1‑72 (MQ=255)
                                                                                                                                        gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                                                                                                    <  2:759659‑M1/180‑111 (MQ=255)
                                                                                                                                        gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                                                                                                    >  1:759659‑M1/1‑70 (MQ=255)
                                                                                                                                              ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaaggaaggtgcgaac  >  1:351043‑M1/1‑64 (MQ=255)
                                                                                                                                                                   gTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttacaaagtttttcgctccgctaacaagcaaaaatgaaacaggaataatcgaaatg                                                                                                                                    <  1:32603‑M1/137‑95 (MQ=255)
                                                                                                                                                                            aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcacaaatgaaacaggaataatcgaaatgggatgttgc                                                                                                                           >  2:32603‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                  ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacaatcaaaataactcacagtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaag              >  1:803933‑M1/1‑12 (MQ=255)
                                                                                                                                                                                                    aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                                                                                                       >  1:705859‑M1/1‑10 (MQ=255)
                                                                                                                                                                                                    aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                                                                                                       <  2:705859‑M1/101‑92 (MQ=255)
                                                                                                                                                                                                             |                                                                                                                                                                                                                                
GTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCA  >  NC_000913/1978297‑1978726

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