Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 72,433 A→G 8.6% R165R (CGT→CGC thiQ ← thiamin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb72,4330AG 8.6% 149.7 / 6.5 58R165R (CGT→CGCthiQthiamin transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (34/19);  new base G (2/3);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

TGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAAC  >  W3110S.gb/72365‑72503
                                                                    |                                                                      
tgCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACg                                                                      <  1:4223655/71‑1 (MQ=255)
      acCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGc                                                                >  1:3777371/1‑71 (MQ=255)
         aTCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCg                                                              >  1:3944458/1‑70 (MQ=255)
            aaTAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGGCGTAACGCCGGATc                                                          >  1:1120134/1‑71 (MQ=255)
            aaTAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATc                                                          >  1:3159930/1‑71 (MQ=255)
               aGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAg                                                       >  1:1216500/1‑71 (MQ=255)
                gCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAgc                                                      >  1:3496626/1‑71 (MQ=255)
                 cGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAgcg                                                     <  1:3557771/71‑1 (MQ=255)
                  gTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAgcgc                                                    <  1:3058750/71‑1 (MQ=255)
                   tCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCa                                                   <  1:4124806/71‑1 (MQ=255)
                    cATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAgcg                                                     >  1:156303/1‑68 (MQ=255)
                     aTTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTaa                                                                   >  1:1278225/1‑53 (MQ=255)
                     aTTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCa                                                   >  1:3621453/1‑69 (MQ=255)
                         tttGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCg                                                              >  1:394458/1‑54 (MQ=255)
                                        tCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCaa                               >  1:3869715/1‑70 (MQ=255)
                                        tCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCaa                               >  1:3788120/1‑70 (MQ=255)
                                          gTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATaa                            <  1:4085840/71‑1 (MQ=255)
                                             cTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATaaaa                          <  1:421986/70‑1 (MQ=255)
                                             cTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAAt                         >  1:3051519/1‑71 (MQ=255)
                                               cACCAGCGTCAACATCTCGTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATaaa                           >  1:4621806/1‑67 (MQ=255)
                                               cACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATaaa                           >  1:268221/1‑67 (MQ=255)
                                               cACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATaaa                           >  1:2096520/1‑67 (MQ=255)
                                                aCCAGCGTCAACATCTCCTGGCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCg                       <  1:2303211/70‑1 (MQ=255)
                                                aCCAGCGTCAACATCTCCTGGCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCg                       <  1:528611/70‑1 (MQ=255)
                                                aCCAGCGTCAACATCTCCTGGCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCg                       <  1:3890783/70‑1 (MQ=255)
                                                   aGCGTCAACATCTCCTGGCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATaaaa                          >  1:1975062/1‑64 (MQ=255)
                                                   aGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCagag                                                <  1:1550832/42‑1 (MQ=255)
                                                   aGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCagag                                                <  1:3250588/42‑1 (MQ=255)
                                                   aGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTg                   <  1:4050771/71‑1 (MQ=255)
                                                      gTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGACCGGTTCATCGAGCAATAAAATCGGCTGTTc                >  1:3560726/1‑71 (MQ=255)
                                                       tCAACATCTCCTGACGTAACGCCGGATCGAGCGCagag                                                <  1:2006627/38‑1 (MQ=255)
                                                       tCAACATCTCCTGACGTAACGCCGGATCGAGCGCagag                                                <  1:2155478/38‑1 (MQ=255)
                                                        cAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATaaa                           <  1:7702/58‑1 (MQ=255)
                                                         aaCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGtt                 >  1:3001641/1‑67 (MQ=255)
                                                         aaCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg             >  1:3393610/1‑71 (MQ=255)
                                                         aaCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg             >  1:2400230/1‑71 (MQ=255)
                                                         aaCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg             >  1:4529395/1‑71 (MQ=255)
                                                         aaCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg             <  1:3153430/71‑1 (MQ=255)
                                                         aaCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg             >  1:1575183/1‑71 (MQ=255)
                                                         aaCATCTCCTGACGTAACGCCGGATCGAGCCCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg             >  1:4256426/1‑71 (MQ=255)
                                                          aCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg             <  1:3525865/70‑1 (MQ=255)
                                                            aTCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGtt                                         >  1:1447705/1‑40 (MQ=255)
                                                            aTCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCAt                                      >  1:2075915/1‑43 (MQ=255)
                                                            aTCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTAc          >  1:3666307/1‑71 (MQ=255)
                                                              ctcCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg             >  1:2267941/1‑66 (MQ=255)
                                                              ctcCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGt            >  1:3304705/1‑67 (MQ=255)
                                                              ctcCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTAc          >  1:869486/1‑69 (MQ=255)
                                                               tcCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTAc          >  1:4124488/1‑68 (MQ=255)
                                                                 ttGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGa      >  1:1209971/2‑70 (MQ=255)
                                                                 cTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCa        <  1:1352313/68‑1 (MQ=255)
                                                                 cTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGAc     <  1:4425659/71‑1 (MQ=255)
                                                                 cTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGAc     <  1:4500690/71‑1 (MQ=255)
                                                                  tGACGTAACGCCGGATCGAGCGCAGAGAACGGTTc                                        <  1:3502740/35‑1 (MQ=255)
                                                                  tGACGTAACGCCGGATCGAGCGCAGAGAACGGTTc                                        <  1:3301700/35‑1 (MQ=255)
                                                                  tGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCTTACCAGAc     >  1:1929027/1‑70 (MQ=255)
                                                                  tGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGAc     >  1:3467945/1‑70 (MQ=255)
                                                                   gACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACaa   >  1:3807570/1‑71 (MQ=255)
                                                                   gACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACaa   >  1:1259765/1‑71 (MQ=255)
                                                                    aCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAAt                         >  1:1233738/1‑48 (MQ=255)
                                                                    aCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAAc  <  1:1241885/71‑1 (MQ=255)
                                                                    |                                                                      
TGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAAC  >  W3110S.gb/72365‑72503

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: