Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,855,581 T→A 15.5% intergenic (+2/+52) yieF → / ← yieG chromate reductase, Class I, flavoprotein/predicted inner membrane protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,855,5810T→A15.5% 2.0 ‑2.7 0.030intergenic (+2/+52)yieF/yieGchromate reductase, Class I, flavoprotein/predicted inner membrane protein
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (12/13):  total (12/18)
Fisher's exact test strand distribution p-value = 6.57e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGTAATCATGCATGGCGG  >  REL606/3855546‑3855611
                                   |                              
aTTTGGTGAGTTTATTCAGCGAGTTAAGATCtaaat                                <  1:2096687/36‑1 (MQ=255)
   tGGTGAGTTTATTCAGCGAGTTAAGATCtaaataaa                             >  1:2065104/1‑36 (MQ=255)
   tGGTGAGTTTATTCAGCGAGTTAAGATCtaaataaa                             >  1:7013548/1‑36 (MQ=255)
     gTGAGTTTATTCAGCGAGTTAAGATCTAAATaaaaa                           <  1:6156515/36‑1 (MQ=255)
      tGAGTTTATTCAGCGAGTTAAGATCTAAATaaaaaa                          >  1:4369634/1‑36 (MQ=255)
      tGAGTTTATTAAGCGAGTTAAGATCTAAATaaaaaa                          <  1:6682160/36‑1 (MQ=255)
           ttATTCAGCGAGTTAAGATCTAAATAAAAACCCCGc                     <  1:3573563/36‑1 (MQ=255)
           ttATTCAGCGAGTTAAGATCTAAATAAAAAACCCGc                     <  1:6147130/36‑1 (MQ=255)
           ttATTCAGCGAGTTAAGATCTAAATAAAAAACCCGc                     <  1:6505346/36‑1 (MQ=255)
           ttATTCAGCGAGTTAAGATCTAAATAAAAAACCCGc                     >  1:56333/1‑36 (MQ=255)
            tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGcc                    <  1:2215953/36‑1 (MQ=255)
            tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGcc                    <  1:7256557/36‑1 (MQ=255)
            tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGcc                    >  1:8203627/1‑36 (MQ=255)
            tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGcc                    <  1:30912/36‑1 (MQ=255)
            tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGcc                    <  1:3877577/36‑1 (MQ=255)
            tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGcc                    <  1:4052800/36‑1 (MQ=255)
            tATTCAGCGAGTTAAGATCTAAATAAAAAACCCGcc                    >  1:6603786/1‑36 (MQ=255)
             aTTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCa                   >  1:3742429/1‑36 (MQ=255)
             aTTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCa                   <  1:2237967/36‑1 (MQ=255)
             aTTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCa                   >  1:981231/1‑36 (MQ=255)
              ttCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAg                  <  1:2985095/36‑1 (MQ=255)
               aaaGCGAGTTAAGATCTAAAAAAAAAACCCGCCAgc                 <  1:611198/34‑2 (MQ=37)
                cAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGTa                >  1:4370630/1‑36 (MQ=255)
                cAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGTa                >  1:1461950/1‑36 (MQ=255)
                 aGCGATTAAAAATATAAAAAAAAACCCCGCCAGTaa               <  1:2659584/36‑1 (MQ=9)
                    gAGTTAAGATATAAATAAAAAACCCGCCAGTAATca            <  1:55078/36‑1 (MQ=255)
                      gTTAAGATCTAAATAAAAAACCCGCCAGTAATcatg          >  1:4078567/1‑36 (MQ=255)
                      attAAGACATAAAAAAAAAACCCGCCAGTAATcatg          <  1:8295123/35‑1 (MQ=21)
                       taaaGATAAAAAAAAAAACCCCGCCAGTAATcatgc         <  1:7092255/34‑1 (MQ=16)
                         aaaaTCAAAAAAAAAATCCCGCCAGTAATcatgcat       <  1:1754842/33‑1 (MQ=19)
                          aGATCTAAATAAAAAACCCGCCAGTAATcatgcatg      >  1:5965710/1‑36 (MQ=255)
                           gATCTAAATAAAAAACCCGCCAGTAATCATGCATgg     >  1:6663456/1‑36 (MQ=255)
                            aTCTAAATAAAAAACCCGCCAGTAATCATGCATGGc    >  1:1368125/1‑36 (MQ=255)
                              ttaaaaaaaAAGCCCGCCAGTAATCATGCATGGCgg  <  1:5252363/35‑1 (MQ=25)
                              aaaaaaaaaAAACCCGCCAGTAATCATGCATGGCgg  <  1:2750420/34‑1 (MQ=37)
                                   |                              
ATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGTAATCATGCATGGCGG  >  REL606/3855546‑3855611

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