Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,474,576 C→T 34.3% T720T (ACC→ACT nrdA → ribonucleoside diphosphate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,474,5760CT34.3% 70.0 / 48.8 90T720T (ACC→ACTnrdAribonucleoside diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (19/40);  new base T (16/15);  total (35/55)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-02

ACTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGC  >  minE/1474509‑1474638
                                                                   |                                                              
aCTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCg                                                               <  1:88424/69‑1 (MQ=255)
aCTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCg                                                               <  1:1390977/69‑1 (MQ=255)
aCTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCg                                                               <  1:1163326/69‑1 (MQ=255)
aCTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCg                                                               <  1:1791737/69‑1 (MQ=255)
aCTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCg                                                               <  1:1244570/69‑1 (MQ=255)
aCTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCg                                                               <  1:1445629/69‑1 (MQ=255)
 cTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCg                                                               <  1:1321509/68‑1 (MQ=255)
   aCGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCGcc                                                             >  1:2604203/1‑68 (MQ=255)
        ccGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCg                                                               <  1:3421585/61‑1 (MQ=255)
        ccGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACaaa                                                       <  1:355689/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:2521474/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:2581494/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:151977/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:1310183/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:250757/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:1725847/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:1753154/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:2614312/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:551899/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:2722382/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:1194577/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:2867767/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:3227849/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCg                                                   <  1:2196227/69‑1 (MQ=255)
            cACGTTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCAACAAATTCg                                                   <  1:3346393/69‑1 (MQ=255)
             aCGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCgg                                                  >  1:2114635/1‑69 (MQ=255)
                            aaaaGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGta                                   <  1:899070/69‑1 (MQ=255)
                             aaaGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGtat                                  >  1:901339/1‑69 (MQ=255)
                             aaaGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGtat                                  >  1:1814208/1‑69 (MQ=255)
                                   ccGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAg                            <  1:3045434/69‑1 (MQ=255)
                                   ccGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAg                            <  1:3315327/69‑1 (MQ=255)
                                   ccGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAg                            <  1:2217467/69‑1 (MQ=255)
                                   ccGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAg                            <  1:2051991/69‑1 (MQ=255)
                                   ccGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAg                            <  1:914971/69‑1 (MQ=255)
                                   ccGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTACCAg                            <  1:2764725/69‑1 (MQ=255)
                                   ccGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGGCAAAACACTGTATTATCAg                            <  1:2388969/69‑1 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGCCCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:2714516/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACaa                                                        >  1:1334953/1‑40 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:2941854/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:575489/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:1419704/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:2510798/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:1642182/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:2478048/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:2391938/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:1168617/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:2301066/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:3065145/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:1138146/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:97155/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:1008425/1‑52 (MQ=255)
                                    cGATGCAGCAGTTGCTGAAAGACCTGCTCACTGCCTACAAATTCGGGGTCaa                                            >  1:86035/1‑52 (MQ=255)
                                        gcagcaGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAacac                       >  1:2362489/1‑69 (MQ=255)
                                        gcagcaGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAacac                       >  1:2645844/1‑69 (MQ=255)
                                         cagcagTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACAcc                      <  1:1249096/69‑1 (MQ=255)
                                            cagTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGt                   <  1:1401865/69‑1 (MQ=255)
                                            cagTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGt                   <  1:1004866/69‑1 (MQ=255)
                                            cagTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGt                   <  1:511700/69‑1 (MQ=255)
                                            cagTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGt                   <  1:2083213/69‑1 (MQ=255)
                                            cagTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGt                   <  1:2426712/69‑1 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGt                                    <  1:2590316/47‑1 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGt                                    <  1:874723/47‑1 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGt                                    <  1:2296319/47‑1 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGt                                    <  1:2357268/47‑1 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGt                                    <  1:2469571/47‑1 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGt                                    <  1:1427862/47‑1 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGAc                >  1:1588760/1‑67 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACgg              >  1:1185558/1‑69 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACgg              >  1:2783432/1‑69 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACgg              >  1:3174526/1‑69 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACgg              >  1:3188399/1‑69 (MQ=255)
                                                 gCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACgg              >  1:648662/1‑69 (MQ=255)
                                                   tGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgc             <  1:342786/68‑1 (MQ=255)
                                                   tGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgc             <  1:1005816/68‑1 (MQ=255)
                                                   tGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgc             <  1:3037925/68‑1 (MQ=255)
                                                   tGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgc             <  1:2833686/68‑1 (MQ=255)
                                                   tGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgc             <  1:2674773/68‑1 (MQ=255)
                                                   tGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgc             <  1:1786178/68‑1 (MQ=255)
                                                   tGAAACACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgc             <  1:3344636/68‑1 (MQ=255)
                                                      aaGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgcgc           >  1:3242438/1‑67 (MQ=255)
                                                      aaGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgcgc           >  1:2684420/1‑67 (MQ=255)
                                                      aaGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgcgc           >  1:1961720/1‑67 (MQ=255)
                                                      aaGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgcgc           >  1:1779499/1‑67 (MQ=255)
                                                      aaGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgcgc           >  1:1292786/1‑67 (MQ=255)
                                                      aaGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGgcg            >  1:1128894/1‑66 (MQ=255)
                                                          ccTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGa     <  1:2624820/69‑1 (MQ=255)
                                                          ccTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGa     <  1:1908977/69‑1 (MQ=255)
                                                          ccTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGa     <  1:1896817/69‑1 (MQ=255)
                                                          ccTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGa     <  1:985943/69‑1 (MQ=255)
                                                              cTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTTAAGACGc  >  1:2190109/1‑68 (MQ=255)
                                                                   |                                                              
ACTACGATCCGTCACGCTTCCCGTCAGGAAAAGTGCCGATGCAGCAGTTGCTGAAAGACCTGCTCACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGC  >  minE/1474509‑1474638

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: