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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503 1979278 776 34 [0] [1] 34 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)25 (0.740) 19/114 0.3 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

GGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGC  >  NC_000913/1978435‑1978569
                                                                   |                                                                   
ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAGACACTCAATTACTTAACATAAATG                                                                     >  2:1801‑M1/1‑68 (MQ=255)
 ggtggctgCGATTTTCTATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                    <  2:64499‑M1/65‑2 (MQ=255)
 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                    <  1:114621‑M1/68‑2 (MQ=255)
 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                    <  2:289067‑M1/68‑2 (MQ=255)
 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                    <  1:648732‑M1/68‑2 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                   >  1:1209477‑M1/1‑66 (MQ=255)
   gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                  >  1:540479‑M1/1‑65 (MQ=255)
          gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAATGtataagtcata                                                          <  2:333168‑M1/68‑12 (MQ=255)
             tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                        <  1:77707‑M1/68‑14 (MQ=255)
             tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                        <  2:1107823‑M1/68‑14 (MQ=255)
                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                 >  2:1345569‑M1/1‑48 (MQ=255)
                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                 >  1:1432734‑M1/1‑48 (MQ=255)
                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                             <  1:642600‑M1/56‑9 (MQ=255)
                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                             >  2:642596‑M1/1‑48 (MQ=255)
                      tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                               >  1:637923‑M1/1‑46 (MQ=255)
                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                          <  1:286174‑M1/68‑28 (MQ=255)
                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                        <  1:85797‑M1/68‑30 (MQ=255)
                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttctatct                                >  1:21556‑M1/1‑31 (MQ=255)
                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgaatttcaatctgt                              >  1:573273‑M1/1‑29 (MQ=255)
                                          aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                           >  2:180726‑M1/1‑26 (MQ=255)
                                                   aTTACTTAGCATAAATGtataagtcatacttttgttttgggt                                            >  2:792551‑M1/1‑17 (MQ=255)
                                                   aTTACTTAACATAAATGtataagtcatacttttgttttgggt                                            <  1:792556‑M1/42‑26 (MQ=255)
                                                    ttACTTAACATAAATGtataagtcaaacttttgttttgggtgtatttcaatctgttaaaaagttattc                 >  2:1373512‑M1/1‑16 (MQ=255)
                                                      aCTTAACATAAATGtataagtcatacttttgttttgg                                              >  2:972408‑M1/1‑14 (MQ=255)
                                                      aCTTAACATAAATGtataagtcatacttttgttttgg                                              <  1:972413‑M1/37‑24 (MQ=255)
                                                         tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac            <  2:1175572‑M1/68‑58 (MQ=255)
                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta        >  2:1172254‑M1/1‑7 (MQ=255)
                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag       >  2:420904‑M1/1‑6 (MQ=255)
                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttattcgctacgctag       <  1:596328‑M1/68‑63 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag   >  2:302412‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag   >  2:346512‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag   >  2:81258‑M1/1‑2 (MQ=255)
                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag   >  1:1057711‑M1/1‑2 (MQ=255)
                                                                   gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc  <  2:116422‑M1/68‑68 (MQ=255)
                                                                   |                                                                   
GGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGC  >  NC_000913/1978435‑1978569

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: