Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,440,923 T→C 19.2% S114S (AGT→AGC thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,440,9230TC19.2% 48.9 / 7.1 26S114S (AGT→AGCthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (13/8);  new base C (3/2);  total (16/10)
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

GCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGCCTCGACGAACTGCGTTTTAGCGGCGTAC  >  W3110S.gb/3440857‑3440987
                                                                  |                                                                
gCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCg                                                              <  1:1737983/71‑1 (MQ=255)
    gCAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGAt                                                            <  1:211067/69‑1 (MQ=255)
    gCAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTa                                                          >  1:1056505/1‑71 (MQ=255)
     cAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGAt                                                            <  1:534710/68‑1 (MQ=255)
     cAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGAt                                                            <  1:1587239/68‑1 (MQ=255)
                     gTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTgg                                         >  1:2404583/1‑71 (MQ=255)
                     gTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTgg                                         >  1:1544631/1‑71 (MQ=255)
                        gATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTgg                                         >  1:1641230/1‑68 (MQ=255)
                          tCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAgat                                    <  1:65014/71‑1 (MQ=255)
                           cGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGCTCCGATTACACTAAAGCGCGGCTGGCAgat                                    >  1:1778768/1‑70 (MQ=255)
                           cGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGCTCCGATTACACTAAAGCGCGGCTGGCAgat                                    >  1:726371/1‑70 (MQ=255)
                           cGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGCTCCGATTACACTAAAGCGCGGCTGGCAgat                                    >  1:2870123/1‑70 (MQ=255)
                            gATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAgatga                                  <  1:2390447/71‑1 (MQ=255)
                            gATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAgatga                                  <  1:3436215/71‑1 (MQ=255)
                            gATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGCTCCGATTACACTAAAGCGCGGCTGGCAgatga                                  <  1:1016240/71‑1 (MQ=255)
                                cgcgcgGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGc                              >  1:1572641/1‑71 (MQ=255)
                                cgcgcgGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGc                              >  1:337449/1‑71 (MQ=255)
                                              gaagaaCTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGCCTCGACGaa                     >  1:3027534/1‑66 (MQ=255)
                                              gaagaaCTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGCCTCGACGAACTGCg                >  1:3385912/1‑71 (MQ=255)
                                                agaaCTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTTGCAGATGATGGCCTCGACGAACTGCGtt              >  1:295235/1‑71 (MQ=255)
                                                 gaaCTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGCCTCGACGAACTGCGttt             >  1:35898/1‑71 (MQ=255)
                                                  aaCTTACCGTGCGCAGCTCCGATTACACTAAAGCGCGGCTGGCAgatga                                  <  1:1606401/49‑1 (MQ=255)
                                                   aCTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGCCTCGACGAACTGCGtttt            <  1:3846827/70‑1 (MQ=255)
                                                    cTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGCCTCGACGAACTGCGTTTTAg          >  1:801672/1‑71 (MQ=255)
                                                          gTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGatgatg                                >  1:2161135/1‑43 (MQ=255)
                                                             cgcAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGCCTCGACGAACTGCGTTTTAGCGGCGTAc  >  1:1786643/1‑70 (MQ=255)
                                                                  |                                                                
GCTGGCAAAACTACGCCAGCCGTGGATCGATGCGCGCGGCGATACCGAAGAACTTACCGTGCGCAGTTCCGATTACACTAAAGCGCGGCTGGCAGATGATGGCCTCGACGAACTGCGTTTTAGCGGCGTAC  >  W3110S.gb/3440857‑3440987

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: