Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,341,364 G→A 19.9% intergenic (+9/+117) yciH → / ← osmB translation intiation factor Sui1/lipoprotein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,341,3640G→A19.9% 4.4 ‑3.2 0.025intergenic (+9/+117)yciH/osmBtranslation intiation factor Sui1/lipoprotein
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (8/11):  total (13/12)
Fisher's exact test strand distribution p-value = 4.11e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GAAGGTAAAACTCGCAGGCGGTTAACATAAAAAGCCACGGATATATCCGTGGCTTTCGAATATTTTACT  >  REL606/1341331‑1341399
                                 |                                   
gAAGGTAAAACTCGCAGGCGGTTAACATAAAAAGcc                                   <  1:3027582/36‑1 (MQ=255)
gAAGGTAAAACTCGCAGGCGGTTAACATAAAAAAcc                                   >  1:5795931/1‑36 (MQ=255)
 aaGGTAAAACTCGCAGGCGGTTTAAATAAAAAAccc                                  >  1:2784269/1‑35 (MQ=25)
  aGGTAAAACTCGCAGGCGGTTAACATAAAAAGCCAc                                 <  1:8091222/36‑1 (MQ=255)
  aGGTAAAACTCGCAGGCGGTTAACATAAAAAGCCAc                                 <  1:2979790/36‑1 (MQ=255)
    gTAAAACTCGCAGGCGGTTAACATAAAAAGCCACgg                               <  1:6108591/36‑1 (MQ=255)
     tAAAACTCGCAGGCGGTTAAAATAAAAAACCACGGa                              >  1:1798385/1‑36 (MQ=37)
         aCTCGCAGGCGGTTAACATAAAAAGCCACGGatata                          <  1:4188613/36‑1 (MQ=255)
           tCGCAGGCGGTTTACATAAAAAACCACCGATATATc                        >  1:4969405/1‑36 (MQ=25)
           tCGCAGGCGGTTAACAAAAAAAACCCCGGGTATATc                        >  1:2410473/1‑36 (MQ=19)
            cGCAGGCGGTTAAAATAAAAAGCCACGGATATATcc                       <  1:3418122/36‑1 (MQ=255)
              cAGGCGGTTAACATAAAAAGCCACGGATATATCCGt                     >  1:3100162/1‑36 (MQ=255)
              cAGGCGGTTAACATAAAAAGCCACGGATATATCCGt                     >  1:3105889/1‑36 (MQ=255)
                ggCGGTTAACATAAAAAGCCACGGATATATCCGTgg                   <  1:3167647/36‑1 (MQ=255)
                ggCGGTTAACATAAAAAGCCACGGATATATCCGTgg                   >  1:6055023/1‑36 (MQ=255)
                    gTTAACATAAAAAGCCACGGATATATCCGTGGCttt               >  1:7935449/1‑36 (MQ=255)
                    gTTAACATAAAAAGCCACGGATATATCCGTGGCttt               <  1:1581457/36‑1 (MQ=255)
                     ttAACATAAAAAGCCACGTATATATCCGTGGCTTTc              <  1:1767626/36‑1 (MQ=255)
                     ttAACATAAAAAGCCACGGATATATCCGTGGCTTTc              >  1:3841723/1‑36 (MQ=255)
                     ttAACATAAAAAGCCACGGATATAGCCGTGGAtttt              >  1:3737652/1‑35 (MQ=18)
                         cATAAAAAGCCACGGATATATCCGTGGCTTTCGAat          <  1:4672045/36‑1 (MQ=255)
                         cATAAAAAGCCACGGATATATCCGTGGCTTTCGAat          >  1:6789345/1‑36 (MQ=255)
                         aaTAAAAATTCACGGATATATCCGTGGCTTTCGAat          <  1:5262213/35‑1 (MQ=18)
                         aaTAAAAAGCCACGGATATATCCGTGGCTTTCGAat          <  1:427274/35‑1 (MQ=255)
                            aaaaaGCCACGGATATATCCGTGGCTTTCGAATAtt       >  1:4667629/1‑36 (MQ=255)
                                 gCCACGGATATATCCGTGGCTTTCGAATATTTTACt  <  1:2178887/36‑1 (MQ=255)
                                 |                                   
GAAGGTAAAACTCGCAGGCGGTTAACATAAAAAGCCACGGATATATCCGTGGCTTTCGAATATTTTACT  >  REL606/1341331‑1341399

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG