Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,931,277 A→G 37.7% A57A (GCT→GCC yjgR ← predicted ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,931,2770AG37.7% 25.3 / 32.5 53A57A (GCT→GCCyjgRpredicted ATPase
Reads supporting (aligned to +/- strand):  ref base A (25/8);  new base G (6/14);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 1.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAC  >  minE/2931212‑2931344
                                                                 |                                                                   
cTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGc                                                                   >  1:200951/1‑68 (MQ=255)
           tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                       <  1:2616172/69‑1 (MQ=255)
           tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                       <  1:2750616/69‑1 (MQ=255)
           tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                       <  1:3228539/69‑1 (MQ=255)
           tttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCgg                                                       <  1:1419529/69‑1 (MQ=255)
            ttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGc                                                      >  1:976880/1‑69 (MQ=255)
            ttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGc                                                      >  1:2612624/1‑69 (MQ=255)
            ttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGc                                                      >  1:714296/1‑69 (MQ=255)
            ttttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGc                                                      >  1:415406/1‑69 (MQ=255)
             tttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCa                                                     >  1:1560987/1‑69 (MQ=255)
             tttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCa                                                     >  1:386608/1‑69 (MQ=255)
             tttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCa                                                     >  1:1205911/1‑69 (MQ=255)
             tttCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCa                                                     >  1:1128361/1‑69 (MQ=255)
                    cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCATATCGGCCATAAACACCGGCACGCCGAt                                              >  1:5826/3‑69 (MQ=255)
                    cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAac                                                            >  1:1318643/3‑55 (MQ=255)
                    cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAACACCGGCACGCCGAt                                              >  1:439266/3‑69 (MQ=255)
                    cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAACACCGGCACGCCGAt                                              >  1:3028717/3‑69 (MQ=255)
                    cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAACACCGGCACGCCGAt                                              >  1:2686343/3‑69 (MQ=255)
                    cgccGTTCCTGCCTGCGCGATACCGGTCAGATCGCCTTTCACATCGGCCATAAACACCGGCACGCCGAt                                              >  1:2249141/3‑69 (MQ=255)
                     acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                             >  1:1095016/1‑69 (MQ=255)
                     acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                             >  1:586359/1‑69 (MQ=255)
                     acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                             >  1:616593/1‑69 (MQ=255)
                     acCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGAtt                                             >  1:1245713/1‑69 (MQ=255)
                               cctgcGCGACGCCGGTCAGATCGCCTCTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAAt                                    <  1:2580459/68‑1 (MQ=255)
                                   cgcgACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACt                                >  1:3219466/1‑68 (MQ=255)
                                   cgcgACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACt                                >  1:64393/1‑68 (MQ=255)
                                   cgcgACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTc                               >  1:1373521/1‑69 (MQ=255)
                                                           cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:19446/1‑69 (MQ=255)
                                                           cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:42645/1‑69 (MQ=255)
                                                           cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:1324241/1‑69 (MQ=255)
                                                           cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:1900123/1‑69 (MQ=255)
                                                           cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:1209585/1‑69 (MQ=255)
                                                           cacaTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGt       >  1:107343/1‑69 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTTACAGtttt    <  1:1170406/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:3273653/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:669974/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:806985/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:476710/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:868472/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:1595308/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2218559/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:1697166/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:1824191/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2687067/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2044463/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2602356/69‑1 (MQ=255)
                                                              aTCGGCCATAAACACCGGCACGCCGATTTCCGATAGAGACTCCGCCAGTTTTTGCAGCGTGACAGtttt    <  1:2353250/69‑1 (MQ=255)
                                                              aTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGttt     <  1:1564567/68‑1 (MQ=255)
                                                              aTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGttt     <  1:253691/68‑1 (MQ=255)
                                                              aTCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGttt     <  1:2084383/68‑1 (MQ=255)
                                                               tCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTa   >  1:2913470/1‑69 (MQ=255)
                                                                cAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAc  >  1:3379532/1‑69 (MQ=255)
                                                                cAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAc  >  1:202409/1‑69 (MQ=255)
                                                                 |                                                                   
CTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTAC  >  minE/2931212‑2931344

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