Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,474,925 T→C 57.9% S19S (AGT→AGC pstS → phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,474,9250TC57.9% 39.8 / 74.0 102S19S (AGT→AGCpstSphosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (18/25);  new base C (35/24);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

ATGAAAGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGGTG  >  minE/2474869‑2474988
                                                        |                                                               
aTGAAAGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCtgtg                                                     <  1:563354/69‑1 (MQ=255)
   aaaGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGc                                                               <  1:1777083/56‑1 (MQ=255)
   aaaGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGc                                                               <  1:3318860/56‑1 (MQ=255)
   aaaGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGc                                                               <  1:3388738/56‑1 (MQ=255)
   aaaGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGc                                                               <  1:2028043/56‑1 (MQ=255)
   aaaGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGc                                                               <  1:671765/56‑1 (MQ=255)
   aaaGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGc                                                               <  1:853068/56‑1 (MQ=255)
   aaaGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGc                                                               <  1:121661/56‑1 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2218400/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1812299/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:3036771/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1846847/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2538008/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1990018/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:3406672/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2489026/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2100204/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2165025/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:325712/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2449104/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:3400088/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:229514/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2313226/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1051076/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2437171/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2443126/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2887033/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:3134756/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1225313/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:3068923/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1135433/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1357892/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:849189/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1473659/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2762372/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1643013/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:512118/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:3005337/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2605595/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:605861/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2593297/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:1088129/1‑68 (MQ=255)
         aTGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAGa                                             >  1:2589327/1‑68 (MQ=255)
                  accGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCt                                    <  1:15940/68‑1 (MQ=255)
                  accGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCt                                    <  1:2536799/68‑1 (MQ=255)
                     gTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAca                                >  1:105547/1‑69 (MQ=255)
                        gCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGaagc                                          >  1:2286350/1‑56 (MQ=255)
                        gCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcagg                              >  1:345084/1‑68 (MQ=255)
                        gCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcagg                              >  1:2447473/1‑68 (MQ=255)
                        gCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcagg                              >  1:2221331/1‑68 (MQ=255)
                        gCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcagg                              >  1:361202/1‑68 (MQ=255)
                        gCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcagg                              >  1:576803/1‑68 (MQ=255)
                        gCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcagg                              >  1:113288/1‑68 (MQ=255)
                        gCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGAAGaagc                                          >  1:1166069/1‑56 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcaggtgca                          >  1:768393/1‑69 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcaggtgca                          >  1:1566035/1‑69 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCTCTTTCTCTGTGTTTGCAg                                              <  1:1729902/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:3401137/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:1047897/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:489790/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:533113/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:885039/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2924818/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:886193/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:1396923/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:1765432/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:1722442/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2888161/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2002587/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2085520/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2293790/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2558873/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:1296792/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2709206/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:1005707/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2802071/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2841520/49‑1 (MQ=255)
                           aCTGTTGTCGCCGCGACCTTAACGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:2379420/49‑1 (MQ=39)
                            cTGTTGTCGCCGCGACCTTATCGATGAGCGCTTTCTCTGTGTTTGCAg                                              <  1:783870/48‑1 (MQ=255)
                              gttgtCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcaggtgcagg                        >  1:3289651/1‑68 (MQ=255)
                              gttgtCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcaggtgcagg                        >  1:1229432/1‑68 (MQ=255)
                              gttgtCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcaggtgcagg                        >  1:2065653/1‑68 (MQ=255)
                              gttgtCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGAcaggtgcagg                        >  1:2815083/1‑68 (MQ=255)
                                 gtCGCCGCGTCCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAggtgcaggtgca                    <  1:906421/69‑1 (MQ=255)
                                 gtCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAggtgcaggtgca                    <  1:1856159/69‑1 (MQ=255)
                                 gtCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAggtgcaggtgca                    <  1:1732408/69‑1 (MQ=255)
                                 gtCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAggtgcaggtgca                    <  1:431795/69‑1 (MQ=255)
                                 gtCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAggtgcaggtgca                    <  1:3301540/69‑1 (MQ=255)
                                                tCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAggtgcaggtgca                    <  1:2923062/54‑1 (MQ=255)
                                                 cGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCgg    >  1:174680/1‑69 (MQ=255)
                                                 cGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCgg    >  1:24470/1‑69 (MQ=255)
                                                 cGATGACTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCgg    >  1:436213/1‑69 (MQ=255)
                                                    tGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGgtg  <  1:1789866/68‑1 (MQ=255)
                                                    tGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGgtg  <  1:1789052/68‑1 (MQ=255)
                                                    tGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGgtg  <  1:1840784/68‑1 (MQ=255)
                                                    tGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGgtg  <  1:61959/68‑1 (MQ=255)
                                                    tGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGgtg  <  1:1701016/68‑1 (MQ=255)
                                                    tGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGgtg  <  1:3241232/68‑1 (MQ=255)
                                                    tGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGgtg  <  1:2778093/68‑1 (MQ=255)
                                                    tGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGgtg  <  1:3061964/68‑1 (MQ=255)
                                                        |                                                               
ATGAAAGTTATGCGTACCACCGTCGCAACTGTTGTCGCCGCGACCTTATCGATGAGTGCTTTCTCTGTGTTTGCAGAAGCAAGCCTGACAGGTGCAGGTGCAACCTTCCCTGCGCCGGTG  >  minE/2474869‑2474988

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: