Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,923 G→A 28.8% E55K (GAG→AAG) ‡ ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9230GA28.8% 115.1 / 43.4 66E55K (GAG→AAG) ‡ycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (36/11);  new base A (6/13);  total (42/24)
Fisher's exact test for biased strand distribution p-value = 1.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAA  >  minE/910857‑910990
                                                                  |                                                                   
aTCGTGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGgaga                                                                  <  1:1087574/70‑1 (MQ=255)
 tCGTGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATc                                                                <  1:5374022/71‑1 (MQ=255)
   gtgtTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAg                                                              >  1:1155616/1‑71 (MQ=255)
   gtgtTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAg                                                              <  1:4314524/71‑1 (MQ=255)
   gtgtTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCac                                                              >  1:888846/1‑70 (MQ=9)
     gtTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCac                                                              >  1:3636701/1‑68 (MQ=2)
      tttGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATc                                                                <  1:394054/66‑1 (MQ=2)
       ttGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAgcgg                                                           >  1:404954/1‑70 (MQ=255)
        tGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAgcggcg                                                         >  1:3345148/1‑71 (MQ=255)
               ccGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCt                                                  >  1:216732/1‑71 (MQ=2)
                cGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTg                                                 <  1:1656270/71‑1 (MQ=9)
                cGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTg                                                 <  1:127381/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:811305/1‑70 (MQ=255)
                    gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                              <  1:5055415/70‑1 (MQ=255)
                    gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGa                                             >  1:1129392/1‑71 (MQ=255)
                    gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGa                                             >  1:2278525/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:4028421/1‑71 (MQ=255)
                                   gACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:1960426/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:2493386/1‑46 (MQ=255)
                                          cTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGt                                                >  1:5182558/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:5325610/1‑52 (MQ=255)
                                           tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGa                                         >  1:4735850/1‑52 (MQ=255)
                                           tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGa                                         >  1:4463217/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:5464613/1‑71 (MQ=255)
                                                   aTTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGc              >  1:2111343/1‑71 (MQ=255)
                                                    ttGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGc              >  1:3299400/1‑70 (MQ=255)
                                                    ttGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGc              >  1:5085129/1‑70 (MQ=255)
                                                     tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaa                              >  1:2537131/1‑53 (MQ=255)
                                                     tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGcc             >  1:3397989/1‑70 (MQ=255)
                                                     tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGcc             >  1:5360278/1‑70 (MQ=255)
                                                     tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGcc             >  1:3391508/1‑70 (MQ=255)
                                                      gCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatg                            <  1:2574722/54‑1 (MQ=255)
                                                          ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatgaat                         >  1:4148776/1‑53 (MQ=255)
                                                          ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatgaat                         >  1:3226819/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:4989657/52‑1 (MQ=255)
                                                                  gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                  <  1:2626302/52‑1 (MQ=255)
                                                                  |                                                                   
ATCGTGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAA  >  minE/910857‑910990

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: