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	100.0%	87.2 / NA	30	A319D (GCC→GAC)	msrP	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/15);  total (15/15)

CGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCAC  >  NC_000913/2040189‑2040667                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           cGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGt                                                                                                                                                                                                                                      >  1:230575/1‑251 (MQ=255)                 ggAAATATGGCTTTAAAGGGAGTAAATCGATAGCAAGTATTAAGCTGAACCGCGAGCGTTCGCCACCCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTGTGTTGATCACACGCGCTGGTCACAGGTGACCGAACGCTTTATTTGTTCAGGCGTCTTCCTCGATGTGCAGCGCCAGCCAACGCTACTGTGTAATGTTTACGCCGACCTGGTGGCATCGCTGTATCg                                                                                                                                                                                                                       <  2:344160/250‑1 (MQ=255)                    aaTATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGAGTCCGCAAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGAATCGTGGc                                                                                                                                                                                                                   >  2:143106/1‑251 (MQ=255)                     attgggCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGcc                                                                                                                                                                                                                  <  2:30332/247‑1 (MQ=255)                           cTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCCGACCCGAGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGAGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGAtt                                                                                                                                                                                                            >  1:185399/1‑251 (MQ=255)                                                   cAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTc                                                                                                                                                                                    >  1:227901/1‑251 (MQ=255)                                                          aaGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAg                                                                                                                                                                              <  1:122376/250‑1 (MQ=255)                                                                gcccGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGTCTGACGAATACGTTTTTTACGCCAACGTTAATCCGTAGGTTGATCACCCGCGCTGGTTACAGTATACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATTGTTACGCCGACCAGTTGGCATCGCTGTATCGTGGCCTGGATTTGCTGGAGAATTTTTAAATGCGTCTGACAGCAAAACa                                                                                                                                                                       <  2:239820/250‑1 (MQ=255)                                                                    gcgAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTg                                                                                                                                                                   <  2:393363/251‑1 (MQ=255)                                                                      gAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGAc                                                                                                                                                                 >  2:348144/1‑251 (MQ=255)                                                                      gAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGGTTTTACGCCAACGTGAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACTCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTGGCGGGATAAGTTCTAAATGCGTCTGACAGCAAAACAGGTGAc                                                                                                                                                                 <  1:39933/251‑1 (MQ=255)                                                                          gTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATgg                                                                                                                                                             <  1:366182/251‑1 (MQ=255)                                                                                cAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGaaa                                                                                                                                                       <  1:232164/251‑1 (MQ=255)                                                                                                       ccTGACGTATACGGTTTTTCCGCCAACGTTAATACGTATGCTGATCACCCGAGCTGGTCCCAGGCTACCGAACGATTTATTGTTTCAGGCGTCATCCTCGCTGTACAGCGCCAGCCAACGCTCCTGTTTAATGGGTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGCCTGCAAAACAGGTGACATGGCTGCACGTTTGCCTGCATCTTGCCgtctt                                                                                                                                >  2:155795/1‑247 (MQ=255)                                                                                                                 aCGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCTCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGCCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGttt                                                                                                                      <  2:230575/251‑1 (MQ=255)                                                                                                                   ggTTTGTACGCCATCGTTAAGCCGTATGTTTACCCCCCGCGCTTGTCACAGGCTCCCGAACGATTTATTTGTTCAGGTGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTTCGCCGACCAGGTGGCATCGCTTTATCGTGTCCTGGATTTGCGGGTGAATTTCTAAATGCGTCTGACAGGAAATCAGGTTACATGGCTGAAAGTTTGCCTGCAGCTTGCCGGATTGTTGCCGTTtct                                                                                                                    <  1:181653/251‑1 (MQ=255)                                                                                                                              cAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCt                                                                                                         >  1:398278/1‑251 (MQ=255)                                                                                                                              cAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCt                                                                                                         >  1:222022/1‑251 (MQ=255)                                                                                                                                    tAATCCGAATGTTGATCACCCGCGCTGGTCACAGGCTACCGACCGATATATTGGTACAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGAATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGGGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCg                                                                                                    >  1:379776/1‑250 (MQ=255)                                                                                                                                           taggTTGATCACCCGCGCTGGTCCCAGGCTACCGAACGATTTCTTGGTTCCGGCCGCAGCCTCGACGTACAGCGCCAGCCAACGCTGCTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTCGCCTGGATTTGCGGGAGAATTTCTACATGCGTCTGACAGCAAAACAGGTGACATGGCTGAACGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGCCGatcattc                                                                                             >  2:121739/4‑247 (MQ=255)                                                                                                                                           tATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGatcaatc                                                                                             >  1:425617/1‑250 (MQ=255)                                                                                                                                                       ccGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTggg                                                                                >  1:191902/1‑251 (MQ=255)                                                                                                                                                        cgcgCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGt                                                                               <  2:21886/251‑1 (MQ=255)                                                                                                                                                           gcTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGcc                                                                            <  2:35291/251‑1 (MQ=255)                                                                                                                                                           gcTGGTCACAGGCTACCGAACGATTTATTGGTGCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGGTTGCCTGCATCTTGCCGGCTTGTTGCCGTTTCTCTGGCTGGTCTGGGCGTTCAATCACGGTGGCCTGGGTGcc                                                                            >  2:214103/1‑251 (MQ=255)                                                                                                                                                                                     aTTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCTGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAGTGTGGCCGTGTCTCCGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGAGCCGGTGAAAGATATTCAGCAttt                                                  <  2:222022/251‑1 (MQ=255)                                                                                                                                                                                                    aTCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgc                                   >  1:85324/1‑251 (MQ=255)                                                                                                                                                                                                                 agctccggccaTCGCTACTGTTTAATGGTTACGCCGACCAGGGTGCATCTCTGGATCGTGGCCTGGATTTGCGGGATAATTTCTAAATGCGTCTGACAGCATAACAGGTGACATTGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTTCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGGGGACTGGGTGCCGAGCCGGTGAAAGATAGTCAGCCTTTTACTGGTCGCCCTGCACTGAAATTTTTg                      <  2:316307/244‑1 (MQ=255)                                                                                                                                                                                                                          caACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAcc             <  1:348144/251‑1 (MQ=255)                                                                                                                                                                                                                                     tttAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGACTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCAc  >  1:53657/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           CGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCAC  >  NC_000913/2040189‑2040667

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