Predicted mutation
evidence seq id position mutation annotation gene description
MC JC AM260479 3,399,144 Δ1,321 bp [leuS]h16_A3140 [leuS], h16_A3140

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ AM260479 3399144–3400451 3400464 14–1321 24 [0] [0] 24 [leuS]–h16_A3140 [leuS],h16_A3140

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AM260479 = 3399143NA (NA)23 (0.910) 20/198 0.3 NA noncoding (2/1313 nt) ISAE1 insertion sequence ISAE1
?AM260479 3400465 = NA (NA)intergenic (+32/+53) h16_A3140/dapB Transposase ISAE1/dihydrodipicolinate reductase

GGTCGGGACAGCAGTCAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAGCCTCTTTTGAATTTCAAGTGGGTGGTATCGTGGCGGGTTTCAGGTAAAGAGCAACAAGAGCGATGCATAGCAAACTGTTTGAAGCAGCCCTTGGCATAAGCGATCCGTGGTTCGTCCGGGAAGTCGACT  >  AM260479/3399015‑3399273
                                                                                                                                |                                                                                                                                  
ggTCGGGACAGCAGTCAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgat                                                                                                                                >  1:1158586‑M1/1‑129 (MQ=255)
      gACAGCAGTCAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataacc                                                                                                                         >  2:128320‑M1/1‑123 (MQ=255)
        cagcagTCAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgt                                                                                                                       <  1:296190‑M1/134‑14 (MQ=255)
         agcagTCAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtg                                                                                                                      <  2:1158586‑M1/134‑15 (MQ=255)
                aGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcag                                                                                                               >  2:77303‑M1/1‑113 (MQ=255)
                    cAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggc                                                                                                             >  2:388939‑M1/1‑109 (MQ=255)
                            gcgcGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggc                                                                                                   >  1:352174‑M1/1‑101 (MQ=255)
                                        gCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgcccc                                                                                       <  2:42700‑M1/134‑46 (MQ=255)
                                        gCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgcccc                                                                                       <  2:352174‑M1/134‑46 (MQ=255)
                                                 cATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggca                                                                                                                >  2:162654‑M1/1‑80 (MQ=255)
                                                 cATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggca                                                                                                                <  1:162654‑M1/100‑21 (MQ=255)
                                                           aaTCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggc                                                                    >  2:235654‑M1/1‑70 (MQ=255)
                                                            aTCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcc                                                                   <  1:77303‑M1/134‑66 (MQ=255)
                                                               gTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggccca                                                                 <  1:388939‑M1/133‑68 (MQ=255)
                                                                   gAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacg                                                            <  1:235654‑M1/134‑73 (MQ=255)
                                                                   gAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacg                                                            <  2:947773‑M1/134‑73 (MQ=255)
                                                                                             gcgcGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcg                                                <  1:608981‑M1/120‑85 (MQ=255)
                                                                                             gcgcGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcg                                                >  2:608981‑M1/1‑36 (MQ=255)
                                                                                                cggcgCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggccgttgggctt                               >  2:427718‑M1/1‑33 (MQ=255)
                                                                                                     gccgccggGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggccgttgggcttgtcgg                          >  2:730918‑M1/2‑28 (MQ=255)
                                                                                                      ccgccgGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgcctta                                                                            <  2:1000579‑M1/83‑57 (MQ=255)
                                                                                                      ccgccgGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgcctga                                                                            >  1:1000579‑M1/1‑27 (MQ=255)
                                                                                                                    aTGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggccgttgggcttgtcggcgaggaaacgggcgg           <  1:730918‑M1/134‑122 (MQ=255)
                                                                                                                              agacgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggccgttgggcttgtcggcgaggaaacgggcggcgcgcacgg  >  1:383711‑M1/1‑3 (MQ=255)
                                                                                                                                |                                                                                                                                  
GGTCGGGACAGCAGTCAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAGCCTCTTTTGAATTTCAAGTGGGTGGTATCGTGGCGGGTTTCAGGTAAAGAGCAACAAGAGCGATGCATAGCAAACTGTTTGAAGCAGCCCTTGGCATAAGCGATCCGTGGTTCGTCCGGGAAGTCGACT  >  AM260479/3399015‑3399273

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