Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,423,616 G→C 31.9% G266G (GGC→GGG yfdV ← predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,423,6160G→C31.9% 2.7 ‑1.9 0.020798yfdVpredicted transporter
Reads supporting (aligned to +/- strand):  new base (1/5):  ref base (7/5):  total (9/11)
Fisher's exact test strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATAATGATTCCGGAGAATGCCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGATGTTCGCTGTT  >  REL606/2423584‑2423651
                                |                                   
ataatGATTCAGGAGAATGCCGGCGGTAATGCGCCt                                  <  1:2666423/36‑1 (MQ=255)
   attaTTCCGGAGAATTCCGGCGGTAATGCGCCTGcc                               <  1:5070240/33‑1 (MQ=37)
   atGATTCCGGAGACGGCCGGCGGGACNNNCCCTGcc                               >  1:1995229/1‑36 (MQ=12)
   atGATTCCGGAGAATGCCGGCGGTAATGCGCCTGcc                               <  1:3993764/36‑1 (MQ=255)
      aTTCCGGAGAATGCCGGCGGTAATTCGCCCGcccgc                            >  1:1446728/1‑33 (MQ=37)
      aTTCCGGAGAATGCCGGCGGTAATTCACCTAccccg                            >  1:2138478/1‑33 (MQ=19)
         ccGGAGAATGCCGGCGGTAAATCGCCTGCCCGTAcc                         >  1:4919506/1‑36 (MQ=21)
            gagaATGCCGGCGGTGATGAGCCTGCCCGCAcccac                      >  1:3253598/1‑33 (MQ=12)
               aaggCCGGCGGTAATGCTCCTGCCAGTACCatcatc                   <  1:1493840/33‑1 (MQ=37)
                aTGCCGGCGGTAATCCCCCTGCCAGTACCATCATCt                  <  1:2520686/36‑1 (MQ=37)
                 tGCCGGCGGTAAGGCGCCTGCCAGTACCATCATCTg                 <  1:3717034/36‑1 (MQ=255)
                  gCCGGCGGTAATGCCCCTGCCAGTACCATCATCTGc                <  1:3790145/36‑1 (MQ=255)
                    cggcggAAATGCCCCTGCCAGTACCATCATCTGCAg              <  1:4433534/36‑1 (MQ=37)
                     ggcggTAATGCGCCTGCCAGTACCATCATCTGCAGa             >  1:5648002/1‑36 (MQ=255)
                        ggTAATGCGCCTGCCAGTACCATCATCTGCAGATGt          >  1:1955017/1‑36 (MQ=255)
                        ggTAATGCGCCTGCCAGTACCATCATCTGCAGATGt          >  1:1940453/1‑36 (MQ=255)
                         gTAAGGCCCCTGCCAGTACCATCATCTGCAGATGtt         <  1:2135762/36‑1 (MQ=37)
                         atAATGCCCCTGCCAGTACCATCATCTGCAGATGtt         <  1:1713914/35‑1 (MQ=255)
                            aTGCGCCTGCCAGTACCATCATCTGCAGATGTTCGc      >  1:5591213/1‑36 (MQ=255)
                             tGCGCCTGCCAGTACCATCATCTGCAGATGTTCGCt     <  1:915079/36‑1 (MQ=255)
                               cgcCTGCCAGTACCATCATCTGCAGATGTTCGCTg  >  1:2356985/1‑35 (MQ=255)
                                gcCTGCCAGTACCATCATCTGCAGATGTTCGCTGtt  >  1:638780/1‑36 (MQ=255)
                                gcCTGCCAGTACCATCATCTGCAGATGTTCGCTGtt  >  1:1440438/1‑36 (MQ=255)
                                |                                   
ATAATGATTCCGGAGAATGCCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGATGTTCGCTGTT  >  REL606/2423584‑2423651

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 10 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG