Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,935,986 C→T 81.5% M111I (ATG→ATA edd ← 6‑phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,935,9860CT81.5% 57.6 / 7.8 27M111I (ATG→ATAedd6‑phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (18/4);  total (20/7)
Fisher's exact test for biased strand distribution p-value = 9.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTT  >  W3110S.gb/1935917‑1936053
                                                                     |                                                                   
cATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATcacaca                                                                    >  1:820331/1‑71 (MQ=255)
cATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATcacaca                                                                    >  1:1454875/1‑71 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:1518262/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:364360/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:71862/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:733323/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:1799427/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:883634/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:1333304/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:1309745/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:886048/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:978370/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGAcc                                        >  1:99084/1‑70 (MQ=255)
                             aCAATTCCATTCCATCCTGCCCCTGGGAGACACCATCACATATCGCCGGAACACCGCCCGCAAc                                              >  1:1909779/1‑64 (MQ=255)
                                 ttccattccaTCCTGCCCCTGGGTGACACCATCACATATCGCCGGAac                                                          <  1:1222805/48‑1 (MQ=255)
                                 ttccattccaTCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCg                                   <  1:1677449/71‑1 (MQ=255)
                                 ttccattccaTCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCg                                   <  1:181570/71‑1 (MQ=255)
                                 ttccattccaTCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCg                                   <  1:714328/71‑1 (MQ=255)
                                        ccatCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGc                            >  1:266261/1‑71 (MQ=255)
                                        ccatCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGc                            >  1:1290289/1‑71 (MQ=255)
                                        acatCCTGCCCCTGGGTGACACCATCACATATCGCCGGAACAccgcc                                                    >  1:1111681/2‑47 (MQ=39)
                                                     gggTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCt               >  1:1807836/1‑71 (MQ=255)
                                                     gggTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCt               >  1:1334405/1‑71 (MQ=255)
                                                     gggTGACACCATCACATATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCt               >  1:1208061/1‑71 (MQ=255)
                                                         gACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTAc           >  1:785309/1‑71 (MQ=255)
                                                         gACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTAc           >  1:929807/1‑71 (MQ=255)
                                                                tCACATATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGAt    <  1:682624/71‑1 (MQ=255)
                                                                tCACATATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGAt    <  1:300438/71‑1 (MQ=255)
                                                                tCACATATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGAt    <  1:1930961/71‑1 (MQ=255)
                                                                  acacaTCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGAttt  >  1:959652/1‑71 (MQ=255)
                                                                  acacaTCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGAtt   >  1:55193/1‑70 (MQ=255)
                                                                  acacaTCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGAtt   >  1:1027778/1‑70 (MQ=255)
                                                                  acacaTCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGAtt   >  1:873492/1‑70 (MQ=255)
                                                                     |                                                                   
CATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTT  >  W3110S.gb/1935917‑1936053

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: