Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,286,913 C→T 68.6% A158A (GCG→GCA argC ← N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,286,9130CT68.6% 43.5 / 25.1 54A158A (GCG→GCAargCN‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (13/4);  new base T (22/15);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCACCGGCCACTGATTGAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTAT  >  minE/2286847‑2286965
                                                                  |                                                    
gATCACCGGCCACTGATTGAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTgccgc                                                    >  1:1078411/1‑69 (MQ=255)
gATCACCGGCCACTGATTGAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTgccgc                                                    >  1:2841435/1‑69 (MQ=255)
gATCACCGGCCACTGATTGAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTgccgc                                                    >  1:592441/1‑69 (MQ=255)
gATCACCGGCCACTGATTGAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTgccgc                                                    >  1:809403/1‑69 (MQ=255)
                  gAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCgg                                  >  1:1103126/1‑69 (MQ=255)
                  gAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCgg                                  >  1:3170284/1‑69 (MQ=255)
                  gAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCgg                                  >  1:3186598/1‑69 (MQ=255)
                        aagaagATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGAt                                            >  1:398776/1‑53 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATaaag                     >  1:1424902/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:1872846/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:718995/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:268451/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:2571436/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:2497330/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:2460240/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:641308/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:872402/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:3428406/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:281663/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:1805164/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:1713827/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:1647529/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:1572715/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:1551126/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                     >  1:2984344/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa                      >  1:482846/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa                      >  1:258202/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa                      >  1:3358785/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa                      >  1:1377326/1‑68 (MQ=255)
                               tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa                      >  1:1112898/1‑68 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:3234805/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:544889/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:532208/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:3117192/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:3095463/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:3078080/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:2929654/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:688170/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:271149/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:245206/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:2115201/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:950053/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:1886862/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:1176309/69‑1 (MQ=255)
                                    atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                <  1:1068179/69‑1 (MQ=255)
                                       aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc             >  1:1071325/1‑69 (MQ=255)
                                       aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc             >  1:3244726/1‑69 (MQ=255)
                                       aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc             >  1:1972490/1‑69 (MQ=255)
                                          caaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCa                       >  1:2194924/1‑56 (MQ=255)
                                          caaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCa                       >  1:918793/1‑56 (MQ=255)
                                                 ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt    >  1:396870/1‑68 (MQ=255)
                                                 ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt    >  1:621768/1‑68 (MQ=255)
                                                 ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt    >  1:995849/1‑68 (MQ=255)
                                                 atCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt    >  1:2722005/2‑68 (MQ=255)
                                                  tCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATttat  <  1:2814048/69‑1 (MQ=255)
                                                  tCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATttat  <  1:2787722/69‑1 (MQ=255)
                                                      cgcCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCtttt         <  1:534484/58‑1 (MQ=255)
                                                      cgcCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCtttt         <  1:2001921/58‑1 (MQ=255)
                                                                  |                                                    
GATCACCGGCCACTGATTGAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTAT  >  minE/2286847‑2286965

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: