Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 172,564 T→G 42.9% G73G (GGT→GGG 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,5640TG42.9% 8.2 / 49.2 49G73G (GGT→GGGdxr1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase
Reads supporting (aligned to +/- strand):  ref base T (13/15);  new base G (3/18);  total (16/33)
Fisher's exact test for biased strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCA  >  minE/172514‑172632
                                                  |                                                                    
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCtt                                                     >  1:2276541/1‑68 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    <  1:1322564/69‑1 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    >  1:1949207/1‑69 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    <  1:178375/69‑1 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    <  1:1727479/69‑1 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    <  1:2179582/69‑1 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    >  1:1413038/1‑69 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    <  1:1222656/69‑1 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    <  1:3348370/69‑1 (MQ=255)
gACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTa                                                    <  1:739433/69‑1 (MQ=255)
     tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               <  1:2391087/69‑1 (MQ=255)
     tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               <  1:2903241/69‑1 (MQ=255)
     tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               <  1:2437067/69‑1 (MQ=255)
     tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               <  1:379890/69‑1 (MQ=255)
     tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               <  1:1994705/69‑1 (MQ=255)
     tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               <  1:3111951/69‑1 (MQ=255)
     tGAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               <  1:3038119/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:3373417/1‑68 (MQ=255)
      gAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               >  1:1070555/1‑68 (MQ=255)
      gAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               >  1:3027469/1‑68 (MQ=255)
      gAAGCGAGTGCGAAACTTCTTAAAATGATGCTACAGCAACAGGGGAGCCGCACCGAAGTCTTAAGTgg                                               <  1:2934539/68‑1 (MQ=255)
                    aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCtt                                 >  1:78098/1‑68 (MQ=255)
                    aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCtt                                 >  1:2761446/1‑68 (MQ=255)
                    aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCtt                                 >  1:2739467/1‑68 (MQ=255)
                    aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCtt                                 >  1:249488/1‑68 (MQ=255)
                      ttcttAAAACGATGCTACAGCAACAGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGa                             <  1:460641/69‑1 (MQ=255)
                       tcttAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGa                             <  1:2884133/69‑1 (MQ=255)
                       tcttAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGa                             <  1:2711709/69‑1 (MQ=255)
                       tcttAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGa                             <  1:361643/69‑1 (MQ=255)
                                   gCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAgcgc                  >  1:2074235/1‑68 (MQ=255)
                                      aCAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTg               >  1:2521897/1‑68 (MQ=255)
                                      aCAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTg               >  1:2441369/1‑68 (MQ=255)
                                      aCAGCAACAGGGTAGCCGCACCGAAGTCTAAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTg               >  1:3267360/1‑68 (MQ=255)
                                         gCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAg             <  1:2351818/67‑1 (MQ=255)
                                         gCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAg             <  1:631952/67‑1 (MQ=255)
                                              aGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGTTTGAGGATGTTg      <  1:2651915/69‑1 (MQ=255)
                                              aGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGttt      <  1:1006114/69‑2 (MQ=255)
                                                 gTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATc   >  1:1038806/1‑69 (MQ=255)
                                                  tAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCa  >  1:2970285/1‑69 (MQ=255)
                                                  |                                                                    
GACGATGAAGCGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCA  >  minE/172514‑172632

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: