Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,164,775 C→T 63.2% A534T (GCA→ACA)  ygiS ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,164,7750CT63.2% 46.4 / 60.0 57A534T (GCA→ACA) ygiSpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (21/15);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGG  >  W3110S.gb/3164707‑3164843
                                                                    |                                                                    
gAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCt                                                                    >  1:50708/1‑71 (MQ=255)
gAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCt                                                                    >  1:269150/1‑71 (MQ=255)
  gTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTg                                                                  <  1:2992823/71‑1 (MQ=255)
     acaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGata                                                               >  1:2565687/1‑71 (MQ=255)
        ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGAtata                                                             >  1:1085939/1‑70 (MQ=255)
        ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGAtata                                                             >  1:868002/1‑70 (MQ=255)
        ggAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGAtata                                                             >  1:1833404/1‑70 (MQ=255)
              cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGata                                                               >  1:997464/1‑62 (MQ=255)
              cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAActct                                                      >  1:2342623/1‑71 (MQ=255)
              cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAActct                                                      >  1:2672024/1‑71 (MQ=255)
              cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAActct                                                      >  1:2476120/1‑71 (MQ=255)
              cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAActct                                                      >  1:2016528/1‑71 (MQ=255)
               aaCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCtt                                                     >  1:1661633/1‑71 (MQ=255)
               aaCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCtt                                                     >  1:603904/1‑71 (MQ=255)
                     aaTAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTg                                                >  1:1629005/1‑70 (MQ=255)
                         gcgcCGCCGGACGTGCTCTGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTa                                              >  1:594048/1‑68 (MQ=255)
                         gcgcCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTa                                              >  1:1540557/1‑68 (MQ=255)
                         gcgcCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGa                                            <  1:2247740/70‑1 (MQ=255)
                                cgGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATcc                                    <  1:417538/71‑1 (MQ=255)
                                cgGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATcc                                    <  1:432149/71‑1 (MQ=255)
                                 ggACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATcc                                    <  1:1600827/70‑1 (MQ=255)
                                     gTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATcc                                    <  1:1560265/66‑1 (MQ=255)
                                      tGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATcc                                    <  1:32097/65‑1 (MQ=255)
                                      tGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATcc                                    <  1:2382981/65‑1 (MQ=255)
                                        cTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAt                             <  1:1283400/70‑1 (MQ=255)
                                        cTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAt                             <  1:1367154/70‑1 (MQ=255)
                                        cTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAt                             <  1:109224/70‑1 (MQ=255)
                                          ccGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTAt                          >  1:455126/1‑71 (MQ=255)
                                           cGACGGCGATAAAGGCATCAATCTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATg                         <  1:2535040/71‑1 (MQ=255)
                                                   aTAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaa                 <  1:2949479/71‑1 (MQ=255)
                                                   aTAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaa                 <  1:2945259/71‑1 (MQ=255)
                                                    tAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATcc                                    <  1:1033919/51‑1 (MQ=255)
                                                    tAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa                <  1:661040/71‑1 (MQ=255)
                                                    tAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa                <  1:1630597/71‑1 (MQ=255)
                                                         gCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATcc                                    <  1:528632/46‑1 (MQ=255)
                                                            tCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATaa                                       <  1:1363151/40‑1 (MQ=255)
                                                            tCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATaa                                       <  1:193028/40‑1 (MQ=255)
                                                            tCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATaa                                       <  1:1605098/40‑1 (MQ=255)
                                                            tCAATGTGCCTTGATATACAACTCTTTGCTGTGACATaa                                       <  1:2652033/39‑1 (MQ=39)
                                                             cAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTTGGGATTATGCAGCGGAAAACCGCCAAc       >  1:71419/1‑71 (MQ=255)
                                                              aaTGTGTCTTGATATACAACTCTTTGCTGTAGACATAATc                                     >  1:153650/1‑40 (MQ=255)
                                                              aaTGTGTCTTGATATACAACTCTTTGCTGTAGACATAATc                                     >  1:2813673/1‑40 (MQ=255)
                                                              aaTGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCaa        >  1:3313290/1‑69 (MQ=255)
                                                              aaTGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACg      >  1:2711871/1‑71 (MQ=255)
                                                              aaTGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACg      >  1:1940043/1‑71 (MQ=255)
                                                                tgtgTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTa    <  1:1385741/71‑1 (MQ=255)
                                                                tgtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGa                              <  1:3097514/45‑1 (MQ=255)
                                                                tgtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGa                              <  1:1382024/45‑1 (MQ=255)
                                                                tgtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGa                              <  1:1821061/45‑1 (MQ=255)
                                                                tgtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGa                              <  1:1753986/45‑1 (MQ=255)
                                                                tgtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGa                              <  1:1748091/45‑1 (MQ=255)
                                                                  tgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAtt                            >  1:341274/1‑45 (MQ=255)
                                                                  tgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAtt                            >  1:3452445/1‑45 (MQ=255)
                                                                  tgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAtt                            >  1:35408/1‑45 (MQ=255)
                                                                  tgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAtt                            >  1:2982395/1‑45 (MQ=255)
                                                                  tgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAgg  >  1:3278657/1‑71 (MQ=255)
                                                                  tgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAgg  >  1:529599/1‑71 (MQ=255)
                                                                    |                                                                    
GAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGG  >  W3110S.gb/3164707‑3164843

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: