Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,045,171 G→A 18.1% S176S (AGC→AGT pdxH ← pyridoxine 5'‑phosphate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,045,1710GA18.1% 112.9 / 12.6 71S176S (AGC→AGTpdxHpyridoxine 5'‑phosphate oxidase
Reads supporting (aligned to +/- strand):  ref base G (11/47);  new base A (6/7);  total (17/54)
Fisher's exact test for biased strand distribution p-value = 6.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CAGGCGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGAATTTACTTTCAAG  >  minE/1045105‑1045233
                                                                  |                                                              
cAGGCGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCt                                                              >  1:2791466/1‑69 (MQ=255)
cAGGCGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCt                                                              >  1:1227840/1‑69 (MQ=255)
cAGGCGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCt                                                              >  1:705737/1‑69 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1572601/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3353595/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3393244/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3350036/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3099511/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1890395/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:254052/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:482255/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2628096/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1731120/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:593830/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1561710/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2806229/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2640957/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:103585/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:766520/69‑1 (MQ=255)
    cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1555102/69‑1 (MQ=255)
     gATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2327447/68‑1 (MQ=255)
      aTGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCaa                                                        <  1:796709/69‑1 (MQ=255)
         cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2958089/64‑1 (MQ=255)
         cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3412105/64‑1 (MQ=255)
         cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2464143/64‑1 (MQ=255)
         cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1472981/64‑1 (MQ=255)
         cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:157188/64‑1 (MQ=255)
             ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGc                                                    >  1:1579363/4‑66 (MQ=255)
             ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:1806786/4‑69 (MQ=255)
             ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:2676949/4‑69 (MQ=255)
             ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:1185249/4‑69 (MQ=255)
             ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:2110193/4‑69 (MQ=255)
             ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:1098924/4‑69 (MQ=255)
                         aaCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACTTCGCCCTGTTGa                                     <  1:34178/69‑1 (MQ=255)
                         aaCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACTTCGCCCTGTTGa                                     <  1:3380551/69‑1 (MQ=255)
                         aaCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACTTCGCCCTGTTGa                                     <  1:127887/69‑1 (MQ=255)
                         aaCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACTTCGCCCTGTTGa                                     <  1:1180043/69‑1 (MQ=255)
                         aaCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACTTCGCCCTGTTGa                                     <  1:109705/69‑1 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          >  1:2068530/1‑47 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          >  1:2408605/1‑47 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          >  1:322261/1‑47 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1964277/47‑1 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:192459/47‑1 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          >  1:1762664/1‑47 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:751152/47‑1 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          >  1:1351261/1‑47 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGaa                                    <  1:851158/69‑1 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGaa                                    <  1:3088854/69‑1 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGaa                                    <  1:1671946/69‑1 (MQ=255)
                          aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGaa                                    <  1:2380671/69‑1 (MQ=255)
                           cTCAATCTGTTCATGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2234914/46‑1 (MQ=255)
                             cAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACt                                 >  1:1369991/1‑69 (MQ=255)
                             cAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACt                                 >  1:2230780/1‑69 (MQ=255)
                             cAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACt                                 >  1:2058866/1‑69 (MQ=255)
                                             gACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACTTCGCCCTGTTGa                                     <  1:564057/49‑1 (MQ=255)
                                              aCGCGAAAACCGCCCCAAAAACTTGGCAATGGCACTTCGCCCTGTTGa                                     <  1:2454692/48‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGaa              <  1:2728306/69‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGaa              <  1:1327080/69‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGaa              <  1:1424527/69‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGaa              <  1:3434671/69‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGaa              <  1:2658547/69‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGaa              <  1:3082917/69‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGa               <  1:1849044/68‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGa               <  1:2504881/68‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGa               <  1:1206628/68‑1 (MQ=255)
                                                gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACCTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGa               <  1:2932699/68‑1 (MQ=255)
                                                 cgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGaa              <  1:1859674/68‑1 (MQ=255)
                                                            ccAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGAATTTACTTTCAAg  <  1:2262556/69‑1 (MQ=255)
                                                            ccAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGAATTTACTTTCAAg  <  1:3188460/69‑1 (MQ=255)
                                                            ccAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGAATTTACTTTCAAg  <  1:1269939/69‑1 (MQ=255)
                                                             cAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGAATTTACTTTCAAg  <  1:526554/68‑1 (MQ=255)
                                                                  |                                                              
CAGGCGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGAATTTACTTTCAAG  >  minE/1045105‑1045233

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: