Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,657,788 C→T 34.4% P687S (CCT→TCT)  ppk → polyphosphate kinase, component of RNA degradosome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,657,7880CT34.4% 39.5 / 27.9 55P687S (CCT→TCT) ppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base C (18/18);  new base T (8/11);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

CGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGA  >  minE/1657731‑1657850
                                                         |                                                              
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:962748/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:2564504/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:1004/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:1929319/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:3182708/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:3274065/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:1653724/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:164154/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:1591926/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:128395/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:1214108/69‑1 (MQ=255)
cgccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAAccc                                                     <  1:738982/69‑1 (MQ=255)
 gccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCt                                                    >  1:2659211/1‑69 (MQ=255)
 gccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCt                                                    >  1:2524963/1‑69 (MQ=255)
 gccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCt                                                    >  1:1515356/1‑69 (MQ=255)
 gccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCt                                                    >  1:243363/1‑69 (MQ=255)
 gccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCt                                                    >  1:2357397/1‑69 (MQ=255)
 gccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCt                                                    >  1:563390/1‑69 (MQ=255)
 gccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCt                                                    >  1:2069128/1‑69 (MQ=255)
 gccgcAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCt                                                    >  1:2017424/1‑69 (MQ=255)
        aGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCaa                                             <  1:3332082/69‑1 (MQ=255)
                     ttGGCGATTTACGACTACATCAAATCACTCGAACAATCTGAATAAc                                                       >  1:735533/1‑46 (MQ=255)
                     ttGGCGATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAAt                                 >  1:3180410/1‑68 (MQ=255)
                     ttGGCGATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAAt                                 >  1:2345523/1‑68 (MQ=255)
                     ttGGCGATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATc                                >  1:2048388/1‑69 (MQ=255)
                     ttGGCGATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATc                                >  1:1725961/1‑69 (MQ=255)
                     ttGGCGATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATc                                >  1:1630093/1‑69 (MQ=255)
                     ttGGCGATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATc                                >  1:141558/1‑69 (MQ=255)
                     ttGGCGATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATc                                >  1:3342183/1‑69 (MQ=255)
                      tGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATcc                               >  1:1261082/1‑69 (MQ=255)
                      tGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATcc                               >  1:2819703/1‑69 (MQ=255)
                      tGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATcc                               >  1:2802205/1‑69 (MQ=255)
                       ggCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATccc                              <  1:3121038/69‑1 (MQ=255)
                       ggCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATccc                              <  1:683416/69‑1 (MQ=255)
                         cGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATccc                              <  1:3186869/67‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:367643/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:2303623/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:263074/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:2747951/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:1700740/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:1110479/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:1695473/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:118082/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:118958/68‑1 (MQ=255)
                          gATTTACGACTACATCAAATCACTCGAACAATCTGAATAACACTATGCCAATACACGATAAATCCCCt                            <  1:2041565/68‑1 (MQ=255)
                             ttATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTc                        >  1:2541728/1‑69 (MQ=255)
                             ttATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTc                        >  1:1996592/1‑69 (MQ=255)
                                  aCTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAg                   <  1:1862127/69‑1 (MQ=255)
                                  aCTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAg                   <  1:1858212/69‑1 (MQ=255)
                                        tCAAATCACTCGAACAATCTGAATAACCCTATGCCAATACACGATAAATCCCCt                            <  1:3077712/54‑1 (MQ=255)
                                              cACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTtgctgc        >  1:1299131/1‑68 (MQ=255)
                                              cACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCg       >  1:1048537/1‑69 (MQ=255)
                                              cACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCg       >  1:47055/1‑69 (MQ=255)
                                              cACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCg       >  1:1340720/1‑69 (MQ=255)
                                              cACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCg       >  1:3303851/1‑69 (MQ=255)
                                                     acaacCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGa  <  1:1101210/67‑1 (MQ=255)
                                                     acaacCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGa  <  1:2814393/67‑1 (MQ=255)
                                                         |                                                              
CGCCGCAAAGTACGGGCGCAGTTGGCGATTTATGACTACATCAAATCACTCGAACAACCTGAATAACCCTATGCCAATACACGATAAATCCCCTCGTCCGCAGGAGTTTGCTGCGGTCGA  >  minE/1657731‑1657850

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: