Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,366,241	C→T	37.6%	K404K (AAG→AAA)	pepQ ←	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,366,241	0	C	T	37.6%	34.9 / 35.7	61	K404K (AAG→AAA)	pepQ	proline dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (27/11);  new base T (9/14);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTA  >  minE/2366175‑2366291                                                                   |                                                   cacGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGtgctt                                                  >  1:2716645/1‑69 (MQ=255) cacGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGtgctt                                                  >  1:670789/1‑69 (MQ=255) cacGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGtgctt                                                  >  1:181641/1‑69 (MQ=255) cacGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGtgctt                                                  >  1:2977741/1‑69 (MQ=255) cacGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGtgctt                                                  >  1:483059/1‑69 (MQ=255) cacGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGtgctt                                                  >  1:1256546/1‑69 (MQ=255) cacGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGtgctt                                                  >  1:698988/1‑69 (MQ=255) cacGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGtgctt                                                  >  1:3247073/1‑69 (MQ=255)     ttGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTg                                              <  1:390636/69‑1 (MQ=255)     ttGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTg                                              <  1:1103245/69‑1 (MQ=255)      tGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGa                                             <  1:1295037/69‑1 (MQ=255)      tGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGa                                             <  1:1637049/69‑1 (MQ=255)       gTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGaa                                            >  1:2851634/1‑69 (MQ=255)              aTACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGc                                        <  1:2692934/66‑1 (MQ=255)              aTACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGc                                        <  1:444875/66‑1 (MQ=255)              aTACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGc                                        <  1:2909033/66‑1 (MQ=255)              aTACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGc                                        <  1:1695670/66‑1 (MQ=255)              aTACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGc                                        <  1:1111010/66‑1 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCGATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:1514545/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:3329310/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:3325080/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:3048583/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:1059001/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:3294537/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:496873/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:2360776/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:1124194/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:742918/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:1551516/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:2225720/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  <  1:1561684/69‑1 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  >  1:559714/1‑69 (MQ=255)                 cGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCa                                  <  1:1689303/69‑1 (MQ=255)                       ccgccgAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCACGCCAc                            >  1:42905/1‑69 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:2310726/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:1418331/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:3063283/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:3048971/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:1989590/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:216225/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:2631599/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:251822/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:2461059/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:2429037/69‑1 (MQ=255)                                          cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:2275311/69‑1 (MQ=255)                                          cTTCAATTTTCTGACAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:2414311/69‑1 (MQ=255)                                           ttCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:622212/68‑1 (MQ=255)                                           ttCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa         <  1:2463305/68‑1 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:67257/1‑68 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:533071/1‑68 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:895384/1‑68 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:2539727/1‑68 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:3373294/1‑68 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:2070778/1‑68 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:2167911/1‑68 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:2260964/1‑68 (MQ=255)                                                  ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa  >  1:1048134/1‑68 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACTGTGCCAGTAGCGATTCGATGAAGTa  >  1:1582592/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGt   >  1:2037359/1‑56 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:2985627/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:3044562/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:3094823/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:382212/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:1592694/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:725986/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:1382995/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:761662/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:1212274/1‑57 (MQ=255)                                                             gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa  >  1:999200/1‑57 (MQ=255)                                                                   |                                                   CACGTTGTCTTCGATACGAATGCCGCCGAACGGTTTCAGTGCTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTA  >  minE/2366175‑2366291

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