Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	647,947	A→G	37.1%	V209V (GTA→GTG)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	647,947	0	A	G	37.1%	38.5 / 42.9	70	V209V (GTA→GTG)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (14/30);  new base G (13/13);  total (27/43)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

TGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAA  >  minE/647881‑648013                                                                   |                                                                   tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAg                                                                   >  1:2184568/1‑68 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:3322482/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:2151529/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:2392034/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  >  1:606134/1‑69 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:1737607/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:1555614/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:2850298/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:1522304/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:2900456/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:922613/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:1145068/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:1149406/69‑1 (MQ=255) tGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGc                                                                  <  1:3264087/69‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:2864050/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:597660/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:3386729/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:701460/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:1861711/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:1761963/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:2906822/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:1297621/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:1537979/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:742985/67‑1 (MQ=255)     ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTAACCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:275477/67‑1 (MQ=255)     gtggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:342843/67‑1 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTgg                                                             >  1:2622762/1‑68 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:258749/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:2484557/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:343062/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:2010221/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:451210/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:1102925/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:2030282/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:1248359/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:1467933/1‑69 (MQ=255)       ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                            >  1:1650299/1‑69 (MQ=255)                  attttGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                               <  1:756278/53‑1 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:448900/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:2915786/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:707418/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:1549525/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:3330291/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:1668907/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:1712456/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:1869964/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:1944714/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:243221/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:2323823/1‑68 (MQ=255)                               gcgcGATACCTTCACTACCCGCTCTGGCCGTGAACTGGCGCTGGAGTTGTATGTCGATCGTGGCaa                                       >  1:2596225/1‑66 (MQ=255)                               gcgcGATACCTTCACTACCCGCTATGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                     >  1:2069169/1‑68 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCTCGGCAACCtt                                   <  1:236491/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:1898983/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:986495/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:1132963/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:1383609/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:1470369/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:1867879/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:612270/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:2703997/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:2214016/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:2219895/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:391927/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:1083147/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                   <  1:3387716/63‑1 (MQ=255)                                      aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGGCGATCGCGGCAACCtt                                   <  1:2221256/63‑1 (MQ=255)                                                      cTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGa             >  1:2510132/1‑69 (MQ=255)                                                                aGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGaa   <  1:197465/69‑1 (MQ=255)                                                                 gTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGaaa  <  1:1191898/69‑1 (MQ=255)                                                                  tAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGaaa  <  1:721454/68‑1 (MQ=255)                                                                   |                                                                   TGCGCTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAA  >  minE/647881‑648013

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