Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	648,655	A→G	36.9%	T445T (ACA→ACG)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	648,655	0	A	G	36.9%	23.5 / 26.5	46	T445T (ACA→ACG)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (5/24);  new base G (10/7);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 7.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCA  >  minE/648590‑648713                                                                  |                                                           gCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCg                                                         >  1:3166668/1‑69 (MQ=255)       gAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATtgtg                                                   <  1:2445973/69‑1 (MQ=255)         aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCTGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:343352/68‑1 (MQ=255)         aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:2622113/68‑1 (MQ=255)         aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:232628/68‑1 (MQ=255)         aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:2317664/68‑1 (MQ=255)         aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:3042759/68‑1 (MQ=255)         aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:2196633/68‑1 (MQ=255)         aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:3339753/68‑1 (MQ=255)          gATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:197479/67‑1 (MQ=255)                 cGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:2850539/60‑1 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCgg                          >  1:3319778/1‑68 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:396484/1‑69 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:775401/1‑69 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:2515633/1‑69 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:2894755/1‑69 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:1695358/1‑69 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:1720884/1‑69 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:1805311/1‑69 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:729724/1‑69 (MQ=255)                                 ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:2102778/1‑69 (MQ=255)                                  cATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaa                        >  1:2400515/1‑69 (MQ=255)                                  cATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaa                        >  1:1662426/1‑69 (MQ=255)                                  cATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaa                        >  1:2692537/1‑69 (MQ=255)                                      tCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCg                    >  1:3227722/1‑69 (MQ=255)                                         gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:288400/68‑1 (MQ=255)                                         gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:552993/68‑1 (MQ=255)                                         gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:1545525/68‑1 (MQ=255)                                         gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:142647/68‑1 (MQ=255)                                         gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:2251609/68‑1 (MQ=255)                                         gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:1822937/68‑1 (MQ=255)                                           cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:3139501/69‑1 (MQ=255)                                           cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:2001885/69‑1 (MQ=255)                                           cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:454686/69‑1 (MQ=255)                                           cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:217851/69‑1 (MQ=255)                                           cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:2008610/69‑1 (MQ=255)                                           cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:2106231/69‑1 (MQ=255)                                            gTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:522255/65‑1 (MQ=255)                                                gTACAGCCAGTCCGGTACAGCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:2475079/64‑1 (MQ=255)                                                       cAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGAcc   <  1:2223693/69‑1 (MQ=255)                                                       cAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGAcc   <  1:2496548/69‑1 (MQ=255)                                                        aGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCa  <  1:3172020/69‑1 (MQ=255)                                                        aGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCa  <  1:2233516/69‑1 (MQ=255)                                                        aGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCa  <  1:3152137/69‑1 (MQ=255)                                                        aGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCAGGAAACCGAGCAGTACACCCTGACCa  <  1:2151700/69‑1 (MQ=255)                                                         gTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCa  <  1:547169/68‑1 (MQ=255)                                                                  |                                                           GCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCA  >  minE/648590‑648713

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑