Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,688,153	T→G	30.9%	M298L (ATG→CTG)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,688,153	0	T	G	30.9%	76.1 / 55.7	81	M298L (ATG→CTG)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base G (8/17);  total (33/48)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTG  >  minE/1688090‑1688220                                                                |                                                                    cgTGATAATCATTACCCAGCGCAGTTTTCGCCTCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCg                                                                >  1:2092537/1‑69 (MQ=255) cgTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGATGGTCGCCGCATCAATGATCATTTCCg                                                                >  1:3132370/1‑69 (MQ=255) cgTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCg                                                                >  1:71754/1‑69 (MQ=255) cgTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCg                                                                >  1:2077912/1‑69 (MQ=255) cgTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCg                                                                >  1:344698/1‑69 (MQ=255) cgTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCg                                                                >  1:2602366/1‑69 (MQ=255)        aTCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCt                                                          >  1:1407237/1‑68 (MQ=255)        aTCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTg                                                         >  1:478777/1‑69 (MQ=255)        aTCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTg                                                         >  1:438499/1‑69 (MQ=255)        aTCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTg                                                         >  1:401540/1‑69 (MQ=255)        aTCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTg                                                         >  1:2528272/1‑69 (MQ=255)           aTTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      >  1:1885612/1‑69 (MQ=255)           aTTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      >  1:201613/1‑69 (MQ=255)           aTTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      >  1:1733454/1‑69 (MQ=255)           aTTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      >  1:3107586/1‑69 (MQ=255)           aTTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      >  1:3365189/1‑69 (MQ=255)           aTTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      >  1:135726/1‑69 (MQ=255)           aTTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      >  1:1329720/1‑69 (MQ=255)           aTTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      >  1:640312/1‑69 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      <  1:2879381/53‑1 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      <  1:3216723/53‑1 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      <  1:2106989/53‑1 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      <  1:969179/53‑1 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      <  1:1872167/53‑1 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      <  1:1683164/53‑1 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTgcgc                                                      <  1:1076503/53‑1 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAg                                       >  1:190517/1‑68 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAg                                       >  1:2684328/1‑68 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAg                                       >  1:2655755/1‑68 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAg                                       >  1:3143179/1‑68 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAg                                       >  1:3227422/1‑68 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAg                                       >  1:1335343/1‑68 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAg                                       >  1:1228893/1‑68 (MQ=255)                           ttCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAg                                       >  1:1989502/1‑68 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:1367479/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:1495517/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:3372554/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:157935/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:887455/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:2930635/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:1680002/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:1901186/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:2863614/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:107009/68‑1 (MQ=255)                              gccgccCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGAcc                                    <  1:1735322/68‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:2499069/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:2023979/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:1747928/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:747689/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:707025/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:535297/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:419710/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:2483353/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:2732250/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:17164/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:1608533/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:3001629/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:2998109/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:2563556/69‑1 (MQ=255)                                   cccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:306553/69‑1 (MQ=255)                                    ccGGTGAGGGTCGCCGCATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:1606754/68‑1 (MQ=255)                                         gAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATTGGCAAg                           >  1:3324862/1‑66 (MQ=255)                                         gAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATTGGCAAGc                          >  1:1657347/1‑67 (MQ=255)                                         gAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATTGGCAAGc                          >  1:2778587/1‑67 (MQ=255)                                         gAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATTGGCAAGc                          >  1:786228/1‑67 (MQ=255)                                         gAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATTGGCAAGc                          >  1:2827475/1‑67 (MQ=255)                                         gAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATTGGCAAGc                          >  1:2303384/1‑67 (MQ=255)                                                gccgcATCATTGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:2667466/69‑1 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGTGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:1569081/69‑1 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGAccc                 <  1:1733233/69‑3 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:2974890/69‑1 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:2998082/69‑1 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:2555895/69‑1 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:2485602/69‑1 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:3262508/69‑1 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:3017788/69‑1 (MQ=255)                                                gccgcATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGc                 <  1:2733365/69‑1 (MQ=255)                                                 ccgcATCAATGATCAGTTCCGGTTTCTGCGCACTGGCATCAATCAGGCCATCGGc                              <  1:217616/55‑1 (MQ=39)                                                    cATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCGGAcc                                    <  1:2357818/46‑1 (MQ=255)                                                               aTTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAgtg  <  1:3202858/69‑1 (MQ=255)                                                               aTTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCATTCGGCATCAgtg  <  1:3299634/69‑1 (MQ=255)                                                                |                                                                    CGTGATAATCATTACCCAGCGCAGTTTTCGCCGCCCCGGTGAGGGTCGCCGCATCAATGATCATTTCCGGTTTCTGCGCACTGGCATCAATCAGACCATCGGCAAGCACCAGACGCCCTTCGGCATCAGTG  >  minE/1688090‑1688220

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