Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,229,470 A→G 14.3% Q684R (CAG→CGG)  metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,229,4700AG14.3% 65.3 / 4.6 28Q684R (CAG→CGG) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base A (10/14);  new base G (2/2);  total (12/16)
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.53e-01

TGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTG  >  W3110S.gb/4229407‑4229540
                                                               |                                                                      
tGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGc                                                                  >  1:1729533/1‑70 (MQ=255)
  gAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGc                                                                  <  1:3217941/68‑1 (MQ=255)
  gAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCt                                                                 <  1:2892257/69‑1 (MQ=255)
  gAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTa                                                                >  1:2632766/1‑70 (MQ=255)
   aaTACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTAcg                                                              <  1:2159486/71‑1 (MQ=255)
          gCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCg                                                        <  1:2082162/70‑1 (MQ=255)
               tCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTg                                                    >  1:1927712/1‑69 (MQ=255)
                cAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTg                                                    >  1:4119698/1‑68 (MQ=255)
                  aaGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTg                                               <  1:1719637/71‑1 (MQ=255)
                      cATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCGGCAGGCTACGCGCCCGATTGAAGTGATTg                                           <  1:357997/71‑1 (MQ=255)
                      cATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTg                                           <  1:311780/71‑1 (MQ=255)
                          aCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAgg                                       >  1:2290444/1‑71 (MQ=255)
                          aCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAgg                                       >  1:865840/1‑71 (MQ=255)
                           ccGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGc                                      <  1:3684204/71‑1 (MQ=255)
                                tttATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTg                                           <  1:4223437/61‑1 (MQ=255)
                                tttATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGtt                                 >  1:2766733/1‑71 (MQ=255)
                                 ttATAGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGtt                                 <  1:4200430/70‑1 (MQ=255)
                                  tATCGAGCAGGATACCGAAGAAGCCCGCCGGCAGGCTACGCGCCCGATTGAAGTGATTGaa                                         >  1:3294542/1‑61 (MQ=255)
                                  tATCGAGCAGGATACCGAAGAAGCCCGCCGGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTg                                <  1:2065994/70‑1 (MQ=255)
                                  tATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGa                               >  1:751862/1‑71 (MQ=255)
                                      gAGCAGGATACCGAAGAAGCCCGCCGGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGa                           >  1:1261536/1‑71 (MQ=255)
                                              tACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGaa                                         >  1:1498222/1‑49 (MQ=255)
                                              tACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATg                    >  1:128152/1‑70 (MQ=255)
                                              tACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGa                   <  1:4440383/71‑1 (MQ=255)
                                                ccGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTg                                           <  1:94760/45‑1 (MQ=255)
                                                     gaagCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGt            <  1:4059660/71‑1 (MQ=255)
                                                      aagCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCa                      <  1:1342999/60‑1 (MQ=255)
                                                       agCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCg          <  1:1354650/71‑1 (MQ=255)
                                                               agcagGCTCCGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTg  >  1:508350/1‑71 (MQ=255)
                                                               agcagGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTg  >  1:1717128/1‑71 (MQ=255)
                                                               |                                                                      
TGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTG  >  W3110S.gb/4229407‑4229540

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