Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	647,950	A→G	45.6%	A210A (GCA→GCG)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	647,950	0	A	G	45.6%	5.2 / 45.9	57	A210A (GCA→GCG)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base G (13/13);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

CTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAAC  >  minE/647885‑648016                                                                  |                                                                   ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:2864050/67‑1 (MQ=255) ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:597660/67‑1 (MQ=255) ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:1861711/67‑1 (MQ=255) ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:1761963/67‑1 (MQ=255) ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:3386729/67‑1 (MQ=255) ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:1537979/67‑1 (MQ=255) ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:701460/67‑1 (MQ=255) ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:1297621/67‑1 (MQ=255) ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:2906822/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:1102925/1‑69 (MQ=255)   ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:2030282/1‑69 (MQ=255)   ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:2010221/1‑69 (MQ=255)   ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:451210/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:707418/1‑68 (MQ=255)                           gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:2915786/1‑68 (MQ=255)                           gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:1549525/1‑68 (MQ=255)                           gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:243221/1‑68 (MQ=255)                           gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:3330291/1‑68 (MQ=255)                           gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:2323823/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)                           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:612270/63‑1 (MQ=255)                                  aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:1867879/63‑1 (MQ=255)                                  aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:2703997/63‑1 (MQ=255)                                  aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:391927/63‑1 (MQ=255)                                  aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:2214016/63‑1 (MQ=255)                                  aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:2219895/63‑1 (MQ=255)                                  aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:3387716/63‑1 (MQ=255)                                  aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:1083147/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)                                                                gCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAAc  >  1:858410/1‑69 (MQ=255)                                                                  |                                                                   CTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAAC  >  minE/647885‑648016

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