Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,580,806 Δ2,622 bp IS1‑mediated ryhByrhA ryhB, yhhY, yhhZ, yrhA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3580806 3584153–3583429 2624–3348 25 [0] [3] 5 ryhB–[insB1] ryhB,yhhY,yhhZ,yrhA,insA,[insB1]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 35808050 (0.000)23 (0.760) 16/130 0.5 100% coding (1/1038 nt) yhhX putative oxidoreductase
?NC_000913 3583428 = NA (NA)noncoding (1/768 nt) IS1 repeat region

AGTACATACGGCAGATGGTAACGGGTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATGATTTCGTCCTTTTTAAGGTGGTTATTTACACCTTAGCGCAAAGCGGACGTGGTTCCTACTGGAGT  >  NC_000913/3580738‑3580871
                                                                   |                                                                  
aGTACATACGGCAGATGGTAACGGGTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCAt                                                                    <  1:65947/68‑1 (MQ=255)
 gTACATACGGCAGATGGTAACGGGTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATg                                                                   <  2:957059‑M1/68‑2 (MQ=255)
      tACGGCAGATGGTAACGGGTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATggtgat                                                              <  1:929297‑M1/68‑7 (MQ=255)
        cGGCAGATGGTAACGGGTGGTGCTTTTGCCGCAGCCAATAAAGGCGCAGATGCTGACCATggtgatgc                                                            >  2:314587‑M1/1‑60 (MQ=255)
        cGGCAGATGGTAACGGGTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATggtgatgc                                                            >  2:864341‑M1/1‑60 (MQ=255)
        cGGCAGATGGTAACGGGTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATggtgatgc                                                            >  1:77219‑M1/1‑60 (MQ=255)
          gCAGATGGTAACGGGTGGTGCTTTTGCCGCAGCCAATGAAGGCGCAGTTGATGCCCATggtggagctg                                                          >  2:7892‑M1/1‑58 (MQ=255)
          gCAGATGGTAAAGGGTGGTGCTTTTGCCGAAGCAAATAAAGGCGCAGATGATGACCATggtgaagctg                                                          >  1:641463‑M1/1‑58 (MQ=255)
            agctggtAACGGGTGGTGCTTTTGCCGAAGCCAATAAATGCGCAGTTTATGACCATggtgatgctgcc                                                        <  1:760857‑M1/65‑13 (MQ=255)
               tGGTAACGGGTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATggtgatgctgccaac                                                     >  2:790628‑M1/1‑53 (MQ=255)
                    aCGGGTGGTGCTTTTGACGAAGCCAATAAAGGCGCAGTTGATGACCATggtgatgctgccaagttact                                                >  1:644599‑M1/1‑48 (MQ=255)
                      gggTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATggtgatgctgccaacttactga                                              <  1:1040945‑M1/68‑23 (MQ=255)
                       ggcggtgcttTTGACGAAGCCAATAAAGGCGCAGTTGATGACCATggtgatgctgccaacatactgat                                             >  2:367590‑M1/4‑45 (MQ=255)
                            tgCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATggtgatgctgccaacttactgatttagt                                        <  2:55034‑M1/68‑29 (MQ=255)
                              cTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATggtgatgctgccaacttactgatttagtgt                                      <  2:1042515‑M1/68‑31 (MQ=255)
                                        aGCCAATAAAGGCGCAGTTGATGACCATggtgatgctgccaacttactgatttagtgtatgatggtgt                            <  1:35588‑M1/68‑41 (MQ=255)
                                         gCCAATAACGGCGCAGTTGATGACCATggcgatgctgccaacttactgatttagtgtatgatggtgtt                           <  2:434004‑M1/68‑42 (MQ=255)
                                         gCCAATAAAGGCGCAGGTGATGACCATggtgatgctgccaacttactgatttagtgtatgatggtgtt                           <  2:511705‑M1/68‑42 (MQ=255)
                                                        gTTGATGACCATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccag            >  1:659730‑M1/1‑12 (MQ=255)
                                                          tgatgaCCATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtg          <  1:282391‑M1/68‑59 (MQ=255)
                                                             tgaCCATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtgtct       <  1:944237‑M1/68‑62 (MQ=255)
                                                               aCCATggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttc     >  2:738833‑M1/1‑5 (MQ=255)
                                                                  ttggagatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgt  <  1:523623‑M1/67‑67 (MQ=255)
                                                                  aTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgt  <  1:864344‑M1/68‑67 (MQ=255)
                                                                  aTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgt  >  1:410330‑M1/1‑2 (MQ=255)
                                                                   |                                                                  
AGTACATACGGCAGATGGTAACGGGTGGTGCTTTTGCCGAAGCCAATAAAGGCGCAGTTGATGACCATGATTTCGTCCTTTTTAAGGTGGTTATTTACACCTTAGCGCAAAGCGGACGTGGTTCCTACTGGAGT  >  NC_000913/3580738‑3580871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ 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.