Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,775,080 T→C 62.9% F285F (TTT→TTC ydiN → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,775,0800TC62.9% 54.5 / 63.8 62F285F (TTT→TTCydiNpredicted transporter
Reads supporting (aligned to +/- strand):  ref base T (5/18);  new base C (17/22);  total (22/40)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

TGTCTTTATTTTTGCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACTGGTGTGCAA  >  W3110S.gb/1775020‑1775145
                                                            |                                                                 
tgtCTTTATTTTTGCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTg                                                            <  1:2628881/68‑1 (MQ=255)
tgtCTTTATTTTTGCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTg                                                            <  1:316845/68‑1 (MQ=255)
          ttagCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGcaacaa                                                <  1:575571/67‑1 (MQ=255)
           ttGCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAAt                                               <  1:315190/70‑1 (MQ=255)
           ttGCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAAt                                               <  1:2453133/70‑1 (MQ=255)
             gCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATaa                                             <  1:2891280/70‑1 (MQ=255)
             gCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATaa                                             <  1:545584/70‑1 (MQ=255)
                   ctactGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGcaacaa                                                <  1:1019559/61‑1 (MQ=255)
                   ctactGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGcaacaa                                                <  1:2946720/61‑1 (MQ=255)
                   ctactGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGcaacaa                                                <  1:2875472/61‑1 (MQ=255)
                   ctactGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGcaacaa                                                <  1:2528973/61‑1 (MQ=255)
                   ctactGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGcaacaa                                                <  1:1077711/61‑1 (MQ=255)
                   ctactGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACagcagc                                      <  1:2829033/71‑1 (MQ=255)
                   ctactGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACagcagc                                      >  1:73693/1‑71 (MQ=255)
                   ctactGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACagcagc                                      <  1:3359442/71‑1 (MQ=255)
                     actGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATa                                              >  1:1498527/1‑61 (MQ=255)
                     actGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCaa                                                                 >  1:2216894/1‑42 (MQ=255)
                     actGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATaa                                             >  1:288804/1‑62 (MQ=255)
                      ctGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCAt                                   >  1:3675817/1‑71 (MQ=255)
                         aaaaaaaTGGTCCGGCCCATCTGGGCTAATGTATTCAActc                                                              >  1:794664/1‑41 (MQ=255)
                         aaaaaaaTGGTCCGGCCCATCTGGGCTAATGTATTCAActc                                                              >  1:405310/1‑41 (MQ=255)
                         aaaaaaaTGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:2027431/69‑1 (MQ=255)
                         aaaaaaaTGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:1650988/69‑1 (MQ=255)
                           aaaaaTGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTTGCAACAATAACAGCAGCCATTATCt                              <  1:1629176/71‑1 (MQ=255)
                           aaaaaTGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCt                              <  1:6467/71‑1 (MQ=255)
                           aaaaaTGGTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCt                              <  1:3201797/71‑1 (MQ=255)
                                 ggTCCGGCCCATCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCt                              <  1:2577405/65‑1 (MQ=255)
                                            tCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:2153422/50‑1 (MQ=255)
                                            tCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:2149247/50‑1 (MQ=255)
                                            tCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:1721801/50‑1 (MQ=255)
                                            tCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:2850391/50‑1 (MQ=255)
                                            tCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:3700146/50‑1 (MQ=255)
                                            tCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:2815976/50‑1 (MQ=255)
                                            tCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:3537718/50‑1 (MQ=255)
                                            tCTGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCAtt                                  <  1:334523/50‑1 (MQ=255)
                                              tGGGCTAATGTATTCAACTCTGCACTGGCAACAATaa                                             >  1:1521317/1‑37 (MQ=255)
                                              tGGGCTAATGTATTCAACTCTGCACTGGCAACAATa                                              >  1:1464864/1‑36 (MQ=255)
                                              tGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACt           >  1:2992422/1‑71 (MQ=255)
                                              tGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACt           >  1:2287028/1‑71 (MQ=255)
                                              tGGGCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACt           >  1:214261/1‑71 (MQ=255)
                                              agggCTAATGTATTCAACTCTGCACTGGCAACAATaa                                             >  1:1197468/2‑37 (MQ=38)
                                                ggCTAATGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACTgg         <  1:2805637/71‑1 (MQ=255)
                                                    aaTGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTa                                 >  1:115113/1‑43 (MQ=255)
                                                    aaTGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATc                               >  1:2799939/1‑45 (MQ=255)
                                                    aaTGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATc                               >  1:1080418/1‑45 (MQ=255)
                                                    aaTGTACTCAACTCTGCACTGGCAACAATAACAGCCGCCAtt                                  <  1:3055688/42‑1 (MQ=37)
                                                     aTGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACTGGTGTGc    <  1:2284257/71‑1 (MQ=255)
                                                      tGTATTTAACTCTGCACTGGCAACAATAACagcagc                                      <  1:2734866/36‑1 (MQ=255)
                                                      tGTATTTAACTCTGCACTGGCAACAATAACagcagc                                      <  1:1429207/36‑1 (MQ=255)
                                                      tGTATTCAACTCTGCACTGGCAACAATAACAGCAGcc                                     <  1:1051884/37‑1 (MQ=255)
                                                      tGTATTCAACTCTGCACTGGCAACAATAACAGCAGcc                                     <  1:2521815/37‑1 (MQ=255)
                                                      tGTATTCAACTCTGCACTGGCAACAATAACAGCAGcc                                     <  1:1484565/37‑1 (MQ=255)
                                                      tGTATTCAACTCTGCACTGGCAACAATAACAGCAGcc                                     <  1:1422639/37‑1 (MQ=255)
                                                      tGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCt                              <  1:3133355/44‑1 (MQ=255)
                                                      tGTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCt                              <  1:822292/44‑1 (MQ=255)
                                                       gTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACTGGTGTGCaa  >  1:2830153/1‑71 (MQ=255)
                                                       gTATTCAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACTGGTGTGCaa  >  1:468684/1‑71 (MQ=255)
                                                        tATTTAACTCTGCACTGGCAATAATAACAGCAGCCa                                    >  1:1297513/1‑36 (MQ=255)
                                                        tATTTAACTCTGCACTGGCAACAATAACAGCAGcc                                     >  1:903345/1‑35 (MQ=255)
                                                        tATTTAACTCTGCACTGGCAACAATAACAGCAGCCa                                    >  1:1328206/1‑36 (MQ=255)
                                                        tATTTAACTCTGCACTGGCAACAATAACAGCAGCCa                                    >  1:1075759/1‑36 (MQ=255)
                                                        tATTTAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTtc                 <  1:842210/55‑1 (MQ=255)
                                                            |                                                                 
TGTCTTTATTTTTGCCGCACTACTGAAAAAAATGGTCCGGCCCATCTGGGCTAATGTATTTAACTCTGCACTGGCAACAATAACAGCAGCCATTATCTACCTGTACCCTTCTCCACTGGTGTGCAA  >  W3110S.gb/1775020‑1775145

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: