Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,946 G→C 10.8% L62F (TTG→TTC ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9460GC10.8% 200.5 / 12.3 74L62F (TTG→TTCycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (38/28);  new base C (2/6);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTACG  >  minE/910877‑911016
                                                                      |                                                                      
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)
   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)
                      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)
                              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:5360278/1‑70 (MQ=255)
                                 tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGcc                                       >  1:3397989/1‑70 (MQ=255)
                                 tGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGcc                                       >  1:3391508/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:3291440/66‑1 (MQ=255)
                                      ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCa                                      <  1:5373925/66‑1 (MQ=255)
                                      ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCa                                 >  1:1218710/1‑71 (MQ=255)
                                         ttGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                                          <  1:2853669/43‑1 (MQ=255)
                                         ttGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                                          <  1:446288/43‑1 (MQ=255)
                                          tGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGc                              >  1:4399760/1‑70 (MQ=255)
                                            aggagATCAGCGGCGGCGATCTGTTGATCGACTGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa                            <  1:1740306/70‑1 (MQ=255)
                                            aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCa                                          >  1:2164008/1‑56 (MQ=255)
                                            aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCa                                          >  1:2190602/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)
                                                       ggcggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGa                    <  1:2257622/67‑1 (MQ=255)
                                                       ggcggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGa                    <  1:3307363/67‑1 (MQ=255)
                                                          ggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACg                               >  1:4398509/1‑53 (MQ=255)
                                                          ggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGtt             >  1:4005116/1‑71 (MQ=255)
                                                           gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa                            <  1:3354471/55‑1 (MQ=255)
                                                           gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa                            <  1:1407766/55‑1 (MQ=255)
                                                           gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa                            <  1:3359202/55‑1 (MQ=255)
                                                           gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc            <  1:158861/71‑1 (MQ=255)
                                                           ccgATCTTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAg                                   >  1:2699680/2‑49 (MQ=255)
                                                              aTCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa                             <  1:2900474/51‑1 (MQ=255)
                                                                  tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGc                              <  1:5012083/47‑1 (MQ=255)
                                                                  tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGc                      <  1:802163/55‑1 (MQ=255)
                                                                  tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGc                      <  1:1088744/55‑1 (MQ=255)
                                                                  tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGc                      <  1:4241606/55‑1 (MQ=255)
                                                                  tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGc                      <  1:2534246/55‑1 (MQ=255)
                                                                  tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGa                    >  1:2297004/1‑57 (MQ=255)
                                                                   gTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAAc      >  1:492713/1‑70 (MQ=255)
                                                                   gTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAAc      >  1:4642395/1‑70 (MQ=255)
                                                                   gTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAAc      >  1:890858/1‑70 (MQ=255)
                                                                   gTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTCTTCCAGAAc      >  1:1172287/1‑70 (MQ=255)
                                                                    ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGAtt                  <  1:793852/57‑2 (MQ=255)
                                                                    ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAg         <  1:566850/66‑1 (MQ=255)
                                                                    ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAAc      <  1:4953420/69‑1 (MQ=255)
                                                                    ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAAc      <  1:1640297/69‑1 (MQ=255)
                                                                    ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTa    >  1:2992337/1‑71 (MQ=255)
                                                                      gATCGACGGCAAACGAATGAATGACGTTCCCGCCa                                      >  1:648233/1‑35 (MQ=255)
                                                                      gATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAata                          >  1:55859/1‑47 (MQ=255)
                                                                      gATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTAcg  >  1:947769/1‑71 (MQ=255)
                                                                      |                                                                      
GGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTACG  >  minE/910877‑911016

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: