Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 464,946 A→G 34.1% F531L (TTC→CTC)  ybaE ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb464,9460AG34.1% 38.1 / 38.4 41F531L (TTC→CTC) ybaEpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (12/15);  new base G (8/6);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

ACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCT  >  W3110S.gb/464877‑465016
                                                                     |                                                                      
aCCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGaa                                                                        >  1:675115/1‑70 (MQ=255)
    gCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAgcg                                                                   >  1:721285/1‑71 (MQ=255)
      cgcgcgGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAgcggc                                                                 >  1:1025417/1‑71 (MQ=255)
         gcgGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTc                                                              <  1:3561936/71‑1 (MQ=255)
         gcgGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGt                                                               >  1:160229/1‑70 (MQ=255)
                     gCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCt                                                  >  1:3681670/1‑71 (MQ=255)
                          ccGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCAt                                             >  1:565963/1‑71 (MQ=255)
                          ccGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCAt                                             >  1:1293347/1‑71 (MQ=255)
                             ttCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTaa                                          >  1:4098264/1‑71 (MQ=255)
                              tCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTc                                                        <  1:186623/56‑1 (MQ=255)
                                 cGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTc                                                              <  1:1804989/47‑1 (MQ=255)
                                  gCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATc                                                   >  1:1566494/1‑57 (MQ=255)
                                  gCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATc                                                   >  1:3159888/1‑57 (MQ=255)
                                  gCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATc                                                   >  1:2663462/1‑57 (MQ=255)
                                    cGGGAGTGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGTCTa                                   <  1:4608075/71‑1 (MQ=255)
                                      ggAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGc                                                           <  1:1271190/45‑1 (MQ=255)
                                      ggAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGc                                                           <  1:3359580/45‑1 (MQ=255)
                                      ggAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTTGCTaaa                                 <  1:656068/71‑1 (MQ=255)
                                      ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                                 <  1:1387529/71‑1 (MQ=255)
                                       gAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAg                                                            >  1:1375460/1‑43 (MQ=255)
                                        aGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaaaa                               >  1:1954426/1‑71 (MQ=255)
                                          gggCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGc                                                      <  1:297255/46‑1 (MQ=255)
                                               cTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCttt                        <  1:1532984/71‑1 (MQ=255)
                                               cTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                                 <  1:1174603/62‑1 (MQ=255)
                                               cTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                                 <  1:1197069/62‑1 (MQ=255)
                                               cTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                                 <  1:4480148/62‑1 (MQ=255)
                                               cTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                                 <  1:555259/62‑1 (MQ=255)
                                                 aaTGCGATAGTGATAGTTGCGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCttt                        >  1:1837782/1‑69 (MQ=255)
                                                 aaTGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCttt                        >  1:1877017/1‑69 (MQ=255)
                                                 aaTGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCttt                        >  1:4652901/1‑69 (MQ=255)
                                                 aaTGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCttt                        >  1:4776180/1‑69 (MQ=255)
                                                  aTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTc                       <  1:3486390/69‑1 (MQ=255)
                                                  aTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTc                       <  1:1432715/69‑1 (MQ=255)
                                                  aTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCa                      <  1:2947177/70‑1 (MQ=255)
                                                   tGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaaa                                >  1:1128610/1‑59 (MQ=255)
                                                   tGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCAGTCATTAAc                                         >  1:3860712/1‑50 (MQ=255)
                                                     cGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAggg                   >  1:2582775/1‑70 (MQ=255)
                                                     cGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAggg                   >  1:429785/1‑70 (MQ=255)
                                                     cGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAggg                   <  1:422225/70‑1 (MQ=255)
                                                             atagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATc          <  1:112908/71‑1 (MQ=255)
                                                                gtTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCgg        <  1:5056301/70‑1 (MQ=255)
                                                                   gAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGttt     <  1:4665469/70‑1 (MQ=255)
                                                                     aCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAg                     >  1:1958167/1‑52 (MQ=255)
                                                                     aCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTttct  <  1:4533826/71‑1 (MQ=255)
                                                                     |                                                                      
ACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCT  >  W3110S.gb/464877‑465016

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: