Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	397,141	C→T	12.1%	K216K (AAG→AAA)	lipA ←	lipoate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	397,141	0	C	T	12.1%	88.4 / 4.5	33	K216K (AAG→AAA)	lipA	lipoate synthase
Reads supporting (aligned to +/- strand):  ref base C (13/16);  new base T (2/2);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GTCCCACCATCAGACCAGACTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTAAATACGCGGTACGTTTTCCAG  >  minE/397077‑397206                                                                 |                                                                  gTCCCACCATCAGACCAGACTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGAGTTCCAGCAGCTTCAGc                                                             <  1:2091732/71‑1 (MQ=255)            agaccagacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGTTTCAGCGACCAGTTGTa                                                  <  1:2157565/71‑1 (MQ=255)            agaccagacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGTTTCAGCGACCAGTTGTa                                                  <  1:813191/71‑1 (MQ=255)            agaccagacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGt                                                   <  1:1025291/70‑1 (MQ=255)            agaccagacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGt                                                   <  1:1963422/70‑1 (MQ=255)            agaccagacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGt                                                   <  1:1565136/70‑1 (MQ=255)            agaccagacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGt                                                   >  1:767438/1‑70 (MQ=255)            agaccagacTTGGTCGGGATTTCCGGATGCGCTCCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGt                                                   >  1:1194652/1‑70 (MQ=255)                cagacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCt                                              >  1:6263/1‑71 (MQ=255)                 agacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTg                                             <  1:250443/71‑1 (MQ=255)                 agacTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTg                                             <  1:1142023/71‑1 (MQ=255)                       ggTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCt                                                                  >  1:796307/1‑44 (MQ=255)                       ggTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACAAGTTGTAAt                                                >  1:2169362/1‑62 (MQ=255)                        gTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAg                                       <  1:929615/70‑1 (MQ=255)                            ggATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCg                                   >  1:71931/1‑70 (MQ=255)                            ggATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGt                                  <  1:1751980/71‑1 (MQ=255)                                  ccGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTa                                 >  1:539650/1‑66 (MQ=255)                                     gATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAg                                       <  1:527195/57‑1 (MQ=255)                                      aTGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACgg                         >  1:2959/1‑70 (MQ=255)                                       tGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTa                       <  1:1954308/71‑1 (MQ=255)                                       tGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTa                       <  1:101450/71‑1 (MQ=255)                                         cgcTTCTTTAAAGCGTTCCAGCAGTTTCAGCGACCAGTTGTAATCTGCAc                                          >  1:492116/1‑50 (MQ=255)                                         cgcTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTaaa                     <  1:982772/71‑1 (MQ=255)                                            ttcttTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGAc                           <  1:921465/62‑1 (MQ=255)                                            ttcttTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTaaa                     <  1:1730274/68‑1 (MQ=255)                                             tcttTAAAGCGTTCCAGCAGTTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCt                              >  1:1070403/1‑58 (MQ=255)                                              cttTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTa                                                  <  1:1931994/37‑1 (MQ=255)                                              cttTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGa                            >  1:155746/1‑59 (MQ=255)                                              cttTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTAAATAcgc                >  1:708877/1‑71 (MQ=255)                                              cttTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTAAATAcgc                >  1:887153/1‑71 (MQ=255)                                              cttTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTAAATAcgc                >  1:901581/1‑71 (MQ=255)                                                ttAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTAAATACGCgg              >  1:502366/1‑71 (MQ=255)                                                 tAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTAAATACGCgg              <  1:2051311/70‑1 (MQ=255)                                                            agcagcTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTAAATACGCGGTACGTTTTCCAg  >  1:285407/1‑71 (MQ=255)                                                              cagcTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTa                       <  1:436333/48‑1 (MQ=255)                                                               agcTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTAc                                >  1:1759482/1‑38 (MQ=255)                                                                 |                                                                  GTCCCACCATCAGACCAGACTTGGTCGGGATTTCCGGATGCGCTTCTTTAAAGCGTTCCAGCAGCTTCAGCGACCAGTTGTAATCTGCACCAGGCCGTACCTGACGGTAAATACGCGGTACGTTTTCCAG  >  minE/397077‑397206

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