Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	467,736	G→A	G109S (GGT→AGT)	cof →	HMP‑PP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	467,736	0	G	A	94.7%	39.4 / NA	19	G109S (GGT→AGT)	cof	HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/6);  minor base C (1/0);  total (13/6)
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.49e-01

CGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACG‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTTCCGCCACGGATTGCCTCGA  >  NC_000913/467505‑467953                                                                                                                                                                                                                                        |                                                                                                                                                                                                                           cGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCTTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTTCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAagn                                                                                                                                                                                                        <  1:219573/251‑2 (MQ=255)                                      ntCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAagag                                                                                                                                                                                                       >  1:65531/2‑215 (MQ=255)                                                      ncagcaTATTCTCGGGGCGCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCa                                                                                                                                                        >  1:66606/2‑246 (MQ=255)                                                                  cGGGGCGATATCGATGGATGCGTATTTGAATACCGGCATAGGAACGCGGGTGCATTCTCTGGAGGGTCGACTTTTACACCGTGATGATTTACCAGCGGATGTCGCGAAGCTGGTGCTGGATCAGCAATGGGAGATCCGAGCCAGCATGCATATCTTCAAAGAC‑ACAAGTTGGTTTAACGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCAGTTTTCGTTATCAGATAATCGATGTCaaa                                                                                                                                        >  1:261285/1‑250 (MQ=255)                                                                             cGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAgn                                                                                                                                                                      <  1:210958/209‑2 (MQ=255)                                                                               cTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCgg                                                                                                                          >  1:265437/1‑251 (MQ=255)                                                                                  gATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGAGGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATc                                                                                                                                                         <  1:289192/217‑1 (MQ=255)                                                                                                       aaCGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATTATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAgn                                                                                                                                                                                    <  1:88173/169‑2 (MQ=255)                                                                                                           naaCGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGAGGCTGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATACTCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc                                                                                              >  1:65659/2‑251 (MQ=255)                                                                                                                      cATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAg                                                                                                                                                                       >  1:94475/1‑167 (MQ=255)                                                                                                                                                  gTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCGTTTGTNTATNNNNNtt                                                                                                                                                                  >  2:214754/1‑144 (MQ=255)                                                                                                                                                             ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAgg                                            >  1:125507/1‑251 (MQ=255)                                                                                                                                                             ccTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCCTTATCAGCTAATCGATGTCACAAAAATg                                                                                                                                  >  1:236294/1‑165 (MQ=255)                                                                                                                                                                   gATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATCTCTTCAATGACGACA‑GTTGGCTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGCTCTGCTTCTGTGGCGATCCCGCCGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                                      >  1:247135/1‑251 (MQ=255)                                                                                                                                                                                                        cccGAGCCAGCATGCATATCTTCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCCGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTTCCGCCACGGATTGCCTCGa  >  1:60598/1‑250 (MQ=255)                                                                                                                                                                                                         ccGAGCCAGCATGCCTATCTTCAATGACGACC‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGTAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTTCCGCCACGGATTGCCTCGa  >  1:93819/1‑249 (MQ=255)                                                                                                                                                                                                                       atatCTNNAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCCTTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATg                                                                                                                                             <  2:23715/96‑1 (MQ=255)                                                                                                                                                                                                                             tCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAg                                                                                                                                                                                     <  2:1453/50‑1 (MQ=255)                                                                                                                                                                                                                             tCAATGACGACA‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAg                                                                                                                                                                                     >  1:1453/1‑50 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                           CGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCGTATTTGATTACCGGCAACGGAACGCGCGTGCATTCTCTGGAAGGTGAACTTTTACATCGTGATGATTTACCTGCGGATGTCGCGGAGCTGGTGCTGTATCAGCAATGGGATACCCGAGCCAGCATGCATATCTTCAATGACGACG‑GTTGGTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTTCCGCCACGGATTGCCTCGA  >  NC_000913/467505‑467953

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