Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,606,942 G→A 30.6% R220R (CGG→CGA hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,606,9420GA30.6% 86.7 / 42.3 91R220R (CGG→CGAhemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base A (17/11);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATG  >  minE/1606878‑1607003
                                                                |                                                             
ttGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAg                                                           <  1:2283270/69‑1 (MQ=255)
ttGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAg                                                           <  1:66009/69‑1 (MQ=255)
  gCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGc                                                          <  1:351703/68‑1 (MQ=255)
  gCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGc                                                          <  1:543176/68‑1 (MQ=255)
  gCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGc                                                          <  1:957913/68‑1 (MQ=255)
  gCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGc                                                          <  1:633288/68‑1 (MQ=255)
  gCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGc                                                          <  1:81126/68‑1 (MQ=255)
  gCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGc                                                          <  1:1994661/68‑1 (MQ=255)
  gCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGc                                                          <  1:757944/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:1030197/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:290982/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:3222704/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:2820939/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:2816897/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:608437/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:707578/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:1858331/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:1284505/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGACGAAAAGCg                                                         <  1:890757/68‑1 (MQ=255)
   ccTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTAGAGCGAAGAAAAGCg                                                         <  1:3341400/68‑1 (MQ=255)
         tGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGcc                                                   <  1:3104179/68‑1 (MQ=255)
         tGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGcc                                                   <  1:3377101/68‑1 (MQ=255)
         tGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGcc                                                   <  1:2707726/68‑1 (MQ=255)
         tGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGcc                                                   <  1:135271/68‑1 (MQ=255)
         tGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGcc                                                   <  1:3170365/68‑1 (MQ=255)
         tGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGcc                                                   <  1:1691117/68‑1 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:2246983/1‑68 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:2173394/1‑68 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:1512528/1‑68 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:1615759/1‑68 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:2772788/1‑68 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:2439202/1‑68 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:384107/1‑68 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:2766688/1‑68 (MQ=255)
             ggCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACg                                               >  1:2938385/1‑68 (MQ=255)
                   aGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGcc                                                   <  1:2104927/58‑1 (MQ=255)
                   aGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGgcgagc                                         >  1:3399082/1‑68 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGAt                                                       >  1:731846/1‑47 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:3314931/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:2048872/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:2898975/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:2835993/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:3392070/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:1304658/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:2816007/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:2471499/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:2550784/1‑69 (MQ=255)
                          ggCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCgagcgcgagcgc                                 >  1:2575324/1‑69 (MQ=255)
                           gCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGc                                                          <  1:2763496/43‑1 (MQ=255)
                                 ccGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTcc                          <  1:928440/69‑1 (MQ=255)
                                 ccGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTcc                          <  1:1157445/69‑1 (MQ=255)
                                         tATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTatat                  >  1:1668813/1‑69 (MQ=255)
                                             tACCAATTGTCGAGCGACGGAAAGGGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTtat                    >  1:138293/1‑63 (MQ=255)
                                             tACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTtat                    >  1:1076341/1‑63 (MQ=255)
                                             tACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTtat                    >  1:1980072/1‑63 (MQ=255)
                                             tACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTtat                    >  1:2452398/1‑63 (MQ=255)
                                             tACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTtat                    >  1:2638562/1‑63 (MQ=255)
                                                cAATTGTAGAGCGATGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgcg           >  1:2667783/1‑69 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCTAGCGCAATTTTCAGCTCTACCGTcgct           >  1:489594/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTAccgtcg             >  1:2841887/1‑67 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTAccgtcg             >  1:883392/1‑67 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTAccgtcg             >  1:235110/1‑67 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTAccgtcg             >  1:2476730/1‑67 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgcg           >  1:321740/1‑69 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgcg           >  1:149185/1‑69 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgc            >  1:1513662/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgc            >  1:917050/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgc            >  1:654318/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgc            >  1:2163848/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgc            >  1:3342627/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgc            >  1:940783/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACCGTcgc            >  1:1929518/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTACAGTcgc            >  1:1403978/1‑68 (MQ=255)
                                                cAATTGTAGAGCGACGAAAAGCGATGCCCTACGGCGAGCGCGAGCGCAATTTTCAGCTCTAccgtcg             >  1:862146/1‑67 (MQ=255)
                                                    tGTCGAGCGACGGAAAGCGATGGGCTACGGCGAGCGAGAGCGCAATTTCCAGTta                     <  1:3344267/55‑1 (MQ=255)
                                                    tGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTta                     >  1:447697/1‑55 (MQ=255)
                                                    tGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTta                     <  1:2616822/55‑1 (MQ=255)
                                                    tGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTta                     <  1:337501/55‑1 (MQ=255)
                                                    tGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTta                     <  1:2818768/55‑1 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      <  1:2375599/69‑1 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      >  1:833656/1‑69 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      >  1:1596271/1‑69 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      <  1:803262/69‑1 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      >  1:1896534/1‑69 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      <  1:2279279/69‑1 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      <  1:558027/69‑1 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      <  1:3132332/69‑1 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      <  1:398538/69‑1 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      <  1:358077/69‑1 (MQ=255)
                                                     gTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGt      <  1:2739858/69‑1 (MQ=255)
                                                          gCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATg  <  1:3012475/68‑1 (MQ=255)
                                                                |                                                             
TTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATG  >  minE/1606878‑1607003

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: