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	96.2%	58.0 / NA	26	A87S (GCG→TCG)	metC	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (15/10);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATC  >  NC_000913/3152247‑3152699                                                                                                                                                                                                                                                        |                                                                                                                                                                                                              aaaGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGgtcgg                                                                                                                                                                                                            >  1:151100/1‑251 (MQ=255) aaaGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGAAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGgtcgg                                                                                                                                                                                                            >  2:206188/1‑251 (MQ=255)     cTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGc                                                                                                                                                                                                        >  2:210559/1‑251 (MQ=255)           aCTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCTCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGGGTGAACTGGAAGGGGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGTTTg                                                                                                                                                                                                   <  2:258138/250‑1 (MQ=255)           aCTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATATCCCTGCGGGTCGGCAGCGGTTGc                                                                                                                                                                                                  >  1:70940/1‑251 (MQ=255)                                cGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGAGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAACAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTAt                                                                                                                                                                             >  1:90545/1‑251 (MQ=255)                                cGCAACAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCTCTTCTTCGCTGTTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGGAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATTTGTGAACTGGAAGGTGGCGCAGGCTTCGTGCTATTTCCCTGCGGGTCGGCAGCGTTTGCTAATACCATTCTTGCGTTTAt                                                                                                                                                                             <  2:242420/251‑1 (MQ=255)                                                       cgcgGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTATATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGGGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtg                                                                                                                                                       >  1:91635/1‑250 (MQ=255)                                                              aaTAGCGGGATGCAGCTCGCGTCTTCGCTGGGCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTGAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGATGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTGTTATCGAACAGGGCGATCATGTGTTGATGACCaa                                                                                                                                               <  2:222900/251‑1 (MQ=255)                                                                     tGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGcc                                                                                                                                        >  1:78898/1‑251 (MQ=255)                                                                           agggcgcTTCTTCGCTGGTCTTTGAAAGTGTAGAAGACAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGa                                                                                                                                   >  2:190430/4‑250 (MQ=255)                                                                                gcTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCg                                                                                                                               >  1:36828/1‑249 (MQ=255)                                                                                     ttcGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGCTGACCAACACCGCCTATGAACCGAGTCAgg                                                                                                                        >  2:101399/1‑251 (MQ=255)                                                                                           ggTCTTGGACAGTGTAGAAGCAAAAAAACCCGCGACCCGTAATCCCGCCAATGGAGATTTGTTCTAGGGACGGCGCGGAACGTTACCCCATTCCTCCTTACAACAAGCGATGTGTGAACTGGAAGGGGGCGCGGGCGGCGTGTTATTCCCCGGCGGGTCGGCAGCGGTGGCTAATCCCATTCTTGCTTTTATCGACCGGGGCGATCTTGTGTTGAGGACCAACACCGCCTAGGAACCGAGTCAGGATTTc                                                                                                                   <  2:50956/250‑1 (MQ=255)                                                                                                     aGTGTAGAAGCCAAAAACCACGCGACACGTCATCGCGCCAATGGGGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCCTTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCACAAt                                                                                                        >  1:21035/1‑251 (MQ=255)                                                                                                           gAAGCCAAAAACAACGGGAAACGAAATAGAGCCAATGGAGGGTTGTTCTATGGCCGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGAGAAGGATGCGTGCTATTTACCTGCGGGTCGGCAGCGGTTGCTAATTCAATTCTTGCTTTTATCGAAAAGGGAGATCATGTGTTGATGACCACCACCGCATATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                                                                    >  2:148918/1‑249 (MQ=255)                                                                                                                           cgACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGACCTGGAAGGTGGCGCCGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACg                                                                                  >  2:22108/1‑251 (MQ=255)                                                                                                                                  tAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACACGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGAAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatcat                                                                           >  1:249315/1‑251 (MQ=255)                                                                                                                                                gagaGTTTTTCTAGGGCCGGCGCGGACCGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGGGGCGCAGGCTGCGTGCTATTTCCCGGGGGGTGGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACGGGGCGATCATGTGTTGAGGACAAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAGCCGCTg                                                             <  1:47407/251‑1 (MQ=255)                                                                                                                                                    gttgttCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTACATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTg                                                         >  1:262057/1‑251 (MQ=255)                                                                                                                                                       gttCTATGGAGGGCCCGGAACTTGAACCCTTTTCTCCTTACAACAAGCGAGGTGTGAACGGGAGGGTGGCGCAGGCTGCGTGCTATTTCCCGGGGGGTGGGCAGCGTTTGCTAATCCCATCTTTTCTTTTACCGAACAGGCCGATCATTTTTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCGGATTGGt                                                       <  2:236781/250‑1 (MQ=255)                                                                                                                                                            agggagggCGCGGAACGTTAACCCTTTCCCCCTTCAAACAGGCGAGGTGGAAACGGAAAGGGGGGGCAGGCTGCGTTTTTTTTCCCTGGGGGGGGGCAGCGGTTGCTAATTCCATTCTGGCTTTTATCAAACAGGGCGATCAGGTGTTGATGACAAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa                                                  <  2:36828/244‑1 (MQ=255)                                                                                                                                                            aTGGAGGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGGGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGTTTTTATCGAACAGGGCGATCATTTGTTAATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGat                                                 <  1:26366/251‑1 (MQ=255)                                                                                                                                                            aTGGACGGCGCGGAACGTGAACCCATGTCTCCTTACAACAAGCGAGGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCTCCTATGAACCGAGTCAGGATTGCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTGGATCCGCTGATTGGTGACGat                                                 <  2:84290/251‑1 (MQ=255)                                                                                                                                                                                                    aGCGATGTGTGAACTGGACGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTc         >  1:81518/1‑251 (MQ=255)                                                                                                                                                                                                     gCGATGTGTGCAGGGGAGGGGGGCGCGGGCTTGGTGCTATGTCCGTGGGGGTCGTCAGCGTGTGTTAATTCCATTCTGGTTTTTATCGAAGGGGGCGATCTTGTTTTGATGACAAACTCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCTATATCGTTAAGCATCTGCAGCCAAACAATAAAATCGTGTTTCt        <  2:50108/251‑1 (MQ=255)                                                                                                                                                                                                           tgtgaaataggaagtagcgcCGGCTGCGTGCTATTTCCCTGCGGGTAGGCAGCGGTTGCTAACTCCATTCATGCCTTTATAGAACCGGGCGATCATGTGTTGATGAACAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGGTTGGTGCCGAGATCGTTACGCATCTGCAGCCAAACACTAATATCGTGTTTCTGGactc  >  2:100324/17‑248 (MQ=255)                                                                                                                                                                                                            gtgaatgggagggTGGCGGGGGCTGCTTGCTTTTTCCCTGGGGGTGGGCTGCGGTTGTTAATACCATTCTTCTTTTTATCGAACGGGGCGATCATGTGTGGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATc  <  2:212802/239‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                              AAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATC  >  NC_000913/3152247‑3152699

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