Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,390,331	C→T	Q10Q (CAG→CAA)	psd ←	phosphatidylserine decarboxylase proenzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,390,331	0	C	T	100.0%	81.6 / NA	27	Q10Q (CAG→CAA)	psd	phosphatidylserine decarboxylase proenzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/11);  total (16/11)

CGCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAA‑AGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTCCCGGAT  >  NC_000913/4390087‑4390576                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cGCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGc                                                                                                                                                                                                                                                  <  2:78235/251‑1 (MQ=255)  gCCATCGGCAGGCATGACCAGTACATTCGGTCCGGTATCGATTGGGCGTTCGTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCAGGTCGTCCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGGCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCg                                                                                                                                                                                                                                                 <  2:186200/251‑1 (MQ=255)                  aCCAGTACATTCGGAGCGTTATCGATTGGGCGTACTTCGTCTCGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAGCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCAAAATTTTAAATGAAtt                                                                                                                                                                                                                                 <  2:23435/251‑1 (MQ=255)                    cAGTACATTCGGATCGGTATCGCGTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGAGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTGAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCAGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAAttt                                                                                                                                                                                                                                >  2:77343/1‑250 (MQ=255)                     aGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTaa                                                                                                                                                                                                                              <  1:34371/251‑1 (MQ=255)                      gTACATTCGGAGCGTTATCGATTGGGCGTACTTCTGCACGCAGCGGACGTACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTGGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTGGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATTGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAAc                                                                                                                                                                                                                             <  2:269107/251‑1 (MQ=255)                        aCATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACaa                                                                                                                                                                                                                           <  1:18436/251‑1 (MQ=255)                          aTTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAgg                                                                                                                                                                                                                         >  1:117912/1‑251 (MQ=255)                                                 gTACTTCGTCGCGCAGCGGGCGGACACAGAATTCGTTAAAGGCGCGGTAACTGGCGGTGTCCGGCTTTTGCGCCTCTTTAATGTCGACCTTGTAGTATTTAACGAAACGATCGATAACCAGTTTTGTCAGCCATCCTACCCGCGTGATTGCGCCCCACCCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGCATGTATTGTAGCGAAAGTTTAACTGAATTTAACAAGGTGGCCTCCAGGCCATTGTTTTGt                                                                                                                                                                                                  >  2:116819/1‑251 (MQ=255)                                                              cAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCg                                                                                                                                                                                     <  1:99291/251‑1 (MQ=255)                                                                      ggacAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACAt                                                                                                                                                                             >  2:258094/1‑251 (MQ=255)                                                                                     tAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTATTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGaaa                                                                                                                                                               <  2:169472/250‑1 (MQ=255)                                                                                       aaGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCg                                                                                                                                                                                                                                                 <  1:169515/166‑1 (MQ=255)                                                                                       aaGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGaaaaa                                                                                                                                                             >  2:189924/1‑250 (MQ=255)                                                                                       aaGGTGCGGTAGCTGGCGGAGTCCGGCTCTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCg                                                                                                                                                                                                                                                 >  2:169515/1‑166 (MQ=255)                                                                                                               ggCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTa                                                                                                                                    <  1:122100/251‑1 (MQ=255)                                                                                                                           tctTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtc                                                                                                                         >  1:7710/1‑250 (MQ=255)                                                                                                                                      gACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTAt                                                                                                             >  2:77412/1‑251 (MQ=255)                                                                                                                                       aCCTTGTAGTATTTAACGAACAGCTCGATAACCCGGTTTGTCCGCCATCCTGCCCGCTTGCTTGCGCCCCAACACGCCCGGCGAGTAAGCCATAGTTTCGGCCGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCACGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATCCCGGATCCTGAACAACAGGGGACGATTCTAACGCCGCTTACCTTAATTGTCAGTAATCCGCATc                                                                                                            >  2:113656/1‑251 (MQ=255)                                                                                                                                                       cGAACAGAGCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGAAACCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTCCATGCCGGAGCCTGAAAAAAAAGGGACGATTCTAAAGACGGTTAGCGTGATTGTCAGTCATCCGTATCAGAAACAGttttacgc                                                                                            >  2:145131/1‑247 (MQ=255)                                                                                                                                                                        aGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAA‑AGTTTTTACGCGTTTTTACCTGCGCCAt                                                                          >  2:11821/1‑251 (MQ=255)                                                                                                                                                                               tCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAA‑AGTTTTTACGCGTTTTTACCTGCGCCATGCTTTcc                                                                   >  1:7619/1‑251 (MQ=255)                                                                                                                                                                                   ccATCCTGCCTGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGTATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTTAAAAAAATGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAA‑AGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGaa                                                               <  1:116819/251‑1 (MQ=255)                                                                                                                                                                                    cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAA‑AGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAAt                                                              >  1:22622/1‑251 (MQ=255)                                                                                                                                                                                                ttgcttgcGCTCCAACCCGCCAGGAGAGTACGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTCCATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAA‑AGTTTTTACGCGTTTTTACCTGAGCCATGCTTTCCAGAATACGGTGCTAGtt                                                  >  2:251080/1‑251 (MQ=255)                                                                                                                                                                                                                                    aGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAA‑AGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAAc              >  1:132355/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCCGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATACGTATCAGAAA‑AGTGTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTCCCGGAt  >  1:121645/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CGCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAA‑AGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTCCCGGAT  >  NC_000913/4390087‑4390576

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