Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,485,691 T→C 15.1% R703R (CGT→CGC ppc → phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,485,6910TC15.1% 78.3 / 7.9 33R703R (CGT→CGCppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (14/14);  new base C (3/2);  total (17/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.58e-01

ACTTCCGCTCCGCTACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCAAAACCGTCTGA  >  W3110S.gb/3485624‑3485761
                                                                   |                                                                      
aCTTCCGCTCCGCTACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAAcgccgc                                                                     >  1:4592662/1‑71 (MQ=255)
 cTTCCGCTCCGCTACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGcc                                                                    >  1:1577716/1‑71 (MQ=255)
 cTTCCGCTCCGCTACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGcc                                                                    >  1:4084631/1‑71 (MQ=255)
 cTTCCGCTCCGCTACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGcc                                                                    >  1:4043643/1‑71 (MQ=255)
    ccgctccgctACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCa                                                                  <  1:3444885/70‑1 (MQ=255)
    ccgctccgctACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCa                                                                  <  1:3438926/70‑1 (MQ=255)
    ccgctccgctACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGCCGccc                                                                   <  1:5116675/69‑1 (MQ=255)
            ctACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGCCGCCCAACCggcggc                                                         >  1:4477438/1‑71 (MQ=255)
            ctACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGCCGCCCAACCggcggc                                                         >  1:3047278/1‑71 (MQ=255)
               cGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTc                                                      <  1:3539932/71‑1 (MQ=255)
               cGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTc                                                      <  1:3533805/71‑1 (MQ=255)
                   ggAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGt                                                  <  1:865792/71‑1 (MQ=255)
                      aCAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCAc                                               <  1:2978920/71‑1 (MQ=255)
                      aCAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCAc                                               <  1:3924421/71‑1 (MQ=255)
                      aCAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCAc                                               <  1:1495854/71‑1 (MQ=255)
                       cAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTc                                                      <  1:1598409/63‑1 (MQ=255)
                            aCTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCCCTAcgcg                                         >  1:4748936/1‑71 (MQ=255)
                              tGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACgcgc                                        >  1:782542/1‑70 (MQ=255)
                              tGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACgcgc                                        >  1:2760396/1‑70 (MQ=255)
                               gggCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGgcggcg                                                        <  1:4630686/53‑1 (MQ=255)
                                   aaaCTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGt                                                  <  1:3114383/55‑1 (MQ=255)
                                   aaaCTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTcc                                  <  1:1607899/71‑1 (MQ=255)
                                         ccGTTGGGTTCACGTCCGGCGAAACGCCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTgg                              <  1:3089197/69‑1 (MQ=255)
                                               ggTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGc                      <  1:4749436/71‑1 (MQ=255)
                                               ggTTCACGTCCGGCGAAACGCCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCg                                 >  1:4333825/1‑60 (MQ=255)
                                                   cACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGc                      >  1:4526917/1‑67 (MQ=255)
                                                    aCGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTg                               >  1:3440241/1‑57 (MQ=255)
                                                    aCGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGa                 >  1:1291185/1‑71 (MQ=255)
                                                           gCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCaaaa          >  1:4697208/1‑71 (MQ=255)
                                                           gCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCaaaa          >  1:1919129/1‑71 (MQ=255)
                                                                aCGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCAAAACCGTc     >  1:4722652/1‑71 (MQ=255)
                                                                aCGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCAAAACCGTc     <  1:3302900/71‑1 (MQ=255)
                                                                aCGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCAAAACCGTc     >  1:1463760/1‑71 (MQ=255)
                                                                aCGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCAAAACCGTc     >  1:743342/1‑71 (MQ=255)
                                                                   tCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCAAAACCGTCTGa  <  1:4695663/71‑1 (MQ=255)
                                                                   |                                                                      
ACTTCCGCTCCGCTACGCCGGAACAAGAACTGGGCAAACTGCCGTTGGGTTCACGTCCGGCGAAACGTCGCCCAACCGGCGGCGTCGAGTCACTACGCGCCATTCCGTGGATCTTCGCCTGGACGCAAAACCGTCTGA  >  W3110S.gb/3485624‑3485761

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