Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,475,665 A→T 20.7% 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,6650AT20.7% 46.7 / 7.3 24T244S (ACC→TCC) nrdBribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base T (2/3);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATT  >  minE/1475595‑1475730
                                                                      |                                                                 
aGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGa                                                                   >  1:227064/1‑71 (MQ=255)
            gATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAg                                                        >  1:1376360/1‑70 (MQ=255)
             aTGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCa                                                      <  1:1460842/71‑1 (MQ=255)
              tGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCat                                                     <  1:46787/71‑1 (MQ=255)
                gAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCata                                                    <  1:37548/70‑1 (MQ=255)
                 aaGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATg                                                  <  1:253029/71‑1 (MQ=255)
                   ggCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGc                                                 <  1:1811072/70‑1 (MQ=255)
                       aCGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCt                                                >  1:1743559/1‑67 (MQ=255)
                               attattCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCtgctgct                                    >  1:1987174/1‑70 (MQ=255)
                                ttattCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCAc                                                            >  1:559007/1‑46 (MQ=255)
                                ttattCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTgcg                                    >  1:198807/1‑70 (MQ=255)
                                 tattCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTgcg                                    >  1:1199465/1‑69 (MQ=255)
                                  attCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCAccc                                                          >  1:416072/1‑46 (MQ=255)
                                   ttCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAg                                 <  1:570890/70‑1 (MQ=255)
                                   ttCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAg                                 <  1:896178/70‑1 (MQ=255)
                                            ttGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATctg                                         <  1:1732539/53‑1 (MQ=255)
                                            ttGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATctg                                         <  1:906883/53‑1 (MQ=255)
                                               cccGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACg                      <  1:114863/69‑1 (MQ=255)
                                                      cgaAGCCCTGCACCTGTCCGGCACCCAGCATATGc                                                 <  1:936692/35‑1 (MQ=255)
                                                       gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt            >  1:1534977/1‑71 (MQ=255)
                                                       gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt            >  1:742168/1‑71 (MQ=255)
                                                               gCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaaa    >  1:52975/1‑71 (MQ=255)
                                                                 aCCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAtt  >  1:1222518/1‑71 (MQ=255)
                                                                  ccTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaa     <  1:358460/67‑1 (MQ=255)
                                                                     gACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTgaga             >  1:2162547/1‑56 (MQ=255)
                                                                      |                                                                 
AGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAATT  >  minE/1475595‑1475730

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: