Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 172,545 C→T 64.3% T67M (ACG→ATG)  dxr → 1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE172,5450CT64.3% 32.5 / 31.4 56T67M (ACG→ATG) dxr1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase
Reads supporting (aligned to +/- strand):  ref base C (8/12);  new base T (3/33);  total (11/45)
Fisher's exact test for biased strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGGAATTCTCTCCCCGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGA  >  minE/172485‑172605
                                                            |                                                            
tGGAATTCTCTCCCCGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTAc                                                      <  1:2907435/69‑1 (MQ=255)
            cccTCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:1137892/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:652980/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:3192986/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:2288931/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:1999038/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:185050/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:1814915/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:1347817/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:2978148/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:1765699/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:1686443/67‑1 (MQ=255)
            cccGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:160392/67‑1 (MQ=255)
                  aTGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:3213463/61‑1 (MQ=255)
                      cGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:2710840/57‑1 (MQ=255)
                        tAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAggg                                            <  1:1393303/55‑1 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCtt                          >  1:2276541/1‑68 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         >  1:1949207/1‑69 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         <  1:739433/69‑1 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         <  1:2179582/69‑1 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         <  1:1222656/69‑1 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         <  1:1322564/69‑1 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         <  1:178375/69‑1 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         <  1:1727479/69‑1 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         >  1:1413038/1‑69 (MQ=255)
                             gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                         <  1:3348370/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:3038119/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:2903241/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:379890/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:3111951/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:1994705/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:2437067/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:2391087/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:1521330/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:1434367/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:1414435/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:1372262/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:124205/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:1174306/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:1161065/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:1019702/69‑1 (MQ=255)
                                  tGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTg                     >  1:3255877/1‑68 (MQ=255)
                                  tGAAGCGAGTGCGAAACTCCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:2993872/69‑1 (MQ=255)
                                   gAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTTAGTgg                    <  1:3400469/68‑1 (MQ=255)
                                   gAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    >  1:1070555/1‑68 (MQ=255)
                                   gAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    >  1:3373417/1‑68 (MQ=255)
                                   gAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    >  1:3027469/1‑68 (MQ=255)
                                   gAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                    <  1:2934539/68‑1 (MQ=255)
                                                 aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCtt      >  1:2761446/1‑68 (MQ=255)
                                                 aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCtt      >  1:2739467/1‑68 (MQ=255)
                                                 aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCtt      >  1:249488/1‑68 (MQ=255)
                                                 aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCtt      >  1:78098/1‑68 (MQ=255)
                                                   ttcttAAAACGATGCTACAGCAACAGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGa  <  1:460641/69‑1 (MQ=255)
                                                    tcttAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGa  <  1:2884133/69‑1 (MQ=255)
                                                    tcttAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGa  <  1:361643/69‑1 (MQ=255)
                                                    tcttAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGa  <  1:2711709/69‑1 (MQ=255)
                                                            |                                                            
TGGAATTCTCTCCCCGCTATGCCGTAATGGACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGA  >  minE/172485‑172605

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: