Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 1,299,499 Δ1,199 bp 100% insH21 insH21

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1299499–1300697 1300697 1–1199 26 [0] [0] 27 insH21 insH21

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12994980 (0.000)20 (1.050) 15/376 NT 100% intergenic (+253/‑69) ychE/insH21 MarC family putative inner membrane protein YchE/IS5 family transposase and trans‑activator
?NC_000913 1300698 = 0 (0.000)intergenic (+151/‑484) insH21/oppA IS5 family transposase and trans‑activator/oligopeptide ABC transporter periplasmic binding protein

CACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGC  >  NC_000913/1299398‑1299723
                                                                                                    |                                                                                                                                                                                                                                 
cACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaaga                                                                                                                                                                        >  2:383216‑M1/1‑101 (MQ=255)
    atcatGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagataga                                                                                                                                                                    <  1:383217‑M1/160‑64 (MQ=255)
                                cTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgttt                                                                                                                          >  2:270226‑M1/1‑69 (MQ=255)
                                   ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggaac                   >  1:358773‑M1/1‑66 (MQ=255)
                                   ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagc                                                                                                      >  2:102227‑M1/1‑66 (MQ=255)
                                   ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaa                                                                                                                                                                            >  1:12321‑M1/1‑66 (MQ=255)
                                    gCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgtta                                                                                                                      <  1:270227‑M1/174‑110 (MQ=255)
                                       tGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctata                                                                                                  <  1:102227‑M1/191‑130 (MQ=255)
                                       tGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaaga                                                                                                                                                                        <  2:12321‑M1/121‑60 (MQ=255)
                                       tGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcaaaaaagatagacatgcttagccaatctctatttgattgaattgaaagat                                                                                                                              <  1:15551‑M1/163‑102 (MQ=255)
                                        gCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaacgatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaaaaaccgcatttggttatttt               >  2:313181‑M1/1‑61 (MQ=255)
                                        gCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgcat                                                                                                                                      >  2:250949‑M1/1‑61 (MQ=255)
                                        gCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagataga                                                                                                                                                                    >  1:35427‑M1/1‑61 (MQ=255)
                                        gCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgt                                                                                                                                                                                                                    >  2:175485‑M1/1‑61 (MQ=255)
                                            aagaaTAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaa                                                                                                                                  <  1:250950‑M1/154‑98 (MQ=255)
                                            aagaaTAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatg                                                                                                                                                                <  2:35427‑M1/124‑68 (MQ=255)
                                            aagaaTAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaa                                                                                                                                                                                                                <  1:175486‑M1/76‑20 (MQ=255)
                                            aagaaTAATTGAAATGATATTATTAATTCCACGTCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattcagatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccggg           <  1:313182‑M1/273‑217 (MQ=255)
                                                      gAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaatt                                                                                                                                     >  1:134666‑M1/1‑47 (MQ=255)
                                                                                        agGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaag                                                                                                       >  1:340468‑M1/1‑13 (MQ=255)
                                                                                            aaGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctat                                                                                                   <  2:340467‑M1/137‑129 (MQ=255)
                                                                                               tGCTAAataatattcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgt                                                                                                                            >  1:58293‑M1/1‑6 (MQ=255)
                                                                                               tGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatg                                                                                                               >  2:203602‑M1/1‑6 (MQ=255)
                                                                                                 cTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaaccgcca  >  2:210141‑M1/1‑4 (MQ=255)
                                                                                                   aaataatattcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgt                                                                                                                        <  2:58293‑M1/109‑108 (MQ=255)
                                                                                                   aaataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatg                                                                                                           <  1:203603‑M1/122‑121 (MQ=255)
                                                                                                    |                                                                                                                                                                                                                                 
CACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGC  >  NC_000913/1299398‑1299723

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