Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	166,699	C→T	67.3%	R259R (CGG→CGA)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	166,699	0	C	T	67.3%	42.5 / 29.4	58	R259R (CGG→CGA)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (12/7);  new base T (16/23);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAC  >  minE/166643‑166763                                                         |                                                                 aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTccgctcgc                                                               >  1:2854763/1‑60 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:3339450/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:1695422/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:2056906/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:2180135/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:2202569/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:220358/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:220851/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:84409/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:1512681/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:2230281/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:2781671/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:1363299/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:433250/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:35002/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCagag                                                      >  1:3164309/1‑69 (MQ=255) aGATGCAGGGCAAAGCGAATGCGCCACAATATATGCAGGCATTCGTTTAATTCCGCTCGCTCCGCAGa                                                       >  1:2369880/1‑68 (MQ=255)           cAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTaagaag                                             >  1:3081394/1‑68 (MQ=255)                cGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCaaa                                        >  1:3357239/1‑68 (MQ=255)                    tACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAAccct                                    >  1:1442123/1‑67 (MQ=255)                    tACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCg                                    >  1:57633/1‑68 (MQ=255)                             atatatGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCa                           >  1:1842785/1‑68 (MQ=255)                             atatatGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCa                           >  1:2415599/1‑68 (MQ=255)                               atatGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATcc                        >  1:532460/1‑69 (MQ=255)                                tatGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCa                       >  1:2608257/1‑69 (MQ=255)                                    gaGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGa                   <  1:2850061/68‑1 (MQ=255)                                    cAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGa                   <  1:1127503/69‑1 (MQ=255)                                    cAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGa                   <  1:1341659/69‑1 (MQ=255)                                    cAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGa                   <  1:2780181/69‑1 (MQ=255)                                    cAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGa                   <  1:2424333/69‑1 (MQ=255)                                    cAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGa                   <  1:2907341/69‑1 (MQ=255)                                     aGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGa                   <  1:454993/68‑1 (MQ=255)                                            cGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTgcg            >  1:1403729/1‑68 (MQ=255)                                            cGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTgcg            >  1:143804/1‑68 (MQ=255)                                            cGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTgcg            >  1:615241/1‑68 (MQ=255)                                            cGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCAAGCGATGTTgcg            >  1:3237728/1‑68 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAACCCCTACCATTTCATc                         <  1:2943535/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:3287942/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:2615607/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:1012619/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:2986425/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:3294642/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:1347166/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:965876/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:3224766/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:1830866/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:2614129/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:574415/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:670012/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:1605706/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:2193214/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:1668397/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:1977081/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:190864/48‑1 (MQ=255)                                                   ttCCGCTCGCTCCGCAGAGGTTAAAAAGCCAAACCCTACCATTTCATc                         <  1:2111717/48‑1 (MQ=37)                                                    tCCGCTCGCTCCTCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:1547236/47‑1 (MQ=37)                                                    tCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:2608393/47‑1 (MQ=255)                                                    tCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:2394335/47‑1 (MQ=255)                                                    tCCGCTCGCTCCGCAGAGGTTAAGAAGCCAAACCCTACCATTTCATc                         <  1:338268/47‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTACGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  <  1:1366367/69‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCaa        <  1:103227/63‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCaa        <  1:515632/63‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCaa        <  1:3114327/63‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCaa        <  1:2859098/63‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCaa        <  1:2273564/63‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  >  1:3257060/1‑69 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  <  1:3134590/69‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  <  1:3101716/69‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  >  1:2656978/1‑69 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  <  1:2190968/69‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  <  1:1579973/69‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  >  1:1537071/1‑69 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  <  1:881847/69‑1 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  >  1:1081444/1‑69 (MQ=255)                                                     ccgcccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAc  <  1:985277/69‑1 (MQ=255)                                                     ccgaccgcTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCaa        <  1:1620760/63‑1 (MQ=255)                                                         |                                                                 AGATGCAGGGCAAAGCGAATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGAC  >  minE/166643‑166763

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