Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	959,697	A→G	35.9%	T203T (ACT→ACC)	pqqL ←	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	959,697	0	A	G	35.9%	45.0 / 44.5	75	T203T (ACT→ACC)	pqqL	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base A (29/19);  new base G (6/21);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGC  >  minE/959634‑959763                                                                |                                                                   attattTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATc                                                               <  1:414269/69‑1 (MQ=255) attattTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATc                                                               <  1:351485/69‑1 (MQ=255) attattTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATc                                                               <  1:2028042/69‑1 (MQ=255)      ttGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACGGTATCCATc                                                           >  1:183218/1‑68 (MQ=255)      ttGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACGGTATCCATc                                                           >  1:2096591/1‑68 (MQ=255)      ttGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACGGTATCCATCa                                                          >  1:985800/1‑69 (MQ=255)      ttGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACGGTATCCATCa                                                          >  1:290570/1‑69 (MQ=255)      ttGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACGGTATCCATCa                                                          >  1:2787294/1‑69 (MQ=255)      ttGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACGGTATCCATCa                                                          >  1:2156837/1‑69 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:1820620/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:1107297/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:2155415/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:2268424/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:2436898/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:2511421/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:3054549/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:3302661/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:737585/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:1150560/68‑1 (MQ=255)               aCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  <  1:1238191/68‑1 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATa                                                  >  1:1519198/1‑67 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAgg                                                >  1:1138373/1‑69 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAgg                                                >  1:2821939/1‑69 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAgg                                                >  1:2582310/1‑69 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAgg                                                >  1:2387391/1‑69 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAgg                                                >  1:332390/1‑69 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAgg                                                >  1:2008759/1‑69 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAg                                                 >  1:841332/1‑68 (MQ=255)                ccAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAg                                                 >  1:140984/1‑68 (MQ=255)                                      aaTTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTa                          <  1:1115877/69‑1 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCg                                                     >  1:1456876/1‑38 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:1737353/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:2678795/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:1006352/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:2847216/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:88685/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:2479221/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:664941/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:2922680/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:1705808/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:760003/1‑68 (MQ=255)                                          gtgCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTAttt                       >  1:779484/1‑68 (MQ=255)                                             ccGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGcc                    >  1:3228667/1‑68 (MQ=255)                                             ccGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAAATAAGGGTATTTGcc                    >  1:201127/1‑68 (MQ=255)                                              cGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGcc                    <  1:1871011/67‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCTGCAGGa             <  1:442429/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:434002/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:480424/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:56727/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:59212/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:3430431/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:1158408/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:3294555/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:2897560/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:3194389/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:1897750/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:2467191/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:2530408/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:2643668/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:2662161/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:127280/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:2809769/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:1779448/69‑1 (MQ=255)                                                   tgACCGTGGCGACGGTATCCATCAGGCCGACAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:2182531/69‑1 (MQ=255)                                                   tgACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCcagcag               >  1:1748228/1‑67 (MQ=255)                                                   tgACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCcagcag               >  1:837641/1‑67 (MQ=255)                                                   tgACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCcagcag               >  1:2027080/1‑67 (MQ=255)                                                   tgACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCcagcag               >  1:3203646/1‑67 (MQ=255)                                                   tgACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCcagca                >  1:3144222/1‑66 (MQ=255)                                                    gACCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:1212041/68‑1 (MQ=255)                                                     aCCGTGGCGACGGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGa             <  1:1724773/67‑1 (MQ=255)                                                           gCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGggcgg     >  1:575189/1‑69 (MQ=255)                                                              aCATTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGgcgc  <  1:409926/69‑1 (MQ=255)                                                              aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGgcgc  <  1:265266/69‑1 (MQ=255)                                                              aCAGTATCCATCAGGCCGATAGCTTCACGGTCTAAATTACGGGTATTCGCCAGCAGGAAAGGGCGgcgc  <  1:1777199/69‑1 (MQ=255)                                                                |                                                                   ATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGC  >  minE/959634‑959763

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