Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 47,828 A→G 85.5% V20V (GTA→GTG kefC → potassium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb47,8280AG85.5% 141.1 / 14.9 55V20V (GTA→GTGkefCpotassium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (2/6);  new base G (31/16);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 4.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

AGGCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGG  >  W3110S.gb/47760‑47893
                                                                    |                                                                 
aggCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTAcc                                                                 <  1:2848987/71‑1 (MQ=255)
 ggCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGcc                                                                 >  1:241478/1‑70 (MQ=255)
 ggCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTAccc                                                                <  1:2653961/71‑1 (MQ=255)
 ggCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTAccc                                                                <  1:2128146/71‑1 (MQ=255)
 ggCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTAccc                                                                <  1:1894405/71‑1 (MQ=255)
 ggCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTAcc                                                                 >  1:1284010/1‑70 (MQ=255)
  gCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTAccc                                                                <  1:1835819/70‑1 (MQ=255)
       tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGcc                                                                 <  1:2682654/64‑1 (MQ=255)
       tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGcc                                                                 <  1:798980/64‑1 (MQ=255)
       tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGcc                                                                 <  1:2066243/64‑1 (MQ=255)
        catGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         >  1:3329159/1‑71 (MQ=255)
        catGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         >  1:2390785/1‑71 (MQ=255)
         atGGGTAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGTAGCGCTGATTGTGCCCATTGCgg                                                         >  1:720102/1‑70 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGGCCATTGCgg                                                         >  1:3000510/1‑70 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         <  1:962950/70‑1 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         >  1:2981026/1‑70 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         >  1:471952/1‑70 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         >  1:2085280/1‑70 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         >  1:669505/1‑70 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         >  1:1593470/1‑70 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCg                                                          >  1:1441601/1‑69 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCg                                                          >  1:1270372/1‑69 (MQ=255)
         atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCg                                                          >  1:1316776/1‑69 (MQ=255)
            gATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTa                                                                   >  1:633071/1‑57 (MQ=255)
            gATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTa                                                                   >  1:3069379/1‑57 (MQ=255)
            gATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTa                                                                   >  1:2448861/1‑57 (MQ=255)
            gATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACg                                                     >  1:1716611/1‑71 (MQ=255)
                   aTACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         <  1:3460444/60‑1 (MQ=255)
                   aTACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCgg                                                         <  1:2999935/60‑1 (MQ=255)
                   aTACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGt                                              <  1:354044/71‑1 (MQ=255)
                   aTACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGt                                              <  1:290377/71‑1 (MQ=255)
                    tACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGcc                                                                 <  1:194102/51‑1 (MQ=255)
                    tACGCTGATTCAGGCGCTGATTTATATCGGTTCGGCAGCGCTGATTGTGcc                                                                 <  1:3719893/51‑1 (MQ=38)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                               >  1:3613651/1‑51 (MQ=255)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                               >  1:1913158/1‑51 (MQ=255)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                               >  1:2918747/1‑51 (MQ=255)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTg                                           >  1:923666/1‑71 (MQ=255)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTg                                           >  1:166287/1‑71 (MQ=255)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTg                                           >  1:1607227/1‑71 (MQ=255)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTg                                           >  1:2123630/1‑71 (MQ=255)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTg                                           >  1:2019876/1‑71 (MQ=255)
                      cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTg                                           >  1:1967706/1‑71 (MQ=255)
                        cTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTggg                                         <  1:2683061/71‑1 (MQ=255)
                        cTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTggg                                         <  1:2612758/71‑1 (MQ=255)
                        cTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTggg                                         <  1:2245335/71‑1 (MQ=255)
                             tCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCgg                                    >  1:1941088/1‑71 (MQ=255)
                              cAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCgg                                    >  1:472748/1‑70 (MQ=255)
                              cAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCgg                                    <  1:790274/70‑1 (MQ=255)
                               aGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGt                                   <  1:3016848/70‑1 (MQ=255)
                                    cTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTgg                             >  1:3167488/1‑71 (MQ=255)
                                    cTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTgg                             >  1:3627080/1‑71 (MQ=255)
                                     tGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTggg                                         <  1:1860492/58‑1 (MQ=255)
                                          tATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCt                       >  1:3245642/1‑71 (MQ=255)
                                          tATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCt                       >  1:2881862/1‑71 (MQ=255)
                                                 gTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGcc                <  1:99513/71‑1 (MQ=255)
                                                        aGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCa         >  1:1359080/1‑71 (MQ=255)
                                                        aGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCa         >  1:333733/1‑71 (MQ=255)
                                                               aTTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTgg  >  1:829829/1‑71 (MQ=255)
                                                                    |                                                                 
AGGCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGG  >  W3110S.gb/47760‑47893

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