Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,290,614 G→A 54.6% P442P (CCG→CCA ppc → phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,290,6140GA54.6% 5.0 / 43.6 55P442P (CCG→CCAppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base G (6/19);  new base A (22/8);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 3.90e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTGATTGCCGAAGC  >  minE/2290564‑2290679
                                                  |                                                                 
aaaCAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAgc                                                 >  1:202168/1‑69 (MQ=255)
aaaCAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAgc                                                 >  1:1415531/1‑69 (MQ=255)
  aCAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAgcg                                                >  1:2392992/1‑68 (MQ=255)
  aCAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAgcg                                                >  1:1543282/1‑68 (MQ=255)
   cAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTACCACGCAACTGGCAACGAAGCGcc                                              <  1:3103149/69‑1 (MQ=255)
   cAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGcc                                              <  1:1218655/69‑1 (MQ=255)
   cAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGcc                                              <  1:1890102/69‑1 (MQ=255)
   cAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGcc                                              <  1:1706367/69‑1 (MQ=255)
   cAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGcc                                              <  1:3308789/69‑1 (MQ=255)
   cAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGcc                                              <  1:526992/69‑1 (MQ=255)
   cAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGcc                                              <  1:737453/69‑1 (MQ=255)
   cAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGcc                                              <  1:997602/69‑1 (MQ=255)
           ccTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAAcgc                                       >  1:152723/1‑68 (MQ=255)
           ccTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAAcgc                                       >  1:579750/1‑68 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:3298691/69‑1 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:598417/69‑1 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:2984537/69‑1 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:2962175/69‑1 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:1232411/69‑1 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:2875254/69‑1 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:2711120/69‑1 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:2413532/69‑1 (MQ=255)
                     gaactgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:2603784/69‑1 (MQ=255)
                        ctgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGAt                         <  1:215277/69‑1 (MQ=255)
                        ctgaactCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGAt                         <  1:3328730/69‑1 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:1493032/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:1358087/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:947325/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:883312/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:864300/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:862442/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:754697/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:1369211/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:674116/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:1396331/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:396917/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:3387916/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:2415583/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:2387202/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:3162365/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:3171644/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:3229781/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:151401/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:3338289/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:344828/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACCACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTAGATACCt                     >  1:593246/1‑68 (MQ=255)
                             ctCCAAACGTCCGCTTCTACACGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCt                     >  1:86565/1‑67 (MQ=255)
                                    cGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTg             <  1:556148/69‑1 (MQ=255)
                                    cGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTg             <  1:946624/69‑1 (MQ=255)
                                    cGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTg             <  1:2830658/69‑1 (MQ=255)
                                    cGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTg             <  1:976137/69‑1 (MQ=255)
                                     gTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTg             <  1:778832/68‑1 (MQ=255)
                                     gTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTg             <  1:1763916/68‑1 (MQ=255)
                                       ccGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTc                            <  1:1700074/51‑1 (MQ=255)
                                                ccGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTGATTGCCGAAGc  <  1:1392793/68‑1 (MQ=255)
                                                  gcgcAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTGATTGCCGAAGc  <  1:839750/66‑1 (MQ=255)
                                                  |                                                                 
AAACAGGCGTTCCTGATCCGCGAACTGAACTCCAAACGTCCGCTTCTGCCGCGCAACTGGCAACCAAGCGCCGAAACGCGCGAAGTGCTCGATACCTGCCAGGTGATTGCCGAAGC  >  minE/2290564‑2290679

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: