Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,601,372	A→G	40.6%	Q341Q (CAA→CAG)	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,601,372	0	A	G	40.6%	21.7 / 44.3	64	Q341Q (CAA→CAG)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base A (28/10);  new base G (12/14);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 3.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGAA  >  minE/2601305‑2601426                                                                    |                                                       ttGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                                                       <  1:1219168/69‑1 (MQ=255) ttGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                                                       <  1:1986448/69‑1 (MQ=255)   gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                                                       >  1:1205804/1‑67 (MQ=255)   gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                                                       >  1:3420056/1‑67 (MQ=255)   gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaaa                                                       >  1:270951/1‑67 (MQ=255)   gCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGaaaaacaa                                                        >  1:1289367/1‑66 (MQ=255)      aaaGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAg                                                  <  1:2330993/69‑1 (MQ=255)      aaaGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAg                                                  <  1:2535038/69‑1 (MQ=255)                     gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc                                   <  1:3356145/69‑1 (MQ=255)                     gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc                                   <  1:2650800/69‑1 (MQ=255)                     gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc                                   <  1:418675/69‑1 (MQ=255)                     gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc                                   <  1:656240/69‑1 (MQ=255)                     gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc                                   <  1:2830903/69‑1 (MQ=255)                     gTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATc                                   <  1:2926715/69‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:2720271/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:3098895/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:1965755/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:2011993/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:1131083/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:264382/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:239376/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:1165911/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:857767/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:1342423/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:610645/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:728344/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTCTACAGCAt                                    <  1:1366775/68‑1 (MQ=255)                     gTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGAAAAACAGAAACAGCACCTACACAGCAt                                    <  1:1137733/68‑1 (MQ=255)                                            tatCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTa            <  1:2274308/69‑1 (MQ=255)                                            tatCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTa            <  1:1664394/69‑1 (MQ=255)                                              tCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACg                            >  1:3253308/1‑51 (MQ=255)                                              tCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACg                            >  1:2442823/1‑51 (MQ=255)                                              tCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACg                            >  1:1269691/1‑51 (MQ=255)                                              tCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACg                            >  1:589283/1‑51 (MQ=255)                                              tCGCTTACAACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACg                            >  1:2266897/1‑51 (MQ=255)                                                gattACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACtt         >  1:172613/3‑68 (MQ=255)                                                gCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGt                                >  1:1790148/1‑45 (MQ=255)                                                gCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACtt         >  1:3060057/1‑68 (MQ=255)                                                gCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACtt         >  1:1841755/1‑68 (MQ=255)                                                gCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTc        >  1:601948/1‑69 (MQ=255)                                                gCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTc        >  1:3121073/1‑69 (MQ=255)                                                gCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTc        >  1:991172/1‑69 (MQ=255)                                                gCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTc        >  1:1740937/1‑69 (MQ=255)                                                gCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTc        >  1:1799758/1‑69 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAAc                         >  1:1844111/1‑48 (MQ=38)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:1593690/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:1547692/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:1369130/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:3097433/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:487586/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:1948693/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:2762908/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:2165388/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:2185909/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:2272660/1‑68 (MQ=255)                                                    actacAGCGAAAAACAGAAACAGCACCTCTACAGCATCAGCGATGAACAACTGCGTCCGTACTTCCCg     >  1:2319263/1‑68 (MQ=255)                                                       acAGCGAAAAACAACAACAGCACCTCTACAGCATCAGTGACg                            >  1:1789537/1‑42 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGt                                >  1:2625086/1‑38 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGa                           >  1:1994049/1‑43 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGa                           >  1:2275809/1‑43 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGa                           >  1:1673849/1‑43 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:1308265/1‑68 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:1077624/1‑68 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:2356719/1‑68 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:1651804/1‑68 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:165716/1‑68 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:3370148/1‑68 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:1852826/1‑68 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:1941860/1‑68 (MQ=255)                                                       acAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGaa  >  1:2497870/1‑68 (MQ=255)                                                                    |                                                       TTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGCATCAGTGACGAACAGCTGCGTCCGTACTTCCCGGAA  >  minE/2601305‑2601426

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