Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 827,597 C→T 82.1% P26P (CCG→CCA ybhS ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb827,5970CT82.1% 86.2 / 14.3 39P26P (CCG→CCAybhSpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base T (16/16);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGG  >  W3110S.gb/827533‑827660
                                                                |                                                               
tAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                             >  1:142238/1‑69 (MQ=255)
tAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                             >  1:2282982/1‑69 (MQ=255)
 aaTGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATc                                                            >  1:354442/1‑69 (MQ=255)
 aaTGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATc                                                            >  1:351852/1‑69 (MQ=255)
   tGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATcg                                                           >  1:2592407/1‑68 (MQ=255)
     ccGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATcgcg                                                         >  1:2290631/1‑68 (MQ=255)
     ccGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGa                                                        >  1:994775/1‑69 (MQ=255)
     ccGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGa                                                        >  1:1759235/1‑69 (MQ=255)
     ccGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGa                                                        >  1:2313667/1‑69 (MQ=255)
     ccGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGACAGCTACTTGGATCGCGa                                                        >  1:144086/1‑69 (MQ=255)
       gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAAc                                                      <  1:1439465/69‑1 (MQ=255)
          aCCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGa                                                    <  1:529313/68‑1 (MQ=255)
                    aaCAGGAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGtc                                         >  1:698081/1‑69 (MQ=255)
                         cagTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGt                                          >  1:1910333/1‑63 (MQ=255)
                         cagTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCttt                                    <  1:1240602/69‑1 (MQ=255)
                         cagTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCttt                                    >  1:312663/1‑69 (MQ=255)
                           gTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTaa                                  <  1:520224/69‑1 (MQ=255)
                            tAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTaa                                  <  1:69480/68‑1 (MQ=255)
                             agcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGa                                                    <  1:1087398/49‑1 (MQ=255)
                             agcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGa                                                    <  1:511295/49‑1 (MQ=255)
                             agcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACg                                <  1:2679172/69‑1 (MQ=255)
                                agcGGGATCACTACCGCGATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                              <  1:2610355/68‑1 (MQ=255)
                                agcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGcac                             <  1:3446070/69‑1 (MQ=255)
                                agcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                              <  1:783591/68‑1 (MQ=255)
                                    ggATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAgcg                         >  1:2390593/1‑69 (MQ=255)
                                       tCACTACCGCAATCAGCCAGCTACTTGGATCGCGAAc                                                      >  1:461553/1‑37 (MQ=255)
                                       tCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACg                                                     >  1:3637927/1‑38 (MQ=255)
                                       tCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACg                                                     >  1:2104048/1‑38 (MQ=255)
                                            aCCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAc                  >  1:268426/1‑68 (MQ=255)
                                               gCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                <  1:587932/67‑1 (MQ=255)
                                                   tCAGCCAGCTACTTGGATCGCGAACGATCTGCGGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGcc          <  1:972728/69‑1 (MQ=255)
                                                    cagccagcTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACg                                <  1:2717081/46‑1 (MQ=255)
                                                     agccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcag        >  1:1557295/1‑69 (MQ=255)
                                                     agccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcag        <  1:1481710/69‑1 (MQ=255)
                                                     agccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcag        <  1:1477478/69‑1 (MQ=255)
                                                        cagcTACTTGGATCGCGAACGATCTGCCGCGGCTCTTTAACGCACAGCGCCCGTACGCGACGCcaggac     <  1:876528/69‑1 (MQ=255)
                                                          gcTACTTGGATCGCGAAGGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcaggacag   >  1:725489/1‑69 (MQ=255)
                                                          gcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcaggacag   >  1:2140620/1‑69 (MQ=255)
                                                           cTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcaggacagg  <  1:1652243/69‑1 (MQ=255)
                                                                |                                                               
TAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGG  >  W3110S.gb/827533‑827660

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: