Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,892,177 G→A 12.6% D278N (GAC→AAC)  ytfQ → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,892,1770GA12.6% 242.2 / 21.6 87D278N (GAC→AAC) ytfQpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (44/32);  new base A (6/5);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AAAGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTA  >  minE/2892108‑2892247
                                                                     |                                                                      
aaaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGa                                                                       >  1:610956/1‑71 (MQ=255)
 aaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCg                                                                        >  1:4627829/1‑69 (MQ=255)
 aaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGa                                                                       >  1:4278533/1‑70 (MQ=255)
 aaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGa                                                                       >  1:733270/1‑70 (MQ=255)
 aaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGa                                                                       >  1:3989326/1‑70 (MQ=255)
 aaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAc                                                                      >  1:350493/1‑71 (MQ=255)
 aaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAc                                                                      >  1:3285118/1‑71 (MQ=255)
 aaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAc                                                                      >  1:515443/1‑71 (MQ=255)
 aaGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAc                                                                      >  1:1643416/1‑71 (MQ=255)
   gCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAcg                                                                     >  1:4703405/1‑70 (MQ=255)
        gatggatgGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTg                                                                <  1:3489938/70‑1 (MQ=255)
         atggatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCAACGCGCTg                                                                >  1:2391769/1‑69 (MQ=255)
         atggatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCAACGCGCTg                                                                >  1:532468/1‑69 (MQ=255)
         atggatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCAACGCGCTg                                                                >  1:1769379/1‑69 (MQ=255)
         atggatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCAACGCGCTg                                                                >  1:313426/1‑69 (MQ=255)
         atggatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCAACGCGCTg                                                                >  1:1549189/1‑69 (MQ=255)
          tggatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTg                                                                >  1:807567/1‑68 (MQ=255)
           ggatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAAGTCCCGCCTTCGACGCGCTGgag                                                             >  1:2969867/1‑70 (MQ=255)
            gatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGgaga                                                            <  1:2898732/70‑1 (MQ=255)
            gatggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGgaga                                                            <  1:347481/70‑1 (MQ=255)
              tggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAcgc                                                                    <  1:3789573/60‑1 (MQ=255)
              tggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAcgc                                                                    <  1:4449095/60‑1 (MQ=255)
              tggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAcgc                                                                    <  1:230630/60‑1 (MQ=255)
              tggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGAcgc                                                                    <  1:2170305/60‑1 (MQ=255)
              tggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGaaa                                                          >  1:133106/1‑70 (MQ=255)
              tggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGaaa                                                          >  1:2293657/1‑70 (MQ=255)
              tggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGaaa                                                          >  1:2665994/1‑70 (MQ=255)
              tggCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGaaa                                                          >  1:543985/1‑70 (MQ=255)
                gcgaagcgaaCGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATa                                                        >  1:3042404/1‑70 (MQ=255)
                gcgaagcgaaCGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATAc                                                       <  1:2407072/71‑1 (MQ=255)
                   aagcgaaCGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCAACGCGCTGGAGAAATACaaa                                                    <  1:4312439/71‑1 (MQ=255)
                     gcgaaCGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACaaaa                                                   >  1:557110/1‑70 (MQ=255)
                     gcgaaCGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACaaaa                                                   >  1:747872/1‑70 (MQ=255)
                      cgaaCGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACaaaaaa                                                 <  1:3338313/71‑1 (MQ=255)
                       gaaCGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGaa                                                           <  1:2465986/60‑1 (MQ=255)
                        aaCGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGaaa                                                          >  1:4677351/1‑60 (MQ=255)
                           gCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACgg                                            <  1:1750154/71‑1 (MQ=255)
                           gCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACgg                                            <  1:1577527/71‑1 (MQ=255)
                            ccAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAg                                                <  1:3439561/66‑1 (MQ=255)
                            ccAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAg                                                <  1:4129331/66‑1 (MQ=255)
                            ccAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAg                                                <  1:4722187/66‑1 (MQ=255)
                            ccAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAg                                                <  1:166448/66‑1 (MQ=255)
                               gtgtTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCAc                                         >  1:3581391/1‑70 (MQ=255)
                                  ttGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACaaaaaa                                                 <  1:4207962/59‑1 (MQ=255)
                                    gAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAg                                                >  1:3243045/1‑58 (MQ=255)
                                    gAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAg                                                >  1:4227343/1‑58 (MQ=255)
                                    gAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACgg                                            >  1:3763214/1‑62 (MQ=255)
                                       cTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGa                                <  1:2633676/71‑1 (MQ=255)
                                        tGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGgaga                                                            <  1:3556524/42‑1 (MQ=255)
                                        tGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTg                                 >  1:1519979/1‑69 (MQ=255)
                                          aCGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGaa                               <  1:4707404/69‑1 (MQ=255)
                                            gCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAACGACGGCTCCAt                                      <  1:3018839/60‑1 (MQ=255)
                                            gCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGaa                               <  1:2670975/67‑1 (MQ=255)
                                            gCCGAATATGGCAGGTCCCGCCTTCAACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAg                            <  1:1860395/70‑1 (MQ=255)
                                             ccGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCt                          <  1:3496803/71‑1 (MQ=255)
                                               gAATATGGCAGGTCCCGCCTTCGACGCGCTGGGGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGa                        >  1:4592261/1‑71 (MQ=255)
                                               gAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGa                        >  1:4433466/1‑71 (MQ=255)
                                               gAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGa                        >  1:1119572/1‑71 (MQ=255)
                                               gAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGa                        >  1:822826/1‑71 (MQ=255)
                                               gAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGa                        >  1:2823237/1‑71 (MQ=255)
                                                aaTATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACaaaa                                                   <  1:705372/43‑1 (MQ=255)
                                                  tatGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGt                     >  1:4757042/1‑71 (MQ=255)
                                                    tGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGtt                    >  1:678880/1‑70 (MQ=255)
                                                    tGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGtt                    >  1:2020452/1‑70 (MQ=255)
                                                    tGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGtt                    <  1:247359/70‑1 (MQ=255)
                                                     ggCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGtt                    <  1:2281236/69‑1 (MQ=255)
                                                      gCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGa                                               >  1:3469718/1‑41 (MQ=255)
                                                      gCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGa                                               >  1:3467987/1‑41 (MQ=255)
                                                       cAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCa                                          >  1:2130979/1‑45 (MQ=255)
                                                       cAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTaa                  >  1:3140150/1‑69 (MQ=255)
                                                       cAGGTCCCGCCTTCAACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGATAAc                 >  1:2456733/1‑70 (MQ=255)
                                                        aGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCa               <  1:417264/71‑1 (MQ=255)
                                                          gTCCCGCCTTCGACGCGCTTGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCaa              <  1:405908/70‑1 (MQ=255)
                                                          gTCCCGCCTTCGACGCGCTGGAGAAATACaaaaaa                                                 <  1:1010086/35‑1 (MQ=255)
                                                          gTCCCGCCTTCAACGCGCTGGAGAAATACaaaaaa                                                 <  1:4727935/35‑1 (MQ=255)
                                                          gTCCCGCCTTCAACGCGCTGGAGAAATACaaaaaa                                                 <  1:3661263/35‑1 (MQ=255)
                                                               gCCTTCGCCGCGCTGGAGAAATACAAAAAAGACGGCACAATGCCTGAAAAGCTGACGTTa                   >  1:697651/1‑60 (MQ=255)
                                                               gCCTTCAACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCAc        <  1:4646513/71‑1 (MQ=255)
                                                                  ttCGACGCGCTGGAGAAATACAAAAAAGACGGCACCa                                       >  1:3192119/1‑37 (MQ=255)
                                                                  ttCGACGCGCTGGAGAAATACAAAAAAGACGGCACCa                                       >  1:4044984/1‑37 (MQ=255)
                                                                  ttCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTa                   <  1:3106509/57‑1 (MQ=255)
                                                                  ttCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTa                   <  1:1594241/57‑1 (MQ=255)
                                                                  ttCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCt     >  1:1649899/1‑71 (MQ=255)
                                                                   tCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCAccc      <  1:4333949/69‑1 (MQ=255)
                                                                    cGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAAc                 >  1:1931326/1‑57 (MQ=255)
                                                                    cGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAAc                 >  1:1586103/1‑57 (MQ=255)
                                                                    cGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAAc                 >  1:4249090/1‑57 (MQ=255)
                                                                    cGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAAc                 >  1:2905249/1‑57 (MQ=255)
                                                                     gACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGaa                               <  1:2676535/42‑1 (MQ=255)
                                                                     gACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCttt   <  1:425307/70‑1 (MQ=255)
                                                                     gACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTa  <  1:708973/71‑1 (MQ=255)
                                                                     |                                                                      
AAAGCGATGATGGATGGCGAAGCGAACGCCAGTGTTGAACTGACGCCGAATATGGCAGGTCCCGCCTTCGACGCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTA  >  minE/2892108‑2892247

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: