Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	672,388	G→T	28.6%	S799Y (TCC→TAC)	leuS ←	leucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	672,388	0	G	T	28.6%	32.8 / 13.3	21	S799Y (TCC→TAC)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (4/11);  new base T (3/3);  total (7/14)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGT  >  NC_000913/672240‑672525                                                                                                                                                     |                                                                                                                                          aaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGa                                                                                                                                          <  2:216060/150‑1 (MQ=255) aaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGa                                                                                                                                          <  2:216050/150‑1 (MQ=255)                     cGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAAACTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCttttt                                                                                                                      >  2:65019/1‑150 (MQ=255)                                          ccAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGgcg                                                                                                 >  1:114458/1‑150 (MQ=255)                                           cAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGgcgc                                                                                                <  1:46733/150‑1 (MQ=255)                                                                        ccTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCGTTCAGtt                                                                   >  1:109242/1‑150 (MQ=255)                                                                        ccTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTGGTCGATATCGCCTTCGCCTTTCAGtt                                                                   >  1:105340/1‑150 (MQ=255)                                                                               ttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGTAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCca                                                            >  2:227330/1‑150 (MQ=255)                                                                               ttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGTAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCca                                                            >  2:227337/1‑150 (MQ=255)                                                                               ttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGTAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCca                                                            >  2:47395/1‑150 (MQ=255)                                                                                                          tttGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGaa                                 <  1:71732/150‑1 (MQ=255)                                                                                                            tGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACg                               <  2:35831/150‑1 (MQ=255)                                                                                                               cACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGt                            <  1:54895/150‑1 (MQ=255)                                                                                                                  ggACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTaa                         <  1:172410/150‑1 (MQ=255)                                                                                                                         aCCGTTAACCCGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACg                  <  2:209127/150‑1 (MQ=255)                                                                                                                          ccGTTAACCTGCACCACGACCAGCGTGTAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACgg                 <  1:227330/150‑1 (MQ=255)                                                                                                                          ccGTTAACCTGCACCACGACCAGCGTGTAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACgg                 <  1:227337/150‑1 (MQ=255)                                                                                                                          ccGTTAACCTGCACCACGACCAGCGTGTAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACgg                 <  1:47395/150‑1 (MQ=255)                                                                                                                           cGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGa                <  1:226570/150‑1 (MQ=255)                                                                                                                                         acGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGt  <  2:101782/150‑1 (MQ=255)                                                                                                                                         acGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGt  <  1:134184/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                          AATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGT  >  NC_000913/672240‑672525

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑