Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,442,061	A→G	V39A (GTA→GCA)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,442,061	0	A	G	100.0%	58.0 / NA	17	V39A (GTA→GCA)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/9);  total (8/9)

AAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAA  >  NC_000913/4441930‑4442193                                                                                                                                    |                                                                                                                                     aaGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGt                                                                                                                               <  1:86532/139‑1 (MQ=255)     cGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCAtt                                                                                                                           >  2:366523/1‑139 (MQ=255)                         cgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATg                                                                                                       <  2:351380/139‑1 (MQ=255)                                 ccTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCTAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTgg                                                                                               <  2:266406/139‑1 (MQ=255)                                   tGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCt                                                                                             >  2:104705/1‑139 (MQ=255)                                     cggcggTGCTGTAAGCGCCGGGTAACTGCCAGAACAGACGATCCACACCCCAGAAACAACCCATCGCAAAAACGGCAATCTCCATTCCGTCAGGTGCATAGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTAc                                                                                           >  1:277638/1‑139 (MQ=255)                                        cggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACggg                                                                                        <  1:104705/139‑1 (MQ=255)                                                             aCTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcagg                                                                   <  2:352198/139‑1 (MQ=255)                                                                          aCGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCAt                                                                                     <  1:324021/108‑1 (MQ=255)                                                                          aCGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCAt                                                                                     >  2:324021/1‑108 (MQ=255)                                                                                               acCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGt                                                                                                                       <  2:48224/53‑1 (MQ=255)                                                                                               acCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGt                                                                                                                       >  1:48224/1‑53 (MQ=255)                                                                                               acCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTa                                 >  2:5754/1‑139 (MQ=255)                                                                                               acCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTa                                 >  2:332999/1‑139 (MQ=255)                                                                                                cccATCGCTAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTAt                                <  2:277638/139‑1 (MQ=255)                                                                                                               gCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATgg                 <  1:332999/139‑1 (MQ=255)                                                                                                                              tCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCACGaa  >  1:108610/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                     AAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAA  >  NC_000913/4441930‑4442193

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