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

ATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTAGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGA  >  NC_000913/4389804‑4390065                                                                                                                                   |                                                                                                                                    aTCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAcgcg                                                                                                                             <  2:144477/139‑1 (MQ=255)            gttACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTaa                                                                                                                  <  1:1468208/139‑1 (MQ=255)                    cAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAg                                                                                                                 <  2:456328/132‑1 (MQ=255)                    cAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAg                                                                                                                 >  1:456328/1‑132 (MQ=255)                    cAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGACCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTcac                                                                                                          >  1:963822/1‑139 (MQ=255)                      aaCAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTGCGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACaa                                                                                                        >  1:820164/1‑139 (MQ=255)                             tCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCgg                                                                                                                            >  2:514408/1‑112 (MQ=255)                             tCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCgg                                                                                                                            <  1:514408/112‑1 (MQ=255)                                     ttCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGa                                                                                         >  1:93183/1‑139 (MQ=255)                                            ccGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                  <  2:963822/139‑1 (MQ=255)                                               tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                               <  1:1069140/139‑1 (MQ=255)                                               tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                               >  2:751579/1‑139 (MQ=255)                                               tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                               >  1:421207/1‑139 (MQ=255)                                                      gTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTaa                                                                                                                  >  2:1014768/1‑97 (MQ=255)                                                      gTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTaa                                                                                                                  <  1:1014768/97‑1 (MQ=255)                                                              agaagaGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTg                                                                                             >  1:353536/1‑110 (MQ=255)                                                              agaagaGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTg                                                                                             <  2:353536/110‑1 (MQ=255)                                                                 agagATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTa                                                                                                                   >  2:569460/1‑85 (MQ=255)                                                                 agagATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTa                                                                                                                   <  1:569460/85‑1 (MQ=255)                                                                    gaTCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCcag                                                          <  2:820164/139‑1 (MQ=255)                                                                     aTCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGTGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCcagc                                                         >  2:85072/1‑139 (MQ=255)                                                                        gCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGga                                                                                                                         >  2:312131/1‑72 (MQ=255)                                                                        gCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGga                                                                                                                         <  1:312131/72‑1 (MQ=255)                                                                        gCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAgg                                                      >  1:236702/1‑139 (MQ=255)                                                                                    gATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTg                                          <  1:85072/139‑1 (MQ=255)                                                                                             aCGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAc                                                                                                                                <  1:1566856/44‑1 (MQ=255)                                                                                             aCGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCAc                                                                                                                                >  2:1566856/1‑44 (MQ=255)                                                                                             aCGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                  >  1:546231/1‑90 (MQ=255)                                                                                             aCGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                  <  2:546231/90‑1 (MQ=255)                                                                                               gCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGcc                               >  1:1399408/1‑139 (MQ=255)                                                                                                cAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCt                              <  2:236702/139‑1 (MQ=255)                                                                                                 aGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCtt                             >  1:1253550/1‑139 (MQ=255)                                                                                                  gAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCttt                            >  2:755074/1‑139 (MQ=255)                                                                                                     tACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATGGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGc                         <  1:673060/139‑1 (MQ=255)                                                                                                         gTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTc                                                 >  2:16927/1‑111 (MQ=255)                                                                                                         gTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTc                                                 <  1:16927/111‑1 (MQ=255)                                                                                                           tGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCCGCCAGCAGGGCTTc                                                 >  1:1197317/1‑109 (MQ=255)                                                                                                           tGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCCGCCAGCAGGGCTTc                                                 <  2:1197317/109‑1 (MQ=255)                                                                                                             cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGa                 >  2:376170/1‑139 (MQ=255)                                                                                                               cGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACATACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCATGCCAGCAGGGCTTCGAGGCTGTAGTtg                                    <  2:999720/118‑1 (MQ=255)                                                                                                               cGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTc                                                 >  2:25739/1‑105 (MQ=255)                                                                                                               cGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTc                                                 <  1:25739/105‑1 (MQ=255)                                                                                                               cGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTtg                                    >  1:999720/1‑118 (MQ=255)                                                                                                                    tgtgTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTAtc          >  1:830809/1‑139 (MQ=255)                                                                                                                            cgGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGACCAGGTCCGCCATCAGATAGTTGCCTGCCACCAGGGCTTCGGGGCTGTCGTTGTGGCCTTTGGCTTGCAGGGTTTTATCTTCTTCgc  >  1:149436/1‑138 (MQ=255)                                                                                                                            cgGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGa  >  1:794627/1‑139 (MQ=255)                                                                                                                            cgGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGCtgtctcttatac                          >  2:690582/1‑106 (MQ=255)                                                                                                                            cgGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGCtgtctcttatac                          <  1:690582/115‑10 (MQ=255)                                                                                                                                ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATcttctt     <  2:1615320/132‑1 (MQ=255)                                                                                                                                ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGCTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAgg                  <  2:993983/119‑1 (MQ=255)                                                                                                                                   |                                                                                                                                    ATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTAGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGA  >  NC_000913/4389804‑4390065

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