Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,389,934	A→G	Y143H (TAC→CAC)	psd ←	phosphatidylserine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,389,934	0	A	G	100.0%	65.6 / NA	19	Y143H (TAC→CAC)	psd	phosphatidylserine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/10);  total (9/10)

AGGCAAATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTAGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCA  >  NC_000913/4389798‑4390047                                                                                                                                         |                                                                                                                  aGGCAAATCACGCGTTTGTTACGGGCAAACAGATTCGGAACGTTCTGAGACGTAAGATGGTTAACGGAGAAGAGATCGCCGGGCACGTAGATGATCTCACGCAGAATACCGTTGCACGGCATGTGCACGCggtggtggt                                                                                                                 <  2:104550/139‑1 (MQ=255) aGGCAAATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCggtggtggt                                                                                                                 <  1:593671/139‑1 (MQ=255)      aaTCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAcgc                                                                                                            >  2:293132/1‑139 (MQ=255)                                                aGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAgg                                                                  >  1:799488/1‑139 (MQ=255)                                                  ccGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                <  1:496179/139‑1 (MQ=255)                                                  ccGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                <  2:54002/139‑1 (MQ=255)                                                     tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTg                                                                                             >  1:345077/1‑107 (MQ=255)                                                     tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTg                                                                                             <  2:345077/107‑1 (MQ=255)                                                       agaTGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTcc                                                               <  1:460478/135‑1 (MQ=255)                                                       agaTGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTcc                                                               >  2:460478/1‑135 (MQ=255)                                                                                     aCGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGa                             <  2:799488/139‑1 (MQ=255)                                                                                          gATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTg                        >  1:335062/1‑139 (MQ=255)                                                                                                 tcACGCAGAATACCGTTNCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTtgt                 <  1:27048/139‑1 (MQ=255)                                                                                                 tcACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTtgt                 >  2:439508/1‑139 (MQ=255)                                                                                                         aaTACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTg         <  1:293132/139‑1 (MQ=255)                                                                                                               gTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGc   <  1:439508/139‑1 (MQ=255)                                                                                                                ttGCACGGCATTTTTCCGCTATGTTCGTCACTCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCCCCATCCGATAGTTCCCTGCCAGCCGTCCTTCGATGCTCTAGTTGTCGCCTTTGGCTTGCa  <  2:560014/139‑1 (MQ=255)                                                                                                                   cACGGCATGTGTACGCGGTGGTGGTCACGCggggc                                                                                                       >  2:764143/1‑34 (MQ=38)                                                                                                                   cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATc                                                         <  1:638154/81‑1 (MQ=255)                                                                                                                   cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATc                                                         >  2:638154/1‑81 (MQ=255)                                                                                                                     cGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTg                                                                                             <  1:834305/43‑1 (MQ=255)                                                                                                                     cGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTg                                                                                             >  2:834305/1‑43 (MQ=255)                                                                                                                                         |                                                                                                                  AGGCAAATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTAGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCA  >  NC_000913/4389798‑4390047

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