Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 196,338 A→T 12.8% intergenic (+67/‑25) frr → / → ispC ribosome releasing factor/1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606196,3380A→T12.8% 2.2 ‑3.3 ‑0.030intergenic (+67/‑25)frr/ispCribosome releasing factor/1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (10/16):  total (14/16)
Fisher's exact test strand distribution p-value = 3.65e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.60e-01

TCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGATGTTTCATGAAGCAAC  >  REL606/196304‑196372
                                  |                                  
tCCTTGCGGATCGGCTGGGGGGGTTTTTTTTTTTAt                                   >  1:128563/1‑36 (MQ=19)
tCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTAt                                   <  1:5647255/36‑1 (MQ=255)
 ccTTGCGGATCGGCTGGGGGGGTTTTTCTTTTTTtt                                  >  1:3003526/1‑36 (MQ=16)
  cTTGCGGATCGGCTGGCGGCGTTTTTCTTTTTTttt                                 >  1:2448777/1‑35 (MQ=25)
   atGCGGATCTGCTGTCGGCGTTTTGCTTTTTATTCt                                <  1:898218/35‑1 (MQ=25)
    tGCGGATCGGCTGGCGGCGTTTTTTTTTTTTttttt                               >  1:4684570/1‑33 (MQ=16)
    tGCGGATCGGCTGGCGGCGTTTTGTTTTTTTttttt                               >  1:6115840/1‑33 (MQ=21)
      cGGATCGGGTGGCGGCGTTTTGCTTTTTATTCTGtc                             >  1:4128091/1‑36 (MQ=255)
        gATCGGCTGGCGGCGTTTTGCTTTTTATTCTGtctc                           <  1:2922398/36‑1 (MQ=255)
         aTCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCa                          <  1:5518501/36‑1 (MQ=255)
            ggCTGGCGGCGTTTTGCTTTTTATTCTTTCTCAAct                       >  1:1692392/1‑36 (MQ=255)
            ggCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAAct                       >  1:685626/1‑36 (MQ=255)
              cTGGCGGCGTTTTGCTTTTTATTCTGTCTCAActct                     <  1:2010475/36‑1 (MQ=255)
                ggcggcGTTTTGCTTTTTATTCTGTCTCAACTCTgg                   <  1:1410385/36‑1 (MQ=255)
                  cggcgTTTTGCTTTTTATTCTGTCTCAACTCTGGAt                 <  1:5570444/36‑1 (MQ=255)
                  cggcgTTTTGCTTTTTATTCTGTCTCAACTCTGGAt                 >  1:6045157/1‑36 (MQ=255)
                  cggcgTTTTGCTTTTTATTCTGTCTCAACTCTGGAt                 <  1:220676/36‑1 (MQ=255)
                       ttttGCTTTTTATTCTTTCTCAACTCTGGATGTTTc            >  1:1947500/1‑36 (MQ=255)
                        tttGCTTTTTATTCTGTCTCAACTCTGGATGTTTCa           >  1:5200119/1‑36 (MQ=255)
                        tttGCTTTTTATTCTGTCTCAACTCTGGATGTTTCa           <  1:6265771/36‑1 (MQ=255)
                         ttGCTTTTTATTCTGTCTCAACTCTGGGTGTTTCAt          <  1:4337612/36‑1 (MQ=255)
                         ttGCTTTTTATTCTGTCTCAACTCTGGATGTTTCAt          >  1:5689887/1‑36 (MQ=255)
                            cTTTTTATTCTGTCTCAACTCTGGATGTTTCATGaa       <  1:4467510/36‑1 (MQ=255)
                            cTTTTTATTCTGTCTCAACTCTGGATGTTTCATGaa       <  1:2741153/36‑1 (MQ=255)
                            cTTTTTATTCTGTCTCAACTCTGGATGTTTCATGaa       <  1:7218079/36‑1 (MQ=255)
                             tttttATTCTGTCTCAACTCTGGATGTTTCATGAAg      >  1:1318535/1‑36 (MQ=255)
                             tttttATTCTGTCTCAACTCTGGATGTTTCATGAAg      >  1:874166/1‑36 (MQ=255)
                              ttttATTCTGTCTCAACTCTGGATGTTTCATGAAGc     <  1:5586507/36‑1 (MQ=255)
                              ttttATTCTGTCTCAACTCTGGATGTTTCATGAAGc     <  1:2818330/36‑1 (MQ=255)
                                 tATTCTGTCTCAACTCTGGATGTTTCATGAAGCAAc  <  1:5672715/36‑1 (MQ=255)
                                  |                                  
TCCTTGCGGATCGGCTGGCGGCGTTTTGCTTTTTATTCTGTCTCAACTCTGGATGTTTCATGAAGCAAC  >  REL606/196304‑196372

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG