Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,152,494	G→T	A87S (GCG→TCG)	metC →	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,152,494	0	G	T	100.0%	75.1 / NA	24	A87S (GCG→TCG)	metC	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/9);  total (15/9)

AAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCG  >  NC_000913/3152285‑3152736                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    aaaaaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAAc                                                                                                                                                                                                              >  1:46375/1‑249 (MQ=255)               ggCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGACCAGGGCGATCATGT‑Gtt                                                                                                                                                                                             >  2:138619/1‑251 (MQ=255)                gcgcGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTtg                                                                                                                                                                                            >  1:112336/1‑251 (MQ=255)                           aGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAacac                                                                                                                                                                                 <  2:132913/251‑1 (MQ=255)                                                        cTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGt                                                                                                                                                     >  1:37549/1‑250 (MQ=255)                                                                  gtAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCt                                                                                                                                          >  1:132830/1‑251 (MQ=255)                                                                            aaaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTg                                                                                                                                >  2:54649/1‑251 (MQ=255)                                                                                       acacGTAATCGCGCCAATGGAGAGTTGATCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAc                                                                                                                     >  2:121789/1‑251 (MQ=255)                                                                                            tAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGCTGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCACCACAGCCTATGCACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatcat                                                                                                                >  2:221756/1‑251 (MQ=255)                                                                                               gcgcgCAAAGGGAGAGTTGTTCAAGGGACGGGGCGGAACGTAAACCCATTCCTCCTTAAAACAAGCGATGGGGGAAGGGGAGGGGGGCGCAGGCTGCGTGTTATTTCCCTG‑GGGGTGGGCAGCGGTTGCTAATTCAATTCTTGCTTTTACCGAACGGGGCGATCATTT‑GTTGAGGACAAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGt                                                                                                             <  2:39424/250‑1 (MQ=255)                                                                                                cgcgCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAAACCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGtt                                                                                                            >  1:137579/1‑251 (MQ=255)                                                                                                 gcgcaaagggAGAGTTGTTCTATGGAGGGCGCGGAACGTTAACCCATTTCTCCTTAAAACAAGCGTTGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGTTATTCCCCGG‑GGGGTGGGCGGCGTTTGCTAATTCCATTCTTGCTTTTATCGAACGGGGCGATCATGT‑TTTGATGACCACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGttt                                                                                                           <  2:94116/243‑1 (MQ=255)                                                                                                    ccAAGGGAGAGTTTTTCTATGGACGGCGCGGAACGTTACCCCTTTTCTCCTTACAACAAGCGTTTTGTGAGCGGGAAGGGGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGGAGCGGTTGCTAATTCAATCCTTGCTTTGATCGAACAGGGCGTTCATTTTGTT‑AGGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                                                                                                        <  1:238496/251‑1 (MQ=255)                                                                                                             agTTGTTATAGGGACGGCGAGGAAAGAGAACCCATTTTTCCTTAAAACAAGAAATTGGTAAACTGGAAGGTGGGGAAGGCTGCGTTATATTTCC‑TGTCGGGTCGGAAGGGGTTGTTAATTCCATTCGTGTTTTTATCAAACAGGGCGATCATGT‑TTTGTTAACCAACTCTGCCTATGAACCGAGTTAGGATTTCTGTAGCAAAATCCTCAGCAAACTTGGCGTAACGACATCATGGTTTGATCCGCTgttt                                                                                               <  1:21162/251‑4 (MQ=255)                                                                                                                    cTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGATATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGAt                                                                                                                                                           >  1:250601/1‑184 (MQ=255)                                                                                                                         gACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATc                                                                                   >  2:216329/1‑251 (MQ=255)                                                                                                                                cgGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGc                                                                            >  1:163246/1‑251 (MQ=255)                                                                                                                                    aCGTTAACCCATTTCTCCTAAAAAAAAGGGATGTGTGAACTGGAGGGTGGCGCAGGCTGCGTGCTATTTCCCGG‑CGGGTCGGCAGCGGTTGCTAATCCCATTCTTGCTTTTATCGAACGGGGCGATAGTGT‑GTTGATGACAACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCt                                                                        <  2:51700/251‑1 (MQ=255)                                                                                                                                    aCGTTAACCCATTTCCCCTTACAACAAGCGATGTGTGAACGGGAAGGGGGCGCGGGCTGCGTGCTATTCCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAAAAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCt                                                                        <  1:54649/251‑1 (MQ=255)                                                                                                                                                         caacaaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgt                                                   >  2:73666/1‑251 (MQ=255)                                                                                                                                                                                      cgcAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCAcc                       >  2:136970/1‑250 (MQ=255)                                                                                                                                                                                        cAGGCTGCGTGCTATTTCCCTG‑CGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCACCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATg                    >  2:41657/1‑251 (MQ=255)                                                                                                                                                                                           gttgcttgttttttCCCTG‑CGGGTCGGCGGCGGTTGTTAATCCCTTTCTTGCTTTTATCGACCAGGGCGATCATGT‑GTGGATGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGaa                 <  2:96195/240‑1 (MQ=255)                                                                                                                                                                                                          ccTG‑CGGGTCGGCAGCGTTTGCTATTCCCATTCTTGCTTTTATCGAACAGGGCTATCATGT‑TTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCg  <  1:138619/251‑1 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    AAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTG‑CGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGT‑GTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCG  >  NC_000913/3152285‑3152736

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