Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,517,661 A→C 42.9% G178G (GGT→GGG folC ← bifunctional folylpolyglutamate synthase and dihydrofolate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,517,6610AC42.9% 10.9 / 74.2 70G178G (GGT→GGGfolCbifunctional folylpolyglutamate synthase and dihydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base A (10/30);  new base C (27/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 4.08e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

CAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATCGGCGTCGACAA  >  minE/1517597‑1517722
                                                                |                                                             
cAATTGCCGGTTTTTCGCTGCGAAAGATACCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGCCcca                                                           <  1:1383689/69‑1 (MQ=255)
cAATTGCCGGTTTTTCGCTGCGAAAGATACCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGCCcca                                                           <  1:3204737/69‑1 (MQ=255)
cAATTGCCGGTTTTTCGCTGCGAAAGATACCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGCCcca                                                           <  1:739248/69‑1 (MQ=255)
         gTTTTTCGCTGCGGTAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:1387789/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:3001886/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:2080945/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:2220418/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:2447068/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:2883466/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:1618575/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:593902/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:1551327/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:679775/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:680834/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:877926/69‑1 (MQ=255)
         gTTTTTCGCTGCGGAAGATGCATGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCgg                                                  <  1:1899391/69‑1 (MQ=255)
                       aaGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAgcgaa                                    <  1:2740467/69‑3 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCATTCGGTATGATCCAGCGCaa                                   <  1:1546370/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:2559700/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:936460/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:2724197/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:255258/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:2549084/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:2529804/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:752175/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:1981374/69‑1 (MQ=255)
                        aGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCaa                                   <  1:1484902/69‑1 (MQ=255)
                         gATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAAt                                  >  1:2018014/1‑69 (MQ=255)
                         gATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAAt                                  >  1:863808/1‑69 (MQ=255)
                         gATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAAt                                  >  1:625475/1‑69 (MQ=255)
                         gATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAAt                                  >  1:3215438/1‑69 (MQ=255)
                         gATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAAt                                  >  1:1686984/1‑69 (MQ=255)
                         gATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAAt                                  >  1:1827224/1‑69 (MQ=255)
                               tGCTTTCTCGCGACCAATACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGt                            >  1:2989444/1‑69 (MQ=255)
                               tGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTgg                             >  1:3023687/1‑68 (MQ=255)
                               tGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGt                            >  1:1594387/1‑69 (MQ=255)
                               tGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGt                            >  1:2596705/1‑69 (MQ=255)
                               tGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGt                            >  1:980727/1‑69 (MQ=255)
                                   ttCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATg                                              <  1:261384/47‑1 (MQ=255)
                                   ttCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTtac                         <  1:1415302/68‑1 (MQ=255)
                                   ttCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTtac                         <  1:1621928/68‑1 (MQ=255)
                                      tcGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTtactac                      <  1:1886334/68‑1 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCTCAATACTGGTTACTACCGCGACATc             >  1:2511531/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATAc                                >  1:1941898/1‑49 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATAc                                >  1:2033952/1‑49 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACAt              >  1:1939426/1‑67 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:853237/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:2214140/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:828071/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:1452718/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:1388487/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:2377907/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:1511991/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:1524887/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:660863/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:985904/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:3385844/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:471089/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:3424142/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:3338583/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:1594877/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:3142613/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:1926984/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:269594/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:2664435/1‑68 (MQ=255)
                                               aTACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:2589828/1‑68 (MQ=255)
                                               aTACTATCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATc             >  1:2682472/1‑68 (MQ=255)
                                                tACTTTCGCGATCTGGCCCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATCgg           >  1:1305285/1‑69 (MQ=255)
                                                         gATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATCGGCGTCGACaa  <  1:85657/69‑1 (MQ=255)
                                                         gATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATCGGCGTCGACaa  <  1:171761/69‑1 (MQ=255)
                                                                |                                                             
CAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCTTTCTCGCGACCAATACTTTCGCGATCTGGACCCAGCCAGTCGGTATGATCCAGCGCAATACTGGTTACTACCGCGACATCGGCGTCGACAA  >  minE/1517597‑1517722

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: