Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,930 G→C 24.6% S57T (AGC→ACC)  ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9300GC24.6% 126.6 / 36.2 65S57T (AGC→ACC) ycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (35/14);  new base C (4/12);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 2.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGA  >  minE/910860‑910998
                                                                      |                                                                    
gtgtTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAg                                                                      >  1:1155616/1‑71 (MQ=255)
gtgtTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAg                                                                      <  1:4314524/71‑1 (MQ=255)
    ttGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAgcgg                                                                   >  1:404954/1‑70 (MQ=255)
     tGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAgcggcg                                                                 >  1:3345148/1‑71 (MQ=255)
            ccGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCt                                                          >  1:216732/1‑71 (MQ=2)
             cGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTg                                                         <  1:127381/71‑1 (MQ=9)
             cGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTg                                                         <  1:1656270/71‑1 (MQ=9)
             cGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTg                                                         <  1:120551/71‑1 (MQ=9)
               tctggCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGtt                                                       <  1:281672/68‑1 (MQ=9)
               tCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGtt                                                       >  1:5202624/1‑71 (MQ=255)
               tCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGtt                                                       >  1:4159067/1‑71 (MQ=255)
                 tggCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGttca                                                     >  1:389365/2‑69 (MQ=16)
                 gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGt                                                        >  1:1951081/1‑68 (MQ=255)
                 gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                      <  1:5055415/70‑1 (MQ=255)
                 gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                      >  1:811305/1‑70 (MQ=255)
                 gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGa                                                     >  1:2278525/1‑71 (MQ=255)
                 gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGa                                                     >  1:1129392/1‑71 (MQ=255)
                   gCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGttca                                                     <  1:1958908/69‑3 (MQ=16)
                    cTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCg                                                  >  1:19106/1‑71 (MQ=255)
                       cGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCggt                                                <  1:245155/70‑3 (MQ=255)
                          tAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCa                                            <  1:3368712/71‑1 (MQ=255)
                           aaGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCa                                            >  1:747362/1‑70 (MQ=255)
                            aGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCaaa                                          >  1:3327204/1‑71 (MQ=255)
                                gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                                      >  1:1960426/1‑71 (MQ=255)
                                gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                                      >  1:4028421/1‑71 (MQ=255)
                                 accaccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                      <  1:315964/54‑1 (MQ=255)
                                   caccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatg                                    >  1:2835923/1‑70 (MQ=255)
                                     tacTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACggatgaat                                 >  1:4872102/3‑71 (MQ=255)
                                      acTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACggatgaat                                 <  1:617487/69‑1 (MQ=255)
                                       cTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGt                                                        >  1:5182558/1‑46 (MQ=255)
                                       cTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGt                                                        >  1:2493386/1‑46 (MQ=255)
                                       cTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGAt                                                    >  1:3697615/1‑50 (MQ=255)
                                       cTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACggatgaat                                 <  1:4600486/69‑1 (MQ=255)
                                       cTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGatgaatga                               <  1:357269/71‑1 (MQ=255)
                                        tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGa                                                 >  1:4463217/1‑52 (MQ=255)
                                        tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGa                                                 >  1:4735850/1‑52 (MQ=255)
                                        tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGa                                                 >  1:5325610/1‑52 (MQ=255)
                                        tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCaaa                                          >  1:1088077/1‑59 (MQ=255)
                                          cgcATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGt                            >  1:3574296/1‑71 (MQ=255)
                                          cgcATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGatgaatga                               <  1:2731504/68‑1 (MQ=255)
                                           gcATGATTGCCGGGCTTGAGGAGATCAGCGGAGGCGATCTGTTGATCGACGGCAAACgaatgaat                                 <  1:233013/65‑1 (MQ=255)
                                             aTGATTGCCGGGCTTGAGACGATCACCAGCGGCGACATGtt                                                       <  1:1951221/41‑1 (MQ=12)
                                               gATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg                       >  1:2038626/1‑71 (MQ=255)
                                               gATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAg                       >  1:4874606/1‑71 (MQ=255)
                                               gATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTTTTCATCggt                                                <  1:2110753/46‑3 (MQ=255)
                                                aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGc                      >  1:2111343/1‑71 (MQ=255)
                                                aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGc                      >  1:5464613/1‑71 (MQ=255)
                                                 ttGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGc                      >  1:5085129/1‑70 (MQ=255)
                                                 ttGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGc                      >  1:3299400/1‑70 (MQ=255)
                                                  tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                                      >  1:2537131/1‑53 (MQ=255)
                                                  tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGcc                     >  1:5360278/1‑70 (MQ=255)
                                                  tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGcc                     >  1:3391508/1‑70 (MQ=255)
                                                  tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGcc                     >  1:3397989/1‑70 (MQ=255)
                                                   gCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatg                                    <  1:2574722/54‑1 (MQ=255)
                                                       ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatgaat                                 >  1:3226819/1‑53 (MQ=255)
                                                       ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatgaat                                 >  1:4148776/1‑53 (MQ=255)
                                                       ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCa                    <  1:5373925/66‑1 (MQ=255)
                                                       ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCa                    <  1:3291440/66‑1 (MQ=255)
                                                       ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCa               >  1:1218710/1‑71 (MQ=255)
                                                          ttGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                        <  1:446288/43‑1 (MQ=255)
                                                          ttGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                        <  1:2853669/43‑1 (MQ=255)
                                                           tGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGc            >  1:4399760/1‑70 (MQ=255)
                                                             aggagATCAGCGGCGGCGATCTGTTGATCGACTGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa          <  1:1740306/70‑1 (MQ=255)
                                                             aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCa                        >  1:2190602/1‑56 (MQ=255)
                                                             aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCa                        >  1:2164008/1‑56 (MQ=255)
                                                               gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                          <  1:2626302/52‑1 (MQ=255)
                                                               gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                          <  1:4989657/52‑1 (MQ=255)
                                                                  aTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAg                 <  1:5521638/58‑1 (MQ=255)
                                                                     aGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGa  <  1:3411660/70‑1 (MQ=255)
                                                                      |                                                                    
GTGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGA  >  minE/910860‑910998

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: