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 3400464 1321 22 [0] [0] 22 [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)21 (0.840) 21/196 0.2 NA noncoding (2/1313 nt) ISAE1 insertion sequence ISAE1
?AM260479 3400465 = NA (NA)intergenic (+32/+53) h16_A3140/dapB Transposase ISAE1/dihydrodipicolinate reductase

CAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAGCCTCTTTTGAATTTCAAGTGGGTGGTATCGTGGCGGGTTTCAGGTAAAGAGCAACAAGAGCGATGCATAGCAAACTGTTTGAAGCAGCCCTTGGCATAAGCGATCCGTGGTTCGTCCGGGAAGTCGACTTC  >  AM260479/3399030‑3399275
                                                                                                                 |                                                                                                                                    
cAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcgg                                                                                                                    >  2:652711‑M1/1‑114 (MQ=255)
 aGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggca                                                                                                                  <  2:418181‑M1/133‑21 (MQ=255)
                 ggCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagttt                                                                                                         >  2:940006‑M1/1‑97 (MQ=255)
                 ggCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagttt                                                                                                         <  1:940006‑M1/126‑30 (MQ=255)
                      cTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgc                                                                                            <  1:652711‑M1/134‑43 (MQ=255)
                              ccGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggc                                                                                    <  2:204255‑M1/134‑51 (MQ=255)
                                   aTCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgcctt                                                                               <  2:428474‑M1/134‑56 (MQ=255)
                                                       tGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcg                                                                                   >  1:634319‑M1/1‑59 (MQ=255)
                                                       tGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcg                                                                                   <  2:634319‑M1/110‑52 (MQ=255)
                                                               gcGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtc                                                   >  1:326751‑M1/1‑51 (MQ=255)
                                                                cGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcg                                                  <  1:601451‑M1/134‑85 (MQ=255)
                                                                 gggCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcg                                                  >  1:714928‑M1/1‑49 (MQ=255)
                                                                   gcgcgcCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaa                                                >  1:120089‑M1/1‑47 (MQ=255)
                                                                        ccgcTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggcc                                          >  1:262558‑M1/1‑42 (MQ=255)
                                                                            ttGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggccgttg                                      >  2:65038‑M1/1‑38 (MQ=255)
                                                                                              gCGTGCAATGCGAAAGAAGAcgattataaccgtggcgg                                                                                                                    <  2:1010622‑M1/38‑19 (MQ=255)
                                                                                              gCGTGCAATGCGAAAGAAGAcgattataaccgtggcgg                                                                                                                    >  1:1010622‑M1/1‑20 (MQ=255)
                                                                                                  gCAATGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggccgttgggcttgtcggcgaggaaacggg                >  2:983679‑M1/1‑16 (MQ=255)
                                                                                                     aTGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacag                                             >  1:770367‑M1/1‑13 (MQ=255)
                                                                                                     aTGCGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacag                                             <  2:770367‑M1/102‑90 (MQ=255)
                                                                                                        cGAAAGAAGAcgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggccgttgggcttgtcggcgaggaaacgggcggcg           <  2:262558‑M1/133‑124 (MQ=255)
                                                                                                                gacgattataaccgtggcggcaggccagtttcggccctgcgcgccccgcggcgccttacttcaggcccagcacgtcctgcatgtcgaacaggccgttgggcttgtcggcgaggaaacgggcggcgcgcacggca  <  2:326751‑M1/134‑133 (MQ=255)
                                                                                                                 |                                                                                                                                    
CAGGACAACGTTTGCGCGGCAACTGGCCTGCCGGCATCCGGGGGAATCGTTGGAATGAGAGCAGCGGGCGCGCCGCTTGCGCGGCGCCCGCCGGGCGTGCAATGCGAAAGAAGAGCCTCTTTTGAATTTCAAGTGGGTGGTATCGTGGCGGGTTTCAGGTAAAGAGCAACAAGAGCGATGCATAGCAAACTGTTTGAAGCAGCCCTTGGCATAAGCGATCCGTGGTTCGTCCGGGAAGTCGACTTC  >  AM260479/3399030‑3399275

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.