Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 556,858 A→T 100% L36Q (CTG→CAG)  folD ← bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb556,8580AT97.6% 119.6 / ‑3.9 42L36Q (CTG→CAG) folDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (21/20);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-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.

CAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCAC  >  W3110S.gb/556792‑556922
                                                                  |                                                                
caAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTcc                                                              <  1:3707845/71‑1 (MQ=255)
caAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTcc                                                              <  1:1708254/71‑1 (MQ=255)
caAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTcc                                                              <  1:294884/71‑1 (MQ=255)
caAGCCTTCCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGAACAACGGCCTGTcc                                                              >  1:1586201/1‑71 (MQ=255)
 aaGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCt                                                             <  1:1521720/71‑1 (MQ=255)
 aaGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCt                                                             <  1:344911/71‑1 (MQ=255)
  aGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTg                                                            >  1:993685/1‑71 (MQ=255)
   gCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTgg                                                           >  1:1703824/1‑71 (MQ=255)
        gCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc                                                               >  1:3100836/1‑62 (MQ=255)
         cGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCg                                                     <  1:4503651/71‑1 (MQ=255)
             tgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCa                                                   >  1:2262465/1‑69 (MQ=255)
             tgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGt                                                 >  1:4431163/1‑71 (MQ=255)
                     aCATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGc                                         <  1:4825801/71‑1 (MQ=255)
                      cATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGc                                         >  1:4346014/1‑70 (MQ=255)
                      cATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGc                                         >  1:4084919/1‑70 (MQ=255)
                         aaaTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATa                                     >  1:2093507/1‑71 (MQ=255)
                          aaTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATAc                                    >  1:2053646/1‑71 (MQ=255)
                                 gATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTg                             <  1:3923251/71‑1 (MQ=255)
                                   tGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGc                                         <  1:3725363/57‑1 (MQ=255)
                                     cAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACt                         <  1:3772065/71‑1 (MQ=255)
                                     cAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACt                         <  1:4252253/71‑1 (MQ=255)
                                        ggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACtttt                      >  1:3179353/1‑71 (MQ=255)
                                         gTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt                     <  1:2931349/71‑1 (MQ=255)
                                         gTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt                     <  1:2971026/71‑1 (MQ=255)
                                          ttACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTg                    <  1:5044687/71‑1 (MQ=255)
                                            actacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACtt                        >  1:3156346/1‑65 (MQ=255)
                                            actacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAg                  >  1:1769919/1‑71 (MQ=255)
                                            actacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTACGGCTGCAATACGCGCCTGAACTTTTTGAg                  >  1:874238/1‑71 (MQ=255)
                                                 ccaccaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACtt                        <  1:583327/60‑1 (MQ=255)
                                                 ccaccaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACt             <  1:400529/71‑1 (MQ=255)
                                                    ccaGCACAACGGCCTGTCCTGGTGCCCGCAGCCCGGCTGCAATACGCGCCTGAACTTTTTGAGc                 <  1:1907953/64‑1 (MQ=255)
                                                     caGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAAt                                      >  1:2992205/1‑42 (MQ=255)
                                                     caGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCag         <  1:4671621/71‑1 (MQ=255)
                                                      aGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCa          >  1:1336162/1‑69 (MQ=255)
                                                      aGCACAAAGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCaa               >  1:367195/1‑64 (MQ=255)
                                                       gCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCagag       <  1:3230937/71‑1 (MQ=255)
                                                        cacaACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCTACTTCAGAgc      >  1:3214373/1‑71 (MQ=255)
                                                          caACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGa                   <  1:4899103/56‑1 (MQ=255)
                                                            aCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCAc  >  1:3877655/1‑71 (MQ=255)
                                                              ggCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCaa               >  1:399764/1‑56 (MQ=255)
                                                              ggCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCaa               >  1:2239426/1‑56 (MQ=255)
                                                               gCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCaa               >  1:3822016/1‑55 (MQ=255)
                                                                  |                                                                
CAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCAC  >  W3110S.gb/556792‑556922

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