Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	39,361	C→T	34.1%	P83P (CCG→CCA)	caiT ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	39,361	0	C	T	34.1%	55.2 / 38.0	73	P83P (CCG→CCA)	caiT	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (37/11);  new base T (14/11);  total (51/22)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ACAGGAGGCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGC  >  minE/39313‑39424                                                 |                                                                aCAGGAGGCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCtt                                             >  1:2762852/1‑69 (MQ=255) aCAGGAGGCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCtt                                             >  1:957650/1‑69 (MQ=255) aCAGGAGGCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCtt                                             >  1:734460/1‑69 (MQ=255)  cAGGAGGCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCttt                                            >  1:1866663/1‑69 (MQ=255)  cAGGAGGCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCttt                                            >  1:892732/1‑69 (MQ=255)       ggCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        <  1:2110901/68‑1 (MQ=255)       ggCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        <  1:785164/68‑1 (MQ=255)        gCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCa                                      <  1:184062/69‑1 (MQ=255)        gCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCa                                      <  1:2677297/69‑1 (MQ=255)           aaCATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATAc                                   >  1:1911695/1‑69 (MQ=255)           aaCATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATAc                                   >  1:1425001/1‑69 (MQ=255)           aaCATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATAc                                   >  1:902047/1‑69 (MQ=255)           aaCATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATAc                                   >  1:2284134/1‑69 (MQ=255)           aaCATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATAc                                   >  1:1149395/1‑69 (MQ=255)           aaCATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATAc                                   >  1:2578466/1‑69 (MQ=255)             catcatAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACgg                                 >  1:2661932/1‑69 (MQ=255)             catcatAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACgg                                 >  1:2425390/1‑69 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTaaa                                                  >  1:3337023/1‑48 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        <  1:3196263/58‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        <  1:3263487/58‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        >  1:1753872/1‑58 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        >  1:1719200/1‑58 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        <  1:2245717/58‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        >  1:676013/1‑58 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        >  1:1101662/1‑58 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        >  1:1122708/1‑58 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             >  1:283759/1‑69 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             >  1:2748847/1‑69 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:971940/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             >  1:2040692/1‑69 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             >  1:1874099/1‑69 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             >  1:1974000/1‑69 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:22487/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGATTTTTGGAATACg                                  >  1:297998/1‑64 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:2273620/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:2490633/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:3302542/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:356383/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:139610/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:1334208/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:1360950/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:1381338/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:267958/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:2249890/69‑1 (MQ=255)                 atAAAGATCCAACTGGCGGTACTAAATTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:789281/69‑1 (MQ=255)                 atAAAGATCCAACTGCCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCa                             <  1:626431/69‑1 (MQ=255)                 atAAAGATCCAACTAGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTgg                                        >  1:29383/1‑58 (MQ=255)                 atAAAGATACAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCttt                                            >  1:576892/1‑54 (MQ=255)                         ccAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACAcca                      >  1:2888254/1‑68 (MQ=255)                             cTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACAccagccag                 <  1:2706154/69‑1 (MQ=255)                                     gCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaac          >  1:740061/1‑68 (MQ=255)                                     gCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaac          >  1:3227859/1‑68 (MQ=255)                                     gCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaac          >  1:1432452/1‑68 (MQ=255)                                     gCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaac          >  1:3027751/1‑68 (MQ=255)                                      cTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaacc         >  1:816987/1‑68 (MQ=255)                                      cTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaacc         >  1:2311138/1‑68 (MQ=255)                                      cTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaacc         >  1:3164186/1‑68 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCaaa                           >  1:1150281/1‑48 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccagc          >  1:850589/1‑63 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:653498/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:3055009/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:1879133/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:2756526/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:2704379/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:3218113/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:1589524/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:835748/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:2531566/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:919731/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:933840/1‑67 (MQ=255)                                        aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg        >  1:1212148/1‑67 (MQ=255)                                        aaaTTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACAcc                       >  1:2494825/1‑52 (MQ=255)                                             cTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGc  >  1:2380840/1‑68 (MQ=255)                                                 |                                                                ACAGGAGGCGAACATCATAAAGATCCAACTGGCGGTGCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGC  >  minE/39313‑39424

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