Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,040,433	C→A	A319D (GCC→GAC)	msrP →	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,433	0	C	A	93.8%	41.1 / ‑6.6	16	A319D (GCC→GAC)	msrP	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (8/7);  minor base T (0/1);  total (8/8)
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.92e-01

AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCG  >  NC_000913/2040207‑2040653                                                                                                                                                                                                                                   |                                                                                                                                                                                                                             aaaTATGGCTTTAAAGGGATTAAATCGATAGGCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTGTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGGTCAGGCGGCATCCTCGATGTATAGCGCCATCCAACGCTACTGGTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTgg                                                                                                                                                                                                      <  1:336802/251‑1 (MQ=255)                                         ngCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGGGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGACTTTCTAAATGCGTCTGACAGc                                                                                                                                                               >  1:308884/2‑250 (MQ=255)                                                                  ccccTGGAATCTGGCAGCGCCTGACGAAGACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCATCCAACGCTACTGTTTAATGGGTACGCCTTCCAGGGGGCAGCGCTGTATCGGGGCCTTGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGc                                                                                                                                              <  1:260911/240‑1 (MQ=255)                                                                                            aaTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGGTTTATTGGTTCATGCGGCAGCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCCGGATTTGCGGGAGAATTTCTAAATGCGGCTGACAGAAAAACAGGTGACATGGc                                                                                                                                              <  1:487228/216‑1 (MQ=255)                                                                                                           ccAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCt                                                                                                                                                                     >  1:391898/1‑178 (MQ=255)                                                                                                                            gTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGACGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGATTAATGGTTACGCCGACCAGGAGGCATAGCTGTATTGTGGCCTGGATTTGCGGGAGAATTTGTAAATGc                                                                                                                                                                         >  1:292949/1‑157 (MQ=255)                                                                                                                                            ggTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTAGCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGTATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGa                                                            <  1:279440/250‑1 (MQ=255)                                                                                                                                                         cGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTc                                                                                                                                                                                >  1:428332/1‑121 (MQ=255)                                                                                                                                                           aaCGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTg                                                                                                                                                                                                              >  1:187883/1‑89 (MQ=255)                                                                                                                                                                       gTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGc                                                                                                             >  1:284715/1‑174 (MQ=255)                                                                                                                                                                            ggcggcATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATAGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCg                           >  1:423522/1‑251 (MQ=255)                                                                                                                                                                                           gTNCAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGc                                                                                                                               <  2:308884/136‑1 (MQ=255)                                                                                                                                                                                                  gccagccaACGCTACTGTTTAATTGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGTATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGATGACATGGCTGAATGTTTGCCTGCATc                                                                                                                            <  1:402208/132‑1 (MQ=255)                                                                                                                                                                                                      gccaACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAACGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCg  >  1:457718/1‑250 (MQ=255)                                                                                                                                                                                                            gCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAGGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGGTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCgg                                                       <  1:281498/191‑1 (MQ=255)                                                                                                                                                                                                                      aaTGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAGGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCg                                                                                                                      <  1:14701/118‑1 (MQ=255)                                                                                                                                                                                                                              cgccgaCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCGGAAAGTTTGTCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGcc                                                              <  1:284039/166‑1 (MQ=255)                                                                                                                                                                                                                               gccgacCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACa                                                                                                                                                                 >  1:117095/1‑66 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                             AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCG  >  NC_000913/2040207‑2040653

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