Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	323,098	G→A	44.3%	P75P (CCG→CCA)	recR →	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	323,098	0	G	A	44.3%	23.8 / 76.7	97	P75P (CCG→CCA)	recR	gap repair protein
Reads supporting (aligned to +/- strand):  ref base G (46/8);  new base A (26/17);  total (72/25)
Fisher's exact test for biased strand distribution p-value = 9.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAG  >  minE/323032‑323161                                                                   |                                                                aaTCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCcgcg                                                               >  1:1647409/1‑69 (MQ=255) aaTCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCcgcg                                                               >  1:2336350/1‑69 (MQ=255)  aTCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGcgt                                                              >  1:1268045/1‑69 (MQ=255)  aTCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGcgt                                                              >  1:3085766/1‑69 (MQ=255)  aTCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGcgt                                                              >  1:144082/1‑69 (MQ=255)     ggCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGCATCCGcgtc                                                             >  1:540613/1‑67 (MQ=255)     ggCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGcgtc                                                             >  1:1556856/1‑67 (MQ=255)     ggCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGcgtc                                                             >  1:423658/1‑67 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2170168/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:1749729/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:1934123/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:3244092/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2004677/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:908049/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2112311/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2126063/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2137161/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:1046638/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2407004/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2426917/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:300628/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2451225/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:698221/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:2516222/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:270141/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:929104/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:1441330/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:3351151/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:3421847/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:150686/1‑69 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGa                                                      >  1:914585/1‑68 (MQ=255)           tGCGCCGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGa                                                      >  1:1942875/1‑68 (MQ=255)           tGCGCCGATTGCCGCACCTTCACAGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:3421598/1‑69 (MQ=255)           tGCGCAGATTGCCGCACCTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGaa                                                     >  1:1035268/1‑69 (MQ=255)            gcgcCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGaaaa                                                    >  1:1506418/1‑69 (MQ=255)            gcgcCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGaaaa                                                    >  1:2404698/1‑69 (MQ=255)            gcgcCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGaaaa                                                    >  1:2112757/1‑69 (MQ=255)            gcgcCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGaaaa                                                    >  1:1462308/1‑69 (MQ=255)            gcgcCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGaaaa                                                    >  1:3109941/1‑69 (MQ=255)                     gCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaaa                                           <  1:2187560/69‑1 (MQ=255)                     gCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaaa                                           <  1:2627081/69‑1 (MQ=255)                     gCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaaa                                           >  1:134769/1‑69 (MQ=255)                     gCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaaa                                           >  1:1566551/1‑69 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCcgc                                                                >  1:3284284/1‑40 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtg                                   >  1:847856/1‑69 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtg                                   >  1:2678981/1‑69 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtg                                   >  1:1864593/1‑69 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtg                                   >  1:399183/1‑69 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtg                                   >  1:1466860/1‑69 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtg                                   >  1:2269961/1‑69 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtg                                   >  1:2112763/1‑69 (MQ=255)                             ttCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgt                                    >  1:3235823/1‑68 (MQ=255)                                 ccGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtggtg                                <  1:1258855/68‑1 (MQ=255)                                 ccGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCgtggtg                                <  1:1477318/68‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:3340012/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:355832/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:648024/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:376987/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:3010228/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:682372/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:1211619/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:964626/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:961734/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:2185361/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:2146028/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:2079518/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:1946894/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:912258/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:1638393/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGagag                           <  1:958347/69‑1 (MQ=255)                                     aCAGGAAGTCTGTAACATCTGTTCGAATCCACGTCGTCAGGAAAACGGTCAAATCAGCGTGGTGGagag                           <  1:2435983/69‑1 (MQ=255)                                              cTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCcggcgg                   >  1:118610/1‑68 (MQ=255)                                              cTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGa                  >  1:2541465/1‑69 (MQ=255)                                              cTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGa                  >  1:1510425/1‑69 (MQ=255)                                              cTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGa                  >  1:85470/1‑69 (MQ=255)                                              cTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGa                  >  1:299692/1‑69 (MQ=255)                                              cTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGa                  >  1:2862224/1‑69 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:620384/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:996078/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:997027/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:3405177/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:120304/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:3274394/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:2898526/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:1241978/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:2471727/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:1571520/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:364551/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCaa                                            >  1:1798106/1‑42 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCa                                             >  1:1389478/1‑41 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGga                              >  1:2236085/1‑56 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGga                              <  1:1183469/56‑1 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGga                              >  1:2074268/1‑56 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGga                              <  1:411714/56‑1 (MQ=255)                                               tGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGAc                 <  1:2817191/69‑1 (MQ=255)                                                     aTCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTAc           >  1:2434067/1‑69 (MQ=255)                                                     aTCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTAc           >  1:2613790/1‑69 (MQ=255)                                                     aTCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTAc           >  1:469566/1‑69 (MQ=255)                                                        tGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGcc        >  1:1236847/1‑69 (MQ=255)                                                               aTCCGCGGCGTCAGGAAAACGGTCAAATTTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAg  <  1:1275595/68‑1 (MQ=255)                                                                   |                                                                AATCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAGTCTGTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAG  >  minE/323032‑323161

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: ‑