Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,096,877 G→T 100% Y232* (TAC→TAA glpF ← glycerol facilitator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,096,8770GT96.8% 98.9 / ‑2.8 32Y232* (TAC→TAAglpFglycerol facilitator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (14/16);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACC  >  REL606/4096733‑4097023
                                                                                                                                                |                                                                                                                                                  
cGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAgg                                                                                                                                               <  1:16411/150‑1 (MQ=255)
cGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAgg                                                                                                                                               <  1:59585/150‑1 (MQ=255)
    gCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATg                                                                                                                                           >  1:678960/1‑150 (MQ=255)
       ttttGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGtct                                                                                                                                        <  1:194065/150‑1 (MQ=255)
           gTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTgc                                                                                                                                    >  1:457215/1‑150 (MQ=255)
                                  ttttCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGccc                                                                                                             <  2:651068/150‑1 (MQ=255)
                                   tttCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGcccc                                                                                                            <  2:680956/150‑1 (MQ=255)
                                       ttcCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGcc                                                                                                        <  2:39462/150‑1 (MQ=255)
                                                             ggCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTACGGAATGACTCTGCCGCCGGTAAAGGCGACATTGCCCCAACCCGCCAGCCAGGCAAACACTTTc                                                                                  >  1:406304/1‑150 (MQ=255)
                                                             ggCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTc                                                                                  >  2:601336/1‑150 (MQ=255)
                                                             ggCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATGGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCg                                                                                 >  1:424335/1‑150 (MQ=255)
                                                                         ccAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTTCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCa                                                                      <  1:336250/150‑1 (MQ=255)
                                                                                gTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTg                                                               >  2:650884/1‑150 (MQ=255)
                                                                                gTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTg                                                               >  2:147751/1‑150 (MQ=255)
                                                                                          agacagtGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGc                                                     <  2:252195/144‑1 (MQ=255)
                                                                                                tGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCAGGGcctgtct                                               >  1:251475/1‑144 (MQ=255)
                                                                                                    ccTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTctgtctcttatacaca                                           >  1:226456/1‑135 (MQ=255)
                                                                                                    ccTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTc                                           >  2:94679/1‑150 (MQ=255)
                                                                                                       acgattgagcCAACTATAGGGCCAACAAGCGGCACCAGGAAGTAAGTAATGTCTCTGCCGCCGGTAAATGTAAAAGTGCCCCAGCCCGCCAGCCCGGCAAAGACTTTCGGAGCGAAATCACGCGCTGGGTTCATCTCAAAACCTGTAAAt                                        <  1:425437/142‑1 (MQ=255)
                                                                                                                 cAACGATAGGGCCAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGa                             <  2:425543/150‑1 (MQ=255)
                                                                                                                  aaCGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCACAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGa                             <  2:429017/150‑1 (MQ=255)
                                                                                                                  aaCGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGa                             <  2:119133/150‑1 (MQ=255)
                                                                                                                  aaCGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGa                             >  2:588285/1‑150 (MQ=255)
                                                                                                                  aaCGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGa                             <  2:239824/150‑1 (MQ=255)
                                                                                                                  aaCGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGa                             <  2:230662/150‑1 (MQ=255)
                                                                                                                  aaCGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGa                             <  2:458517/150‑1 (MQ=255)
                                                                                                                  aaCGAAAGGGCATAAAAGCGGCACCAGGAATTAAGGAATGTTTCTGCCGCCGGTAAAGGCGACATTGCCCCAGACCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGa                             <  1:215981/150‑1 (MQ=255)
                                                                                                                    cGATAGGGCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATg                           >  2:28288/1‑150 (MQ=255)
                                                                                                                           gCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATGCGCCAAt                    >  1:498311/1‑150 (MQ=255)
                                                                                                                                            ggAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGGAAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGAc   <  1:147327/150‑1 (MQ=255)
                                                                                                                                            ggAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGAc   >  1:250303/1‑150 (MQ=255)
                                                                                                                                             gAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGAcc  >  1:377588/1‑150 (MQ=255)
                                                                                                                                             gAAATAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCACCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGACGTCACGCGCTGGGCAAATGGCAAAACCCGTCAAAGGCCCCACCATGGCGCCAATGAGCGCCAGCAGTaaccc  >  1:311006/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACC  >  REL606/4096733‑4097023

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: