Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,314,009 A→G 26.7% intergenic (‑380/+1) yciG ← / ← trpA hypothetical protein/tryptophan synthase subunit alpha

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,314,0090A→G26.7% 2.2 ‑3.7 0.027intergenic (‑380/+1)yciG/trpAhypothetical protein/tryptophan synthase subunit alpha
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (10/9):  total (18/9)
Fisher's exact test strand distribution p-value = 2.58e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TAAAATGCCGCCAGCGGAACTGGCGGCTGTGGGATTAACTGCGCGTCGCCGCTTTCATCGGTTGTACAA  >  REL606/1313976‑1314044
                                 |                                   
tAAAATGCCGCCAGCGGAACTGGCGGCTGTGGGAtt                                   <  1:4244489/36‑1 (MQ=255)
tAAAAGTCCGCCAGCGGATCGGGCGGCTGTGGGAtt                                   <  1:3257085/36‑1 (MQ=16)
    aTGCCGCCAGCGGAACTGGCGGCTGTGGGGTTAACt                               >  1:6227052/1‑36 (MQ=255)
    aTGCCGCCAGCGGAACTGGCGGCTGTGGGGTTAACt                               >  1:4747106/1‑36 (MQ=255)
      gccgccAGCGGAACTGGCGGCTGTGGGGTTTACTgc                             >  1:5024701/1‑36 (MQ=25)
        cgccAGCGGAACTGGCGGCTGTGGGATTAACTgcgc                           >  1:7472865/1‑36 (MQ=255)
        cgccAGCGGAACTGGCGGCTATGGGGTTAAATgcgg                           >  1:1416614/1‑35 (MQ=21)
          ccAGCGGAACTGGCGGCTGTGGGATTAACTGCGCGt                         >  1:2359436/1‑36 (MQ=255)
           cAGCGGAACTGGCGGCTGTGGGGTTAACTGCGCGTc                        >  1:3565226/1‑36 (MQ=255)
           cAGCGGAACTGGCGGCTGTGGGGTTAACTGCGCGTc                        >  1:7231196/1‑36 (MQ=255)
              cGGAACTGGCGGCTGTGGGATTAACTGCGCGTcgcc                     <  1:3206557/36‑1 (MQ=255)
              cGGAACTGGCGGCTGTGGGATTAACTGCGCGTcgcc                     <  1:6806812/36‑1 (MQ=255)
               ggAACTGGCGGCTGTGGGGTTTACTGCGCGTcgccg                    >  1:3204166/1‑36 (MQ=37)
               ggAACTGGCGGCTGTGGGGTTAACTGCGCGTcgccc                    >  1:6253740/1‑35 (MQ=255)
                gAACTGGCGGCTGTGGGATTAACTGCGCGTcggcgg                   >  1:7283510/1‑35 (MQ=255)
                   cTGGCGGCTGTGGGATTAACTGCGCGTCGCCGCttt                <  1:5804996/36‑1 (MQ=255)
                   cTGGCGGCTGTGGGATTAACTGCGCGTCGCAGCttt                >  1:3300056/1‑36 (MQ=255)
                    tGGCGGCTGTGGGATTAACTTCGCGTCGCCGCTTTc               >  1:288401/1‑36 (MQ=255)
                    tGGCGGCTGTGGGATTAACTGCGCGTCGCCGCTTTc               >  1:3636878/1‑36 (MQ=255)
                     ggcggcTGTGGGATTAACCGCGCGTCGCCGCTTTCa              >  1:7439859/1‑36 (MQ=255)
                        ggcTGTGGGATTAACTGCGCGTCGCCGCTTTCATCg           >  1:924423/1‑36 (MQ=255)
                          cTGTGGGATTAACTGCGCGTCGCCGCTTTCATCGGt         <  1:5392514/36‑1 (MQ=255)
                             tgGGATTAACTGCGCGTCGCCGCTTTCATCGGTTGt      <  1:6394122/36‑1 (MQ=255)
                             tgGGATTAACTGCGCGTCGCCGCTTTCATCGGTTGt      <  1:3575391/36‑1 (MQ=255)
                              gggATTAACTGCGCGTCGCCGCTTTCATCGGTTGTa     >  1:6323282/1‑36 (MQ=255)
                               ggATTAACTGCGCGTCGCCGCTTTCATCGGTTGTAc    <  1:3050786/36‑1 (MQ=255)
                               ggATTAACTGCGCGTCGCCGCTTTCATCGGTTGTAc    >  1:3134184/1‑36 (MQ=255)
                                 aTTAACTGCGCGTCGCCGCTTTCATCGGTTGTACaa  >  1:3437527/1‑36 (MQ=255)
                                 |                                   
TAAAATGCCGCCAGCGGAACTGGCGGCTGTGGGATTAACTGCGCGTCGCCGCTTTCATCGGTTGTACAA  >  REL606/1313976‑1314044

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 15 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG