Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,475,665 A→T 34.4% 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,6650AT34.4% 59.9 / 51.9 58T244S (ACC→TCC) nrdBribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like
Reads supporting (aligned to +/- strand):  ref base A (29/9);  new base T (9/11);  total (38/20)
Fisher's exact test for biased strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

AGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAT  >  minE/1475595‑1475729
                                                                      |                                                                
aGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCAGCGACGAAGCCCTGCACCTGa                                                                  <  1:1696308/71‑1 (MQ=255)
        aaTTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCAcc                                                          >  1:4050549/1‑71 (MQ=255)
              tGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCa                                                     <  1:1933643/70‑1 (MQ=255)
                gAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCatat                                                  >  1:2415058/1‑71 (MQ=255)
                gAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCatat                                                  >  1:4019687/1‑71 (MQ=255)
                gAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCatat                                                  >  1:1856817/1‑71 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCTGCACCCAGCATATGCTGa                                             <  1:3435938/71‑1 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCata                                                   >  1:27364/1‑65 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCata                                                   >  1:2409515/1‑65 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGa                                             <  1:2581722/71‑1 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGa                                             <  1:2782624/71‑1 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGa                                             <  1:2153528/71‑1 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGc                                                >  1:2760631/1‑68 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGc                                                >  1:2381027/1‑68 (MQ=255)
                     cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGc                                                >  1:3496697/1‑68 (MQ=255)
                      aaCGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGaa                                            >  1:479337/1‑71 (MQ=255)
                      aaCGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGaa                                            >  1:544233/1‑71 (MQ=255)
                        cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAAt                                           >  1:3446092/1‑70 (MQ=255)
                        cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATc                                          >  1:2997529/1‑71 (MQ=255)
                        cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATc                                          >  1:2786424/1‑71 (MQ=255)
                        cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATc                                          >  1:3876872/1‑71 (MQ=255)
                         gCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATct                                         >  1:3116808/1‑71 (MQ=255)
                          ccAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAAt                                           >  1:3847034/1‑68 (MQ=255)
                          ccAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAAt                                           >  1:854641/1‑68 (MQ=255)
                          ccAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAAt                                           >  1:903938/1‑68 (MQ=255)
                               attattCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTgcg                                   >  1:3695916/1‑71 (MQ=255)
                                ttattCGCCTGATTGCCAGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATctgc                                       >  1:2123061/1‑66 (MQ=255)
                                    tCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAgcg                              >  1:136550/1‑71 (MQ=255)
                                       ccTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAgcggcg                           <  1:2991391/71‑1 (MQ=255)
                                         tGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:3193270/1‑71 (MQ=255)
                                         tGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTg                                              <  1:697335/50‑1 (MQ=255)
                                         tGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTg                                              <  1:3276048/50‑1 (MQ=255)
                                          gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:3034946/1‑70 (MQ=255)
                                          gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         >  1:2236627/1‑70 (MQ=255)
                                           aTTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         <  1:1052986/69‑1 (MQ=255)
                                           aTTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         <  1:1296140/69‑1 (MQ=255)
                                           aTTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg                         <  1:328083/69‑1 (MQ=255)
                                              gCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCCGCTGCGCAGCGGCGCGGACg                     <  1:2263600/70‑1 (MQ=255)
                                               cccGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGa                    <  1:3145287/70‑1 (MQ=255)
                                                    gacgaAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              <  1:785782/71‑1 (MQ=255)
                                                    gacgaAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              <  1:781035/71‑1 (MQ=255)
                                                    gacgaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga              <  1:1982986/71‑1 (MQ=255)
                                                       gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt           >  1:3000900/1‑71 (MQ=255)
                                                       gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt           >  1:3113611/1‑71 (MQ=255)
                                                       gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt           >  1:3783330/1‑71 (MQ=255)
                                                       gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt           >  1:3561668/1‑71 (MQ=255)
                                                         aGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATctgc                                       >  1:1085509/1‑41 (MQ=255)
                                                         aGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATga         >  1:3807341/1‑70 (MQ=255)
                                                             cTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa     <  1:3885098/71‑1 (MQ=255)
                                                             cTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa     <  1:1212125/71‑1 (MQ=255)
                                                             cTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa     <  1:1929052/71‑1 (MQ=255)
                                                             cTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa     <  1:3312226/71‑1 (MQ=255)
                                                               gCACCTGACCGGCACCCAGCATATGCTGAATCTTCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaaa   >  1:2761962/1‑71 (MQ=255)
                                                               gCACATGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGAGCCTGAGATGGCGGaaa   >  1:1708579/1‑71 (MQ=255)
                                                                cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaaa   >  1:670780/1‑70 (MQ=255)
                                                                cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt  >  1:3538549/1‑71 (MQ=255)
                                                                cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt  >  1:3282385/1‑71 (MQ=255)
                                                                cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt  >  1:897421/1‑71 (MQ=255)
                                                                cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt  >  1:3325665/1‑71 (MQ=255)
                                                                      |                                                                
AGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAT  >  minE/1475595‑1475729

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: