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%	169.8 / NA	48	Y143H (TAC→CAC)	psd	phosphatidylserine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/20);  total (28/20)

ATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTAGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGATTTT  >  NC_000913/4389804‑4390069                                                                                                                                   |                                                                                                                                        aTCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAcgcg                                                                                                                                 <  1:72612/139‑1 (MQ=255) aTCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAcgcg                                                                                                                                 <  2:728262/139‑1 (MQ=255)     cgcgTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGga                                                                                                                             <  2:122436/139‑1 (MQ=255)                    cAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCggtggtggt                                                                                                                                       >  1:820066/1‑114 (MQ=255)                                            ccGTGAGATGGTTAACGGAGAAGAGATCGCCCCGCACGTAGATAATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      <  1:1124011/139‑1 (MQ=255)                                              gTGAGATGGTTAACGGAGAAGAGAACGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGCGTGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCg                                                                                    >  2:766137/1‑139 (MQ=255)                                               tGAGTTGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                                   >  2:1124011/1‑139 (MQ=255)                                               tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      <  2:409119/136‑1 (MQ=255)                                               tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      >  1:409119/1‑136 (MQ=255)                                                      gTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTg                                                                                                 >  1:604423/1‑118 (MQ=255)                                                      gTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTg                                                                                                 <  2:604423/118‑1 (MQ=255)                                                           cGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGa                                                                            >  1:508817/1‑134 (MQ=255)                                                           cGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGa                                                                            <  2:508817/134‑1 (MQ=255)                                                                aagaGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtgcctg                                                                  >  1:1443832/1‑139 (MQ=255)                                                                 agagATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAc                                                                                                                                    <  1:1225174/72‑1 (MQ=255)                                                                 agagATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAc                                                                                                                                    >  2:1225174/1‑72 (MQ=255)                                                                     aTCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCcagc                                                             >  1:1044336/1‑139 (MQ=255)                                                                          ccGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGc                                                        >  2:1001338/1‑139 (MQ=255)                                                                            ggCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCtt                                                      <  1:594172/139‑1 (MQ=255)                                                                            ggCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCtt                                                      <  1:1001338/139‑1 (MQ=255)                                                                                   aGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCt                                               >  1:1514313/1‑139 (MQ=255)                                                                                               gCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAgg                                                                                        <  2:1210752/86‑1 (MQ=255)                                                                                               gCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAgg                                                                                        >  1:1210752/1‑86 (MQ=255)                                                                                                   aaTACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAcg                                                                                                                                   <  2:1297928/39‑1 (MQ=255)                                                                                                   aaTACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAcg                                                                                                                                   >  1:1297928/1‑39 (MQ=255)                                                                                                   aaTACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTg                               >  1:697067/1‑139 (MQ=255)                                                                                                   aaTACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTg                               >  1:1433043/1‑139 (MQ=255)                                                                                                     tACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGc                             <  1:538632/139‑1 (MQ=255)                                                                                                        cGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCgg                                                                                                                                <  1:821218/37‑1 (MQ=255)                                                                                                        cGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCgg                                                                                                                                >  2:821218/1‑37 (MQ=255)                                                                                                          ttGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCa                        >  1:308441/1‑139 (MQ=255)                                                                                                           tgaaCGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGAAg                       >  1:348554/4‑139 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      >  2:858446/1‑58 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      <  1:858446/58‑1 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCg                                                                                    >  2:1111500/1‑76 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCg                                                                                    <  1:1111500/76‑1 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                                   >  2:1598880/1‑77 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                                   <  1:1598880/77‑1 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCAt                                                                                <  1:754809/80‑1 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCAt                                                                                >  2:754809/1‑80 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGAt                                                                           >  1:919333/1‑85 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGAt                                                                           <  2:919333/85‑1 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGc                             <  2:1562150/131‑1 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGc                             >  1:1562150/1‑131 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGa                     >  2:1459619/1‑139 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGa                     >  2:65643/1‑139 (MQ=255)                                                                                                               cGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGAtt                   >  2:521575/1‑139 (MQ=255)                                                                                                                    tgtgTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCAGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTAtc              >  1:37687/1‑139 (MQ=255)                                                                                                                          cgcgGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCttcttc        >  2:1469764/1‑139 (MQ=255)                                                                                                                                ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGAtttt  >  2:1509959/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                        ATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTAGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGATTTT  >  NC_000913/4389804‑4390069

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: ‑