Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,688,994	G→A	20.4%	R17R (CGC→CGT)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,688,994	0	G	A	20.4%	113.4 / 19.5	78	R17R (CGC→CGT)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base G (26/36);  new base A (5/11);  total (31/47)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTC  >  minE/1688935‑1689058                                                            |                                                                 ccGTTCAAATGCAGTGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                                          <  1:3381662/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:3067695/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                                          <  1:3138779/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                                          <  1:3198250/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                                          <  1:1540226/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                                          <  1:1471215/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                                          <  1:648952/68‑1 (MQ=255) ccGTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTc                                                          <  1:1392527/68‑1 (MQ=255)   gTTCAAATGCAGGGTAATGCCATCATTATTAATGCTGTAAGTCGCTTTTTCTCCCCAACGTGCGTCGGc                                                       >  1:378631/1‑69 (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)      cAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAgg                                                    >  1:1321602/1‑69 (MQ=255)      cAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAgg                                                    >  1:1893135/1‑69 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:246518/69‑1 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:3415629/69‑1 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:3056761/69‑1 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:1847801/69‑1 (MQ=255)                  aaTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:331068/69‑1 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:3183074/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:73599/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:2950381/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:2805305/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:270042/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:741060/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:1017206/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:320665/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:2656525/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:2441184/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:1654394/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:209667/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:2103341/1‑54 (MQ=255)                   aTGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCa                                                      >  1:2175684/1‑54 (MQ=255)                     gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAg                                                     <  1:62848/53‑1 (MQ=255)                     gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt                                     <  1:710337/69‑1 (MQ=255)                     gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt                                     <  1:1630082/69‑1 (MQ=255)                     gCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGt                                     <  1:2822012/69‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:679838/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:2909458/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:2946993/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:3200204/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:2248216/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTagagag                                        <  1:1692665/65‑1 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTa                                    >  1:259869/1‑69 (MQ=255)                      ccGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTa                                    >  1:2860118/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:666993/50‑1 (MQ=255)                            ttattaATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTg                                                 <  1:1617780/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:2443790/68‑1 (MQ=255)                                   tGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCt                        <  1:2568826/68‑1 (MQ=255)                                   tGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCt                        <  1:2638328/68‑1 (MQ=255)                                     cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTc                      >  1:1735524/1‑68 (MQ=255)                                     cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTc                      >  1:1434499/1‑68 (MQ=255)                                     cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCt                     >  1:1363567/1‑69 (MQ=255)                                     cTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCt                     >  1:3306493/1‑69 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt             <  1:1407937/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt             <  1:1259078/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt             <  1:737156/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt             <  1:2801964/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt             <  1:1970577/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt             <  1:3435730/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt             <  1:1809201/69‑1 (MQ=255)                                             tGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCAtttt             <  1:3167207/69‑1 (MQ=255)                                              gCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTc                              <  1:2097572/51‑1 (MQ=255)                                                          cAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCcttc  <  1:656897/67‑1 (MQ=255)                                                          cAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCcttc  <  1:2882544/67‑1 (MQ=255)                                                          cAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCcttc  <  1:626540/67‑1 (MQ=255)                                                            |                                                                 CCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCGCGCGTCGGCAGGTTGGGTAGAGAGGGTAATCTTCATCGCTTCTGTCATTTTAGTTATCCTTC  >  minE/1688935‑1689058

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