Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,040,289	C→T	A271V (GCG→GTG)	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,289	0	C	T	100.0%	82.8 / NA	26	A271V (GCG→GTG)	msrP	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/13);  total (13/13)

ACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGG‑TGTTTATGGCAAGGCGTT‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTT  >  NC_000913/2040055‑2040533                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     cccgttttATCGGCGGGGGGTTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTAGG‑TGTTTAGGGCAAGGCTTT‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATAc                                                                                                                                                                                                                                         <  1:146838/243‑1 (MQ=255)        tATCGGCGGGGGGCGAAAAAATCCTTAGGTCAAAGGATGGTGCCTCGACGAAGCAAGGCCTCCCCTCACACGAATGACCGTGGG‑TGTTTATGGCAAGG‑GTTGACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTAc                                                                                                                                                                                                                                <  1:83864/250‑1 (MQ=255)                            atCCTTATCTAGAAGGATTCCGTCTCGACGAAGCACTGAACCCCCTAACACTGATGACCGTAGC‑AGTTTATGGCAAGGAGTT‑ACCGACACAACATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGACTTTAAAGGGATTAAATCGATAGTCCGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGCATACGGTTTTTACGCCACCGTTAATCCGTATGtt                                                                                                                                                                                                           >  2:366630/1‑251 (MQ=255)                                                       acgaAGCAATGCATCCCCTCACACTGATG‑CCGGAGGGTGTTAAGGGCAGGGCGTT‑CCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTAcc                                                                                                                                                                                <  1:207235/251‑1 (MQ=255)                                                               agccaccccctcaagcgaaagagcGTAGG‑TTTTTAGGGCAGGGTGTT‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCGGGAATCGGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGGTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGAtt                                                                                                                                                                        <  1:42604/228‑1 (MQ=255)                                                                                    aCCGTAGG‑TGTTTATGGCAAGGCGTT‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATcc                                                                                                                                                                                                                  >  2:356267/1‑188 (MQ=255)                                                                                             tgtTTATGGCAAGGCGTT‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTAc                                                                                                                                           >  2:425170/1‑251 (MQ=255)                                                                                                            gTT‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGCTTGTGCCGTGGAAATATGGCGTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCACCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGCACGTTTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGc                                                                                                                            >  2:51201/1‑251 (MQ=255)                                                                                                            gTT‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTGGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGc                                                                                                                            >  2:170109/1‑251 (MQ=255)                                                                                                             tt‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCt                                                                                                                           >  1:141973/1‑251 (MQ=255)                                                                                                                       aaaaaTGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTaa                                                                                                                  <  2:198833/250‑1 (MQ=255)                                                                                                                        aaaaTGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAAt                                                                                                                 >  2:159875/1‑251 (MQ=255)                                                                                                                                  gcCGGTGGGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTAcgccg                                                                                                       <  2:30543/251‑1 (MQ=255)                                                                                                                                   ccGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTAcgccga                                                                                                      <  2:141973/251‑2 (MQ=255)                                                                                                                                     ggTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGCACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTCATGGTTAcgccgacc                                                                                                    >  2:568480/1‑248 (MQ=255)                                                                                                                                             tGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGcc                                                                                            >  1:129210/1‑250 (MQ=255)                                                                                                                                                  gtgCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGTCATCGCt                                                                                       <  1:138803/251‑1 (MQ=255)                                                                                                                                                  gggCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGAGCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCt                                                                                       <  2:190756/249‑1 (MQ=255)                                                                                                                                                                tGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTgg                                                                         >  1:79892/1‑251 (MQ=255)                                                                                                                                                                       aaaGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGTCAGTGCTTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCAGCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGc                                                                   <  1:49532/250‑1 (MQ=255)                                                                                                                                                                       aaaGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGc                                                                   <  1:319903/250‑1 (MQ=255)                                                                                                                                                                                      aTAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAt                                                   >  1:259729/1‑251 (MQ=255)                                                                                                                                                                                                              gAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGAc                           >  1:29405/1‑251 (MQ=255)                                                                                                                                                                                                                             accTGGATTCTTGCAGTGCCTGTCGAATCCGGTATTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGTTCACAGGCTACCGAACTATTTATTGGTTGATGCGGCATCCTCGATGTTCAGCGCCCGCCAACGCTACTGTGTAATGGTTGCGACGTCCAGGTGGCATCTCTGTATCGGGGCCTGGATTTGCGGGACAATTTCTAAATTCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTg            <  1:48058/251‑1 (MQ=255)                                                                                                                                                                                                                                tGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCt         >  2:67987/1‑251 (MQ=255)                                                                                                                                                                                                                                       gggCAGTGCCTGACGAATACGGTTTGTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTAGTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGGTTAATGGTTACGCCTACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTGTGCCTGCATCtt  <  2:558152/250‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     ACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGG‑TGTTTATGGCAAGGCGTT‑ACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTT  >  NC_000913/2040055‑2040533

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