Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,622,543 C→T 29.2% A362A (GCG→GCA)  metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,622,5430CT29.2% 33.2 / 15.2 24A362A (GCG→GCA) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (3/14);  new base T (4/3);  total (7/17)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTT  >  W3110S.gb/3622474‑3622601
                                                                     |                                                          
ggTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAgcg                                                           <  1:990767/71‑1 (MQ=255)
ggTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAgcg                                                           <  1:1215304/71‑1 (MQ=255)
ggTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAgcg                                                           <  1:219763/71‑1 (MQ=255)
ggTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAgcg                                                           <  1:1958938/71‑1 (MQ=255)
ggTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAgcg                                                           <  1:1586376/71‑1 (MQ=255)
       tGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATcgcg                                                    >  1:1055621/1‑71 (MQ=255)
       tGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATcgcg                                                    >  1:711168/1‑71 (MQ=255)
       tGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATcgcg                                                    >  1:1687055/1‑71 (MQ=255)
       tGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATcgc                                                     >  1:1079209/1‑70 (MQ=255)
       tGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                                                    <  1:1055596/71‑1 (MQ=255)
       tGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcgcg                                                    <  1:755445/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                <  1:1958622/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                <  1:768288/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                <  1:439215/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                <  1:1224280/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                <  1:1226544/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                <  1:1945376/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                <  1:1943012/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                <  1:1822641/71‑1 (MQ=255)
                     cttCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGtt                                      <  1:1945785/69‑1 (MQ=255)
                     cttCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGtt                                      <  1:194205/69‑1 (MQ=255)
                     cttCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGtt                                      <  1:1804910/69‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                         <  1:147376/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                         <  1:299743/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                         <  1:665509/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                         <  1:1604601/71‑1 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt  >  1:125634/1‑71 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt  >  1:1682254/1‑71 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt  >  1:903609/1‑71 (MQ=255)
                                                                     |                                                          
GGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTT  >  W3110S.gb/3622474‑3622601

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