Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	647,833	C→T	79.7%	R171R (CGC→CGT)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	647,833	0	C	T	79.7%	58.4 / 8.9	44	R171R (CGC→CGT)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base C (2/7);  new base T (20/15);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGCTACCTGTTTG  >  minE/647766‑647882                                                                    |                                                  aaaTCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGcc                                                  <  1:820234/69‑1 (MQ=255)            cccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGt                                           >  1:25797/1‑65 (MQ=255)            cccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAg                                       >  1:1956816/1‑69 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:1102814/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:906952/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:814729/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:643951/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:498258/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:2921744/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:2583520/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:240108/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:1938750/56‑1 (MQ=255)             ccTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:1230292/56‑1 (MQ=255)             cattCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:1815941/54‑1 (MQ=255)              cTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGc                                                   <  1:3209380/55‑1 (MQ=255)                tCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGc                                   <  1:1705924/69‑1 (MQ=255)                tCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGc                                   <  1:2955244/69‑1 (MQ=255)                tCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGc                                   <  1:2246093/69‑1 (MQ=255)                tCCTGCTTTCCAATGGTAACCGCGTTCCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGc                                   <  1:2618291/69‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:3433876/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:429611/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:465663/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:521732/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:2886721/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:272652/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:1558120/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:2017471/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:1848547/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:1895332/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:1998910/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:1994211/68‑1 (MQ=255)                                    cgcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:1599459/68‑1 (MQ=255)                                     gcgTTGCGCAAGGTGAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:2800216/67‑1 (MQ=255)                                     gcgTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGt              <  1:3168895/69‑1 (MQ=255)                                     gcgTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGt              <  1:1403547/69‑1 (MQ=255)                                     gcgTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGt              <  1:781137/69‑1 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGTGTACAGTGGCAGGACCCGTTCCCGAAAcc                <  1:2927220/53‑1 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGa                               >  1:653831/3‑40 (MQ=37)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGAccc                            >  1:3303540/3‑43 (MQ=37)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGt              >  1:1607075/3‑57 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTg      >  1:1399207/3‑65 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGttt   >  1:1332252/3‑68 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGttt   >  1:254699/3‑68 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:630514/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:1354543/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:1067385/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:3404908/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:3108088/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:3105006/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:3046500/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:1472889/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:14860/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:2030630/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:1925336/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:1862712/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:177317/3‑69 (MQ=255)                                                 gtgAACTGGAAAATGGACGTCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGTTACCTGTTTg  >  1:1747362/3‑69 (MQ=255)                                                                    |                                                  AAATCAAATATCCCTTCCTGCTTTCCAATGGTAACCGCGTTGCGCAAGGCGAACTGGAAAACGGACGCCATTGGGTACAGTGGCAGGACCCGTTCCCGAAACCGTGCTACCTGTTTG  >  minE/647766‑647882

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