Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,688,958	G→A	47.6%	D29D (GAC→GAT)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,688,958	0	G	A	47.6%	3.7 / 52.5	63	D29D (GAC→GAT)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base G (18/15);  new base A (6/24);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 8.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTA  >  minE/1688900‑1689024                                                           |                                                                   gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:2428300/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:1097160/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:1147551/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:2491666/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:1319057/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:50047/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:3384049/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:2693846/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:1519399/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:2975988/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:3004207/65‑1 (MQ=255) gcgcCGCACGCTGGATCAACCCCAGATCGTCTGCCCCGTGCAAATGCAGTGTAATGCCATCatta                                                              <  1:1092794/65‑1 (MQ=255)  cgccgcACGCTGGATCAACCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCATCatta                                                              <  1:919042/64‑1 (MQ=255)                                    ccGTTCAAATGCAGTGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:3381662/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:3067695/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:3138779/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:3198250/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:648952/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:1041602/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:2236463/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:2156153/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:1540226/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:1471215/68‑1 (MQ=255)                                    ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                        <  1:1392527/68‑1 (MQ=255)                                      gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCt                                            >  1:2389301/1‑46 (MQ=255)                                      gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                     >  1:1627588/1‑69 (MQ=255)                                      gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                     >  1:2230820/1‑69 (MQ=255)                                      gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                     >  1:2061063/1‑69 (MQ=255)                                      gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                     >  1:1304247/1‑69 (MQ=255)                                      gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                     >  1:378631/1‑69 (MQ=255)                                         cAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAgg                  >  1:1893135/1‑69 (MQ=255)                                         cAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAgg                  >  1:1321602/1‑69 (MQ=255)                                                     aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:3056761/69‑1 (MQ=255)                                                     aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:331068/69‑1 (MQ=255)                                                     aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:3415629/69‑1 (MQ=255)                                                     aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:1847801/69‑1 (MQ=255)                                                     aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:246518/69‑1 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:73599/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:320665/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:741060/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:3183074/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:1017206/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:1654394/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:270042/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:2175684/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:209667/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:2441184/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:2656525/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:2103341/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:2805305/1‑54 (MQ=255)                                                      aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                    >  1:2950381/1‑54 (MQ=255)                                                        gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAg                   <  1:62848/53‑1 (MQ=255)                                                        gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt   <  1:2822012/69‑1 (MQ=255)                                                        gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt   <  1:710337/69‑1 (MQ=255)                                                        gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt   <  1:1630082/69‑1 (MQ=255)                                                         ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:2909458/65‑1 (MQ=255)                                                         ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:3200204/65‑1 (MQ=255)                                                         ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:2946993/65‑1 (MQ=255)                                                         ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:679838/65‑1 (MQ=255)                                                         ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:2248216/65‑1 (MQ=255)                                                         ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag      <  1:1692665/65‑1 (MQ=255)                                                         ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTa  >  1:2860118/1‑69 (MQ=255)                                                         ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTa  >  1:259869/1‑69 (MQ=255)                                                           |                                                                   GCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTA  >  minE/1688900‑1689024

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