Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	649,279	G→A	48.4%	A653A (GCG→GCA)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	649,279	0	G	A	48.4%	‑0.1 / 53.3	62	A653A (GCG→GCA)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base G (22/10);  new base A (14/16);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGAT  >  minE/649212‑649342                                                                    |                                                                tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:631093/1‑69 (MQ=255) tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:61049/1‑69 (MQ=255) tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3365004/1‑69 (MQ=255) tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3184040/1‑69 (MQ=255) tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:870526/1‑69 (MQ=255) tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2833225/1‑69 (MQ=255) tcatcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:162298/1‑69 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:717633/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:1889601/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2159660/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2202309/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2270281/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2418521/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2684704/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2699246/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:2773834/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:1039972/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3146874/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3196030/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3233232/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:3425462/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:410802/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:686701/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:1640649/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:1679176/1‑67 (MQ=255)   atcGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCgg                                                                >  1:167989/1‑67 (MQ=255)         ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        <  1:1787745/69‑1 (MQ=255)         ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        >  1:1256431/1‑69 (MQ=255)         ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        >  1:3089956/1‑69 (MQ=255)         ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        >  1:1384253/1‑69 (MQ=255)         ccGATTGCTATTGCCGAAGTACGCGAAGCACTCAATCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                        >  1:1746206/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCCCTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2375520/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGTCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:777253/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:1439487/1‑56 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:2977839/1‑56 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:2860771/1‑56 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:2843449/1‑56 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCt                                                         >  1:1914475/1‑56 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1687972/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1081225/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2371486/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2697056/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1731785/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2700809/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1939681/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:511292/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:2929220/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:3314504/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1781457/1‑69 (MQ=255)                     gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGAGACTGAACTGGCGGATGAGCt                                                         >  1:2751336/1‑56 (MQ=255)                     gCCGAAGTACGCGAAGCACTAACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                            >  1:1170248/1‑69 (MQ=255)                        gAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAc                                         <  1:1209284/69‑1 (MQ=255)                           gTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAcgcg                                      <  1:1241995/69‑1 (MQ=255)                           gTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAcgcg                                      <  1:72064/69‑1 (MQ=255)                              cgcgAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAAt                                   >  1:2025204/1‑69 (MQ=255)                              cgcgAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAAt                                   >  1:2768018/1‑69 (MQ=255)                                 gAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTAc                                <  1:3148309/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:700612/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2572033/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:1000810/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:411763/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:3156675/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:529008/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:3154873/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2879092/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2122453/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:742896/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2822558/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:27993/69‑1 (MQ=255)                                    gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2684228/69‑1 (MQ=255)                                     cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2983182/68‑1 (MQ=255)                                     cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:1350446/68‑1 (MQ=255)                                     cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                             <  1:2602838/68‑1 (MQ=255)                                            tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGTTATTTACAACGCGAATTACCAGAGCGAGTa                     <  1:1477375/69‑1 (MQ=255)                                            tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                     <  1:1304859/69‑1 (MQ=255)                                            tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                     <  1:1753305/69‑1 (MQ=255)                                            tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                     <  1:2162719/69‑1 (MQ=255)                                             cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                     >  1:1364016/1‑68 (MQ=255)                                             cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:890897/1‑69 (MQ=255)                                             cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:954750/1‑69 (MQ=255)                                             cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:3324481/1‑69 (MQ=255)                                             cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:2700229/1‑69 (MQ=255)                                             cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:2444096/1‑69 (MQ=255)                                             cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                    >  1:1752492/1‑69 (MQ=255)                                                               cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat  >  1:2892202/1‑69 (MQ=255)                                                               cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat  >  1:2880047/1‑69 (MQ=255)                                                               cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat  >  1:2843020/1‑69 (MQ=255)                                                               cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGa   <  1:1300499/68‑1 (MQ=255)                                                                    |                                                                TCATCGACCCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGAT  >  minE/649212‑649342

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