Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,366,274	T→A	43.7%	L393L (CTA→CTT)	pepQ ←	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,366,274	0	T	A	43.7%	22.9 / 75.3	71	L393L (CTA→CTT)	pepQ	proline dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (27/13);  new base A (9/22);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACGGG  >  minE/2366216‑2366338                                                           |                                                                 cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:251822/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:216225/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:2275311/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:2310726/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:2429037/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:2461059/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:1989590/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:2631599/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:3048971/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:1418331/69‑1 (MQ=255) cTTCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:3063283/69‑1 (MQ=255) cTTCAATTTTCTGACAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:2414311/69‑1 (MQ=255)  ttCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:2463305/68‑1 (MQ=255)  ttCAATTTTCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGa                                                        <  1:622212/68‑1 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:2539727/1‑68 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:1048134/1‑68 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:3373294/1‑68 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:533071/1‑68 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:67257/1‑68 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:2260964/1‑68 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:2167911/1‑68 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:895384/1‑68 (MQ=255)         ttCTGCCAGTTGAAGTGTTTGCTGAACTGCCCTTCACGCCACGGTGCCAGAAGCGATTCGATGAAGTa                                                 >  1:2070778/1‑68 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACTGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:1582592/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGt                                                  >  1:2037359/1‑56 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:3044562/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:999200/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:761662/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:725986/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:1212274/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:382212/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:1382995/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:3094823/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:2985627/1‑57 (MQ=255)                    gAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTa                                                 >  1:1592694/1‑57 (MQ=255)                            tgctGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATg                                                      >  1:1137633/1‑44 (MQ=255)                            tgctGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGtt                              >  1:452053/1‑68 (MQ=255)                            tgctGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGtt                              >  1:34605/1‑68 (MQ=255)                            tgctGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGtt                              >  1:2681551/1‑68 (MQ=255)                            tgctGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGtt                              >  1:2696750/1‑68 (MQ=255)                                gAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAacac                         <  1:2208743/69‑1 (MQ=255)                                gAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAacac                         <  1:453729/69‑1 (MQ=255)                                gAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAacac                         <  1:3395203/69‑1 (MQ=255)                                gAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAacac                         <  1:205077/69‑1 (MQ=255)                                gAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAacac                         <  1:1399478/69‑1 (MQ=255)                                gAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAacac                         <  1:773574/69‑1 (MQ=255)                                gAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAacac                         <  1:1046654/69‑1 (MQ=255)                                gAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAacac                         <  1:2586729/69‑1 (MQ=255)                                           cACTCCACGTTGCCAGTAGAGATTAGATGAAGTAGATACACGGTTCGATGGTTAACACCAAGCCCGGc               >  1:3317573/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:1860337/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:2961619/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:3122861/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:2888587/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:2876980/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               <  1:2477988/68‑1 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:356526/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:2311158/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:463792/1‑68 (MQ=255)                                           cACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               >  1:975859/1‑68 (MQ=255)                                                 aCGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGaa        <  1:2114346/69‑1 (MQ=255)                                                 aCGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGaa        <  1:450482/69‑1 (MQ=255)                                                  cGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGaa        <  1:1370544/68‑1 (MQ=255)                                                       ccAGAAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACggg  <  1:532530/69‑1 (MQ=255)                                                       ccAGAAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACggg  <  1:2315008/69‑1 (MQ=255)                                                       ccAGAAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACggg  <  1:2035174/69‑1 (MQ=255)                                                       ccAGAAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACggg  <  1:2062100/69‑1 (MQ=255)                                                       ccAGAAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACggg  <  1:1995524/69‑1 (MQ=255)                                                       ccAGAAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACggg  <  1:2117234/69‑1 (MQ=255)                                                       ccAGAAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACggg  <  1:3020432/69‑1 (MQ=255)                                                       ccAGAAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACggg  <  1:1797654/69‑1 (MQ=255)                                                         aGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGc               <  1:1613071/54‑1 (MQ=255)                                                           |                                                                 CTTCAATTTTCTGCCAGTTGAAGTGCTTGCTGAACTGCCCTTCACGCCACGGTGCCAGTAGCGATTCGATGAAGTAGATACCCGGTTCGATGGTTAACACCATGCCCGGCTGGAGAATACGGG  >  minE/2366216‑2366338

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