Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,688,979	T→C	25.3%	A22A (GCA→GCG)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,688,979	0	T	C	25.3%	76.7 / 22.9	67	A22A (GCA→GCG)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base C (6/11);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTT  >  minE/1688935‑1689047                                             |                                                                     ccGTTCAAATGCAGTGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:3381662/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:1041602/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:2156153/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:2236463/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:648952/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:1540226/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:1471215/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:3067695/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:3138779/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:1392527/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                               <  1:3198250/68‑1 (MQ=255)   gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCt                                                                   >  1:2389301/1‑46 (MQ=255)   gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                                            >  1:378631/1‑69 (MQ=255)   gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                                            >  1:1627588/1‑69 (MQ=255)   gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                                            >  1:1304247/1‑69 (MQ=255)   gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                                            >  1:2230820/1‑69 (MQ=255)   gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                                            >  1:2061063/1‑69 (MQ=255)      cAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAgg                                         >  1:1321602/1‑69 (MQ=255)      cAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAgg                                         >  1:1893135/1‑69 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:3415629/69‑1 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:3056761/69‑1 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:246518/69‑1 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:331068/69‑1 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:1847801/69‑1 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:209667/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:741060/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:73599/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:3183074/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:2950381/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:2805305/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:320665/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:270042/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:2656525/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:1017206/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:2103341/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:1654394/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:2175684/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                           >  1:2441184/1‑54 (MQ=255)                     gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAg                                          <  1:62848/53‑1 (MQ=255)                     gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt                          <  1:710337/69‑1 (MQ=255)                     gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt                          <  1:2822012/69‑1 (MQ=255)                     gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt                          <  1:1630082/69‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:1692665/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:2248216/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:2909458/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:2946993/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:679838/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                             <  1:3200204/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTa                         >  1:2860118/1‑69 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTa                         >  1:259869/1‑69 (MQ=255)                         tCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAAt                       >  1:303032/1‑68 (MQ=255)                          cATTATTAATGCTGTAAGTTGCTTTTTCTCCCCCGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATc                      >  1:1453251/1‑68 (MQ=255)                           attattAATGCTGTAAGTTGCTATTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTa                         >  1:2642593/1‑64 (MQ=255)                            ttattaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTg                                      <  1:1617780/50‑1 (MQ=255)                            ttattaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTg                                      <  1:666993/50‑1 (MQ=255)                             tattaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTc                   <  1:3107580/68‑1 (MQ=255)                             tattaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTc                   <  1:920545/68‑1 (MQ=255)                              attaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCAt                 >  1:1445877/1‑69 (MQ=255)                                   tGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCt             <  1:59696/68‑1 (MQ=255)                                   tGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCt             <  1:939100/68‑1 (MQ=255)                                   tGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCt             <  1:2638328/68‑1 (MQ=255)                                   tGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCt             <  1:2568826/68‑1 (MQ=255)                                   tGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCt             <  1:2443790/68‑1 (MQ=255)                                     cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTc           >  1:1735524/1‑68 (MQ=255)                                     cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTc           >  1:1434499/1‑68 (MQ=255)                                     cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCt          >  1:3306493/1‑69 (MQ=255)                                     cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCt          >  1:1363567/1‑69 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt  <  1:3167207/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt  <  1:3435730/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt  <  1:2801964/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt  <  1:1970577/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt  <  1:1809201/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt  <  1:737156/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt  <  1:1407937/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt  <  1:1259078/69‑1 (MQ=255)                                             |                                                                     CCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTT  >  minE/1688935‑1689047

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