Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,579,592	A→G	7.1%	*288W (TGA→TGG)	psd →	phosphatidylserine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,579,592	0	A	G	7.1%	163.3 / 2.4	56	*288W (TGA→TGG)	psd	phosphatidylserine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

CTGGAGAAAGGTGCCGAACTAGGCCGCTTCAAGCTGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCG  >  NC_002947/5579443‑5579738                                                                                                                                                      |                                                                                                                                                   cTGGAGAAAGGTGCCGAACTAGGCCGCTTCAAGCTGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGa                                                                                                                                                    >  2:1100408/1‑150 (MQ=255)              ccGAACTAGGCCGCTTCAAGCTGGGTTCGACCGCCCACGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATg                                                                                                                                       <  1:67129/150‑1 (MQ=255)                    tAGGCCGCTTCAAGCTGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAgg                                                                                                                                 >  2:1120463/1‑150 (MQ=255)                            ttCAAGCTGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCat                                                                                                                         >  1:55612/1‑150 (MQ=255)                                 gCTGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGAtt                                                                                                                                                  <  2:402101/120‑1 (MQ=255)                                 gCTGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGAtt                                                                                                                                                  >  1:402101/1‑120 (MQ=255)                                 gCTGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCa                                                                                                                    >  2:734700/1‑150 (MQ=255)                                   tGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAgg                                                                                                                  <  1:734700/150‑1 (MQ=255)                                    gggTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAggg                                                                                                                 <  2:55612/150‑1 (MQ=255)                                      gTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGcc                                                                                                               >  1:435695/1‑150 (MQ=255)                                      gTTCGACCGCCAGCGTGCTGTTCGTGCCAGCGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCCTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGcc                                                                                                               <  2:1015311/150‑1 (MQ=255)                                      gTTAGACCGCCATCGTGCTGTTCGGGACAGAGCAGGTGAAATGGGGTGACACCCAAAAAGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATAAAAAAAAAAAAAAAAGGCCCTGCATATGCAGGGcc                                                                                                               <  2:358660/150‑1 (MQ=255)                                       ttCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCt                                                                                                              >  2:121037/1‑150 (MQ=255)                                        tCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCtt                                                                                                             <  2:1250730/150‑1 (MQ=255)                                         cGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCttt                                                                                                            <  1:1100408/150‑1 (MQ=255)                                           aCCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGGTTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGcc                                                                                                               >  1:878476/1‑145 (MQ=255)                                           aCCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGGTTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGcc                                                                                                               >  1:877879/1‑145 (MQ=255)                                           aCCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGGTTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGcc                                                                                                               <  2:878476/145‑1 (MQ=255)                                           aCCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAACGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGGTTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGcc                                                                                                               <  2:877879/145‑1 (MQ=255)                                            ccGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTc                                                                                                         >  1:1045894/1‑150 (MQ=255)                                               ccATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGtt                                                                                                      <  2:885780/150‑1 (MQ=255)                                                 aTCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTgg                                                                                                    <  2:435695/150‑1 (MQ=255)                                                      gCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGgccgc                                                                                               >  1:1107417/1‑150 (MQ=255)                                                          ttCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGa                                                                                                                                               >  1:679727/1‑98 (MQ=255)                                                          ttCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGa                                                                                                                                               <  2:679727/98‑1 (MQ=255)                                                           tCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGc                                                                                          <  1:40487/150‑1 (MQ=255)                                                            cGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCg                                                                                         <  2:341120/150‑1 (MQ=255)                                                              ggCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCat                                                                                                                                          >  1:408394/1‑99 (MQ=255)                                                              ggCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCat                                                                                                                                          <  2:408394/99‑1 (MQ=255)                                                                    agCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGgccgc                                                                                               <  2:991318/136‑1 (MQ=255)                                                                    agCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGgccgc                                                                                               >  1:991318/1‑136 (MQ=255)                                                                       aGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCat                                                                                                                         >  2:667653/1‑107 (MQ=255)                                                                       aGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCat                                                                                                                         <  1:667653/107‑1 (MQ=255)                                                                         gTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTgg                                                                                                    <  1:356602/126‑1 (MQ=255)                                                                         gTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTgg                                                                                                    >  2:356602/1‑126 (MQ=255)                                                                         gTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCAc                                                                            >  2:2676/1‑150 (MQ=255)                                                                          tGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTg                                                                                                     >  2:835995/1‑124 (MQ=255)                                                                          tGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTg                                                                                                     <  1:835995/124‑1 (MQ=255)                                                                           gAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCt                                                                                                              <  1:353259/114‑1 (MQ=255)                                                                           gAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCt                                                                                                              >  2:353259/1‑114 (MQ=255)                                                                              aTGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCAGCCTTGCGGCCCCGGATTCACGGAtg                                                                       >  1:266528/1‑149 (MQ=255)                                                                                 ggCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGc                                                                    <  2:1107417/150‑1 (MQ=255)                                                                                    tGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCcgc                                                                 >  1:930325/1‑150 (MQ=255)                                                                                            tGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGtg                                                          >  1:661123/1‑148 (MQ=255)                                                                                            tGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTc                                                          <  2:661123/149‑1 (MQ=255)                                                                                                   aGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCAc                                                  >  2:774221/1‑150 (MQ=255)                                                                                                         tGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGc                                            >  1:1255503/1‑150 (MQ=255)                                                                                                              ttCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGAGCAGCCGCGGGCGTCGCGATCACGCAGGCCcagc                                       <  1:1264648/150‑1 (MQ=255)                                                                                                                   aTGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTa                                  <  2:1045894/150‑1 (MQ=255)                                                                                                                     ggggCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTACa                                >  2:293822/1‑150 (MQ=255)                                                                                                                        gcagcaGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTACAGCa                             >  1:768814/1‑150 (MQ=255)                                                                                                                                  gcggcgCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTg                                                                                           >  2:833843/1‑78 (MQ=255)                                                                                                                                  gcggcgCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTg                                                                                           <  1:833843/78‑1 (MQ=255)                                                                                                                                   cggcgCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGc                  <  1:774221/150‑1 (MQ=255)                                                                                                                                                 gCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTAc                                 >  2:823513/1‑121 (MQ=255)                                                                                                                                                 gCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTAc                                 <  1:823513/121‑1 (MQ=255)                                                                                                                                                   ttGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCg  >  1:1259998/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CTGGAGAAAGGTGCCGAACTAGGCCGCTTCAAGCTGGGTTCGACCGCCATCGTGCTGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCG  >  NC_002947/5579443‑5579738

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