Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 3,147,765 Δ27 bp [ygjF] [ygjF]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3147765 3147791 27 33 [0] [0] 33 [ygjF] [ygjF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 31477640 (0.000)30 (0.800) 25/142 0.5 100% intergenic (+59/+20) rpoD/ygjF RNA polymerase sigma factor/G/U mismatch‑specific DNA glycosylase
?REL606 3147792 = 0 (0.000)coding (499/507 nt) ygjF G/U mismatch‑specific DNA glycosylase

TTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGT  >  REL606/3147691‑3147838
                                                                         |                                                                          
ttCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCca                                                                          <  1:129757‑M1/76‑3 (MQ=255)
      gACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcact                                                                    <  2:880016‑M1/76‑9 (MQ=255)
        cGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactac                                                                  >  1:1131661‑M1/1‑66 (MQ=255)
        cGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactac                                                                  >  1:19584‑M1/1‑66 (MQ=255)
             aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcg                                                             <  1:126842‑M1/76‑16 (MQ=255)
              aTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCGCTGCACAAACGCCcacgcactacccgcgc                                                            >  1:979284‑M1/1‑60 (MQ=255)
                    aGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtc                                                      <  1:597279‑M1/76‑23 (MQ=255)
                       ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccag                                                   >  1:325887‑M1/1‑51 (MQ=255)
                        cGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagc                                                  >  1:154970‑M1/1‑50 (MQ=255)
                           aTCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcg                                               <  1:221518‑M1/76‑30 (MQ=255)
                           aTCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcg                                               <  1:335388‑M1/76‑30 (MQ=255)
                           aTCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcg                                               >  1:521037‑M1/1‑47 (MQ=255)
                            tCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAAGACCgacgcactagcagcggctggtccagctagt                                              >  2:77451‑M1/1‑46 (MQ=255)
                             cAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcg                                             >  1:1124021‑M1/1‑45 (MQ=255)
                             cAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcg                                             <  1:358322‑M1/76‑32 (MQ=255)
                              aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcga                                            >  2:1082587‑M1/1‑44 (MQ=255)
                                  gTTACGCCGCCCCCGGCACTAGCGCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacg                                        <  1:549‑M1/76‑37 (MQ=255)
                                  gTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacg                                        <  1:472022‑M1/76‑37 (MQ=255)
                                    tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacgct                                      <  1:898864‑M1/76‑39 (MQ=255)
                                    tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacgct                                      >  2:493570‑M1/1‑38 (MQ=255)
                                    tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagcccccgatacgct                                      >  1:56854‑M1/1‑38 (MQ=255)
                                     aCGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacgctt                                     <  1:1159835‑M1/76‑40 (MQ=255)
                                      cgccgcACCCGGCACTAGGCCCTCTGCACAAACTCCcccgcactatcagcgcctgctccagcgcgcggtatgcttc                                    >  2:578069‑M1/1‑36 (MQ=255)
                                       gccgcACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacgcttca                                   <  1:452974‑M1/76‑42 (MQ=255)
                                        ccgcACCCGGCACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacgcttcaa                                  >  2:1039886‑M1/1‑34 (MQ=255)
                                                  cACTAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacgcttcaaccagtttctc                        <  2:1188778‑M1/76‑53 (MQ=255)
                                                     tAGGCCCTCTGCACAAACGCCcacgcactaccagcgcctggtccagctcgcgatacgcttcaaccagtttctccag                     <  2:459786‑M1/76‑56 (MQ=255)
                                                           ctctGCACAAACGCCcacgcactaccagcgcctggtccagctcacgatacgcttcaaccagtttctccagtgaaac               >  1:101773‑M1/1‑15 (MQ=255)
                                                                  caAACGCCcacgcactaccagcgcctggtccagctcgcgatacgcttcaaccagtttctccagtgaaacgcgactt        >  2:376459‑M1/1‑8 (MQ=255)
                                                                    aaCGCCcacgcactaccagcgcctggtccagctcgcgatacgcttcaaccagtttctccagtgaaacgcgacttaa      >  2:730630‑M1/1‑6 (MQ=255)
                                                                      cGCCcccgccctaccagcgcctggtccagctcgcgatacgcttcaaccagtttctccagtgaaacgcgacttaaac    <  2:101773‑M1/76‑73 (MQ=255)
                                                                       gCCcacgcactaccagcgcctggtccagctcgcgatacgcttcaaccagtttctccagtgaaacgcgacttaaacc   >  1:1244509‑M1/1‑3 (MQ=255)
                                                                        cccacgcactaccagcgcctggtccagctcgcgatacgcttcaaccagtttctccagtgaaacgcgacttaaaccg  >  1:652491‑M1/1‑2 (MQ=255)
                                                                         |                                                                          
TTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGT  >  REL606/3147691‑3147838

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