Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,724,095	C→A	7.4%	intergenic (+13/‑33)	pssA → / → yfiM	phosphatidylserine synthase; CDP‑diacylglycerol‑serine O‑phosphatidyltransferase/putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,724,095	0	C	A	7.4%	181.3 / 4.0	54	intergenic (+13/‑33)	pssA/yfiM	phosphatidylserine synthase; CDP‑diacylglycerol‑serine O‑phosphatidyltransferase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (27/23);  new base A (2/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTT  >  NC_000913/2723961‑2724242                                                                                                                                        |                                                                                                                                                    ccaccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTACAGGACGGGgttgt                                                                                                                                       <  1:35240/149‑3 (MQ=255)  caccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAAACCCCGTACAGGACGGGgttg                                                                                                                                        <  3:144392/148‑2 (MQ=255)  caccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGgtttgt                                                                                                                                      >  2:61128/1‑149 (MQ=255)     caTCGTTAAGCACTATCGCGATCTGCAAGGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGtttt                                                                                                                                   >  1:88993/1‑149 (MQ=255)      aTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGtttt                                                                                                                                   >  4:166825/1‑148 (MQ=255)       tCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGttttt                                                                                                                                  <  4:329531/148‑1 (MQ=255)           tAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGa                                                                                                                              <  3:48862/148‑1 (MQ=255)             aGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAgg                                                                                                                            <  4:114706/148‑1 (MQ=255)                                aaGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttt                                                                                                        >  2:101596/1‑149 (MQ=255)                                          gATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGtttttt                                                                                                                                 <  2:331783/114‑1 (MQ=255)                                          gATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGtttttt                                                                                                                                 >  1:331783/1‑114 (MQ=255)                                             tATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCt                                                                                           <  2:309728/149‑1 (MQ=255)                                               tCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                      <  2:38812/136‑1 (MQ=255)                                               tCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                      >  1:38812/1‑136 (MQ=255)                                                ccGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTAt                                                                                         <  4:331061/148‑1 (MQ=255)                                                    tGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTa                                                                                    >  2:51598/1‑149 (MQ=255)                                                          ttCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTc                                                                               <  3:297349/148‑1 (MQ=255)                                                          ttCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTc                                                                               >  4:16695/1‑148 (MQ=255)                                                          ttCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTc                                                                               >  4:329110/1‑148 (MQ=255)                                                               aaaCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                      <  2:208936/120‑1 (MQ=255)                                                               aaaCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTGTTGCGTATTCttttt                                                                                                      >  1:208936/1‑120 (MQ=255)                                                                   tCATCCGCCGTTTGCGACGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGc                                                                     >  1:341373/1‑149 (MQ=255)                                                                     aTCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCa                                                                   >  2:307822/1‑149 (MQ=255)                                                                      tCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCat                                                                  >  1:282734/1‑149 (MQ=255)                                                                        cgccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCata                                                                 >  3:195544/1‑148 (MQ=255)                                                                          ccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATg                                                               <  4:185315/148‑1 (MQ=255)                                                                              ttGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTaa                                                          <  2:282734/149‑1 (MQ=255)                                                                                      tATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGc                                                   <  3:21399/148‑1 (MQ=255)                                                                                                gACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCa                                        >  2:212146/1‑149 (MQ=255)                                                                                                gACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCa                                        >  2:209710/1‑149 (MQ=255)                                                                                                   cGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGa                                     <  1:61128/149‑1 (MQ=255)                                                                                                    gATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGa                                                                           >  1:76954/1‑110 (MQ=255)                                                                                                    gATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGa                                                                           >  1:367660/1‑110 (MQ=255)                                                                                                    gATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGa                                                                           <  2:76954/110‑1 (MQ=255)                                                                                                    gATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGa                                                                           <  2:367660/110‑1 (MQ=255)                                                                                                     attaattaGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATa                                   >  1:334181/1‑149 (MQ=255)                                                                                                       taattaGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATaaa                                 <  2:400063/149‑1 (MQ=255)                                                                                                         attaGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGc                               >  2:283110/1‑149 (MQ=255)                                                                                                          ttaGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCt                              >  1:336084/1‑149 (MQ=255)                                                                                                            aGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTATGGGGTTTGTTTTTTCGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCTGGCAGGTTAAAGCTCa                            >  2:352777/1‑149 (MQ=255)                                                                                                               cGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAac                          <  4:195544/148‑1 (MQ=255)                                                                                                               cGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAac                          <  4:231416/148‑1 (MQ=255)                                                                                                                gCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAaca                         <  3:329110/148‑1 (MQ=255)                                                                                                                      tGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATc                  >  2:343013/1‑149 (MQ=255)                                                                                                                           tCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTc             >  1:390601/1‑149 (MQ=255)                                                                                                                            cac‑aACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTc             >  4:211325/1‑148 (MQ=255)                                                                                                                             ac‑aACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCg            >  3:200337/1‑148 (MQ=255)                                                                                                                             ac‑aACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCAGGCAGGATAAAGCTCAACACTTTATCGCCTCgg           >  2:217490/1‑149 (MQ=255)                                                                                                                             ac‑aACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTAATCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCgg           <  2:336084/149‑1 (MQ=255)                                                                                                                              c‑aACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGc          <  1:101596/149‑1 (MQ=255)                                                                                                                                aaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGc          >  4:353591/1‑148 (MQ=255)                                                                                                                                  ccccGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGAt       >  2:98484/1‑149 (MQ=255)                                                                                                                                   cccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATg      <  1:51598/149‑1 (MQ=255)                                                                                                                                       tCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCttt  <  1:137758/149‑1 (MQ=255)                                                                                                                                        |                                                                                                                                                    CCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCAC‑AACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGCTTT  >  NC_000913/2723961‑2724242

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