Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,442,070 A→G 18.4% T497A (ACG→GCG)  thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,442,0700AG18.4% 106.1 / 19.4 49T497A (ACG→GCG) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base G (2/7);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

ACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATT  >  W3110S.gb/3442002‑3442138
                                                                    |                                                                    
aCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                                                                    >  1:1135255/1‑71 (MQ=255)
      tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                >  1:194183/1‑69 (MQ=255)
       cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                >  1:2425015/1‑68 (MQ=255)
       cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                >  1:3422977/1‑68 (MQ=255)
       cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                >  1:3487759/1‑68 (MQ=255)
       cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                >  1:3706619/1‑68 (MQ=255)
       cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                >  1:3758739/1‑68 (MQ=255)
       cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                                >  1:4031233/1‑68 (MQ=255)
          gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc                                                          <  1:1495588/71‑1 (MQ=255)
          gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc                                                          <  1:2464096/71‑1 (MQ=255)
           tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCa                                                         <  1:3371667/71‑1 (MQ=255)
           tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCa                                                         <  1:2957858/71‑1 (MQ=255)
                gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTg                                                     >  1:2507763/1‑70 (MQ=255)
                 aTTGGTGCGGCGATGATTGGCTGGTTTGGTTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTggg                                                   <  1:3788940/71‑1 (MQ=255)
                 aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTggg                                                   >  1:3732918/1‑71 (MQ=255)
                 aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTggg                                                   <  1:1366099/71‑1 (MQ=255)
                 aTTGGTGCGGCGATGATTGGCTGGTTTGGCT‑‑GCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTc                                                 >  1:1249628/1‑71 (MQ=255)
                    ggTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCt                                                >  1:1400025/1‑71 (MQ=255)
                    ggTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCt                                                >  1:915467/1‑71 (MQ=255)
                             atgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                        >  1:2094733/1‑70 (MQ=255)
                             atgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                        >  1:2009242/1‑70 (MQ=255)
                             atgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                        >  1:2004602/1‑70 (MQ=255)
                                 ttGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                   <  1:398588/71‑1 (MQ=255)
                                  tGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCAt                                                        >  1:2644865/1‑49 (MQ=255)
                                     cTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATg                               <  1:270778/71‑1 (MQ=255)
                                        gTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                                           <  1:1389619/40‑1 (MQ=255)
                                        gTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                                           <  1:3866531/40‑1 (MQ=255)
                                        gTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                                           <  1:2178243/40‑1 (MQ=255)
                                        gTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                                           <  1:1859788/40‑1 (MQ=255)
                                            ggCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc                         <  1:1622368/70‑1 (MQ=255)
                                            ggCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc                         <  1:626410/70‑1 (MQ=255)
                                                gcgcGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc                         <  1:1834339/66‑1 (MQ=255)
                                                      tGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                                            <  1:565884/41‑1 (MQ=255)
                                                         tctGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt                             <  1:820362/53‑1 (MQ=255)
                                                         tctGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt                             <  1:756229/53‑1 (MQ=255)
                                                             ttACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtata       <  1:1095693/71‑1 (MQ=255)
                                                             ttACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaag                                  <  1:2898300/44‑1 (MQ=255)
                                                             ttACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                   <  1:4004442/43‑1 (MQ=255)
                                                              tACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtatt       <  1:2791400/70‑2 (MQ=255)
                                                               aCGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                                 >  1:261041/1‑43 (MQ=255)
                                                               aCGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                                 >  1:1423003/1‑43 (MQ=255)
                                                               aCGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                                 >  1:2056239/1‑43 (MQ=255)
                                                               aCGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtat        >  1:4045430/1‑68 (MQ=255)
                                                                  tAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                            <  1:1475069/45‑1 (MQ=255)
                                                                  tAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATaa      <  1:2574235/67‑1 (MQ=255)
                                                                  tAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATaa      <  1:1978463/67‑1 (MQ=255)
                                                                  tAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt  <  1:2709178/71‑1 (MQ=255)
                                                                  tAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt  <  1:4041536/71‑1 (MQ=255)
                                                                  tAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt  <  1:502046/71‑1 (MQ=255)
                                                                    |                                                                    
ACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATT  >  W3110S.gb/3442002‑3442138

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: