Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,045,174 C→T 18.9% P175P (CCG→CCA pdxH ← pyridoxine 5'‑phosphate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,045,1740CT18.9% 103.0 / 14.0 68P175P (CCG→CCApdxHpyridoxine 5'‑phosphate oxidase
Reads supporting (aligned to +/- strand):  ref base C (8/47);  new base T (6/7);  total (14/54)
Fisher's exact test for biased strand distribution p-value = 2.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGAATTTACTTTCAAG  >  minE/1045109‑1045233
                                                                 |                                                           
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2806229/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3353595/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3350036/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3393244/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:1731120/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:482255/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1572601/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1561710/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1555102/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:593830/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2628096/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2640957/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:766520/69‑1 (MQ=255)
cGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:103585/69‑1 (MQ=255)
 gATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2327447/68‑1 (MQ=255)
  aTGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCaa                                                        <  1:796709/69‑1 (MQ=255)
     cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:157188/64‑1 (MQ=255)
     cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2958089/64‑1 (MQ=255)
     cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:1472981/64‑1 (MQ=255)
     cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:2464143/64‑1 (MQ=255)
     cTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          <  1:3412105/64‑1 (MQ=255)
         ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGc                                                    >  1:1579363/4‑66 (MQ=255)
         ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:1806786/4‑69 (MQ=255)
         ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:2110193/4‑69 (MQ=255)
         ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:1098924/4‑69 (MQ=255)
         ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:2676949/4‑69 (MQ=255)
         ccgcccTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAACTTGGCAATGGCACt                                                 >  1:1185249/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:1964277/47‑1 (MQ=255)
                      aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          >  1:322261/1‑47 (MQ=255)
                      aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGc                                                          >  1:2408605/1‑47 (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:2380671/69‑1 (MQ=255)
                      aCTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGaa                                    <  1:1671946/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:2658547/69‑1 (MQ=255)
                                            gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGaa              <  1:3082917/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)
                                            gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGa               <  1:1206628/68‑1 (MQ=255)
                                            gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGa               <  1:1849044/68‑1 (MQ=255)
                                            gcgAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGa               <  1:2504881/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)
                                                                 |                                                           
CGATGCTCACCACCCTGCCAGAACTCAATCTGTTCAAGGCTGACGCGAAAACCGCCCCAAAAGCTCGGCAATGGCACTTCGCCCTGTTGAAACTTCTGCTTCAGCTCCAGGAATTTACTTTCAAG  >  minE/1045109‑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: