Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,450,176 A→G 66.2% L253L (CTA→CTG atoC → fused response regulator of ato opeon, in two‑component system with AtoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,450,1760AG66.2% 41.1 / 28.2 56L253L (CTA→CTGatoCfused response regulator of ato opeon, in two‑component system with AtoS
Reads supporting (aligned to +/- strand):  ref base A (8/11);  new base G (26/11);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 4.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATTATTTGAACGAGCCAACGAAGGTACTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGACCA  >  minE/1450110‑1450240
                                                                  |                                                                
attattTGAACGAGCCAACGAAGGTACTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGtactac                                                                 >  1:649542/1‑68 (MQ=255)
attattTGAACGAGCCAACGAAGGTACTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACa                                                                >  1:1941169/1‑69 (MQ=255)
attattTGAACGAGCCAACGAAGGTACTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACa                                                                >  1:1155503/1‑69 (MQ=255)
               cAACGAAGGTACTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACg                                                 <  1:2887980/69‑1 (MQ=255)
                    aaGGTACTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAAt                                                        >  1:1525087/1‑57 (MQ=255)
                         aCTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACa                                        >  1:1941243/1‑68 (MQ=255)
                         aCTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACa                                        >  1:1337969/1‑68 (MQ=255)
                         aCTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACa                                        >  1:657925/1‑68 (MQ=255)
                         aCTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACa                                        >  1:1734185/1‑68 (MQ=255)
                         aCTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACa                                        >  1:3435191/1‑68 (MQ=255)
                                    cgacgaAATTGGCGAAATGCCGCTGTTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:2055339/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:3358287/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:575331/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:497006/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:2322474/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:3094802/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:3179195/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:2737874/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:2790609/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:2873910/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:2908188/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:2978434/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:1354389/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:173102/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:97698/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:1420513/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:1516857/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:1808695/1‑68 (MQ=39)
                                    cgacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAt                             >  1:811498/1‑68 (MQ=39)
                                     gacgaAATTGGGGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAtt                            >  1:2306005/1‑68 (MQ=39)
                                     gacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAtt                            >  1:80111/1‑68 (MQ=39)
                                     gacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAtt                            >  1:1288381/1‑68 (MQ=39)
                                     gacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAtt                            >  1:1398607/1‑68 (MQ=39)
                                     gacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAtt                            >  1:3201237/1‑68 (MQ=39)
                                     gacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAtt                            >  1:1722745/1‑68 (MQ=39)
                                     gacgaAATTGGCGAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAAtt                            >  1:713037/1‑68 (MQ=39)
                                       cgaAATTGGCGAAATGCCGCTGGTACTATAAACCAAATTACTACGCATTCTACAGGAACGGGAATTTg                          <  1:1770673/68‑1 (MQ=255)
                                       cgaAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTg                          <  1:2639238/68‑1 (MQ=255)
                                          aaTTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACg                      <  1:2536322/69‑1 (MQ=255)
                                           aTTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACgg                     <  1:736194/69‑1 (MQ=255)
                                           aTTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACgg                     <  1:258645/69‑1 (MQ=255)
                                           aTTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACgg                     <  1:2522916/69‑1 (MQ=255)
                                           aTTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACgg                     <  1:2101434/69‑1 (MQ=255)
                                                 gAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:2498925/68‑1 (MQ=255)
                                                 gAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:1716491/68‑1 (MQ=255)
                                                 gAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:245377/68‑1 (MQ=255)
                                                 gAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:968904/68‑1 (MQ=255)
                                                 gAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:1770064/68‑1 (MQ=255)
                                                 gAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:880609/68‑1 (MQ=255)
                                                 gAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:1820421/68‑1 (MQ=255)
                                                 gAAATGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:1576882/68‑1 (MQ=255)
                                                  aaaTGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:20985/67‑1 (MQ=255)
                                                  aaaTGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:433910/67‑1 (MQ=255)
                                                   aaTGCCGCTGGTACTGCAAGCCAAATTACTACGCATTCTGCAGGAACGGGAATTTGAACGGATTgg                <  1:808940/66‑1 (MQ=255)
                                                       ccGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCAt         <  1:1901966/69‑1 (MQ=255)
                                                       ccGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCAt         <  1:744592/69‑1 (MQ=255)
                                                       ccGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCAt         <  1:1659157/69‑1 (MQ=255)
                                                               actacAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGAc    >  1:2081799/1‑66 (MQ=255)
                                                               actacAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGACCa  >  1:1216601/1‑68 (MQ=255)
                                                               actacAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGACCa  >  1:3018968/1‑68 (MQ=255)
                                                               actacAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGACCa  >  1:641836/1‑68 (MQ=255)
                                                               actacAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGACCa  >  1:490052/1‑68 (MQ=255)
                                                               actacAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGACCa  >  1:347517/1‑68 (MQ=255)
                                                               actacAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGACCa  >  1:3325253/1‑68 (MQ=255)
                                                                  |                                                                
ATTATTTGAACGAGCCAACGAAGGTACTCTGCTCCTCGACGAAATTGGCGAAATGCCGCTGGTACTACAAGCCAAATTACTACGCATTCTACAGGAACGGGAATTTGAACGGATTGGCGGCCATCAGACCA  >  minE/1450110‑1450240

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: