Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,475,665 A→T 48.6% T244S (ACC→TCC)  nrdB → ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,475,6650AT48.6% ‑5.9 / 103.4 70T244S (ACC→TCC) nrdBribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like
Reads supporting (aligned to +/- strand):  ref base A (27/9);  new base T (18/16);  total (45/25)
Fisher's exact test for biased strand distribution p-value = 8.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

GCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAT  >  minE/1475600‑1475729
                                                                 |                                                                
gcgAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGc                                                             <  1:2520440/71‑1 (MQ=255)
gcgAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGc                                                             <  1:2700037/71‑1 (MQ=255)
       gATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAg                                                       <  1:277064/70‑1 (MQ=255)
        aTGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCa                                                     <  1:415052/71‑1 (MQ=255)
        aTGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCa                                                     <  1:4196738/71‑1 (MQ=255)
        aTGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCa                                                     <  1:3090789/71‑1 (MQ=255)
               gCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCt                                               >  1:2192071/1‑70 (MQ=255)
                   cGCCGAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTTTCCGGCACCCAGc                                                      <  1:3071682/59‑1 (MQ=255)
                   cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGc                                                      <  1:1201954/59‑1 (MQ=255)
                    gCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCCGCata                                                   >  1:2737361/1‑61 (MQ=255)
                    gCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCata                                                   >  1:3540688/1‑61 (MQ=255)
                    gCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCata                                                   >  1:2479399/1‑61 (MQ=255)
                    gCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCata                                                   >  1:369456/1‑61 (MQ=255)
                    gCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATct                                         >  1:2227613/1‑71 (MQ=255)
                     ccAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATctg                                        >  1:351621/1‑71 (MQ=255)
                     ccAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATctg                                        >  1:364272/1‑71 (MQ=255)
                       aaaaTTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATctgct                                      <  1:956886/71‑1 (MQ=255)
                           ttattCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTgcgc                                  <  1:2074900/71‑1 (MQ=255)
                              ttCGCCTGATTGCCTGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAgc                               <  1:1520514/71‑1 (MQ=255)
                              ttCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAg                                <  1:2620987/70‑1 (MQ=255)
                                cGCCTGATTGCCCGCGCCGAAGCCCTGCACCTGTCCGGCCCCCAGCATATGCTGAATCTGCTGCGCAg                                >  1:1852713/1‑68 (MQ=255)
                                cGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAg                                >  1:1526752/1‑68 (MQ=255)
                                cGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAg                                >  1:70061/1‑68 (MQ=255)
                                cGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAg                                >  1:1289253/1‑68 (MQ=255)
                                  ccTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGc                                                             <  1:3936258/37‑1 (MQ=255)
                                    tGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTgcgc                                  <  1:2369457/62‑1 (MQ=255)
                                    tGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:3310229/1‑71 (MQ=255)
                                    tGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:1212557/1‑71 (MQ=255)
                                     gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAg                                <  1:3012560/63‑1 (MQ=255)
                                     gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:2250382/1‑70 (MQ=255)
                                     gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:455992/1‑70 (MQ=255)
                                     gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:3875364/1‑70 (MQ=255)
                                     gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:2743611/1‑70 (MQ=255)
                                          cccGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGa                    <  1:200174/70‑1 (MQ=255)
                                          cccGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGa                    <  1:144582/70‑1 (MQ=255)
                                          cccGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGAt                   <  1:2978092/71‑1 (MQ=255)
                                          cccGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGa                    <  1:3876784/70‑1 (MQ=255)
                                            cgcgACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAg                                <  1:222416/56‑1 (MQ=255)
                                             gcgACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCt                >  1:4066903/1‑71 (MQ=255)
                                               gacgaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCt                >  1:324433/1‑69 (MQ=255)
                                               gacgaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCt                >  1:922770/1‑69 (MQ=255)
                                               gacgaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCt                >  1:541663/1‑69 (MQ=255)
                                               gacgaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              <  1:1198786/71‑1 (MQ=255)
                                               gacgaAGCCCTGCACCTGACCGGCACCCAGCATAAGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTg               >  1:2690541/1‑70 (MQ=255)
                                                 cgaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACg                     >  1:3437178/1‑62 (MQ=255)
                                                 cgaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCt                >  1:2324019/1‑67 (MQ=255)
                                                  gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              >  1:923235/1‑68 (MQ=255)
                                                  gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              >  1:1727885/1‑68 (MQ=255)
                                                  gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              >  1:2170450/1‑68 (MQ=255)
                                                    aGCCCTGCACCTGTCCGGCACCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGc        >  1:492917/1‑71 (MQ=255)
                                                      cccTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              >  1:2071647/1‑64 (MQ=255)
                                                      cccTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              >  1:2413883/1‑64 (MQ=255)
                                                      cccTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              >  1:3816489/1‑64 (MQ=255)
                                                      cccTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTg               >  1:2017773/1‑63 (MQ=255)
                                                       ccTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCg       >  1:975543/1‑70 (MQ=255)
                                                        cTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa     <  1:1839362/71‑1 (MQ=255)
                                                        cTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa     <  1:290611/71‑1 (MQ=255)
                                                        cTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa     <  1:295652/71‑1 (MQ=255)
                                                        cTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa     <  1:3033548/71‑1 (MQ=255)
                                                          gCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaaa   >  1:1932230/1‑71 (MQ=255)
                                                          gCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaaa   >  1:119628/1‑71 (MQ=255)
                                                          gCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaaa   >  1:30341/1‑71 (MQ=255)
                                                           cACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt  >  1:2701638/1‑71 (MQ=255)
                                                           cACCTGTCCGGCAACCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt  >  1:2737911/1‑71 (MQ=255)
                                                           cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGa                       >  1:3289723/1‑50 (MQ=255)
                                                           cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGa                       >  1:1225759/1‑50 (MQ=255)
                                                           cACCTGACCGGCACCCAGCATATGCAGAATCTGCTGCGCAGCGGCGCgg                        >  1:1281017/1‑49 (MQ=255)
                                                               tGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTgaga            >  1:2846701/1‑57 (MQ=255)
                                                                gACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTg               >  1:1899434/1‑53 (MQ=255)
                                                                gACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTg               >  1:423402/1‑53 (MQ=255)
                                                                 |                                                                
GCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAT  >  minE/1475600‑1475729

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: