Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,286,931 C→T 62.8% P152P (CCG→CCA 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,9310CT62.8% 31.1 / 35.5 59P152P (CCG→CCAargCN‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (17/5);  new base T (22/15);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAG  >  minE/2286865‑2286991
                                                                  |                                                            
gAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCgg                                                            >  1:1103126/1‑69 (MQ=255)
gAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCgg                                                            >  1:3186598/1‑69 (MQ=255)
gAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCgg                                                            >  1:3170284/1‑69 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATaaag                                               >  1:1424902/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:1713827/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:2460240/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:718995/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:872402/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:2497330/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:2571436/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:268451/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:641308/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:1872846/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:281663/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:1805164/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag                                               >  1:2984344/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:3428406/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:1112898/1‑68 (MQ=255)
             tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa                                                >  1:1377326/1‑68 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:245206/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:3234805/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:1176309/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:1068179/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:2929654/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:950053/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:544889/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:1886862/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:688170/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:271149/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:532208/69‑1 (MQ=255)
                  atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc                                          <  1:2115201/69‑1 (MQ=255)
                     aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc                                       >  1:1972490/1‑69 (MQ=255)
                     aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc                                       >  1:3244726/1‑69 (MQ=255)
                     aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc                                       >  1:1071325/1‑69 (MQ=255)
                        caaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCa                                                 >  1:2194924/1‑56 (MQ=255)
                        caaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCa                                                 >  1:918793/1‑56 (MQ=255)
                               ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt                              >  1:621768/1‑68 (MQ=255)
                               ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt                              >  1:995849/1‑68 (MQ=255)
                               ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt                              >  1:396870/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)
                                                      cGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCtttt                                   <  1:1813406/40‑1 (MQ=255)
                                                          tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGAcc      >  1:38202/1‑65 (MQ=255)
                                                          tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAg  >  1:223440/1‑69 (MQ=255)
                                                          tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAg  >  1:940097/1‑69 (MQ=255)
                                                          tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAg  >  1:1699649/1‑69 (MQ=255)
                                                          tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCAGCCAGACCGTAg  >  1:2363374/1‑69 (MQ=255)
                                                                  |                                                            
GAGGTCAAGAAGATCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAG  >  minE/2286865‑2286991

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: