Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	466,444	A→C	6.6%	D148A (GAC→GCC)	PP_0385 →	serine acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	466,444	0	A	C	6.6%	177.6 / 4.8	61	D148A (GAC→GCC)	PP_0385	serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base C (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

CCCGGTACGTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCAGAGGGGA  >  NC_002947/466297‑466593                                                                                                                                                    |                                                                                                                                                      cccGGTACGTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCGCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGaga                                                                                                                                                     >  1:1583174/1‑148 (MQ=255)      tACGTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGTCCCTGGGGTTGCAGGAGTACGCACTGTGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGc                                                                                                                                                <  2:300267/150‑1 (MQ=255)         gTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGAt                                                                                                                                             <  2:531650/150‑1 (MQ=255)                      gcAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCg                                                                                                                                <  2:158313/150‑1 (MQ=255)                        aaCTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGgc                                                                                                                              >  2:225499/1‑150 (MQ=255)                             cttcATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGAt                                                                                                                         >  1:654179/1‑150 (MQ=255)                             cttcATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGAGATCGGCGCGAt                                                                                                                         <  2:1126119/150‑1 (MQ=255)                                   ccGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                   <  2:788840/150‑1 (MQ=255)                                    cGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTg                                                                                                                  <  2:646942/150‑1 (MQ=255)                                    cGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTg                                                                                                                  >  2:971838/1‑150 (MQ=255)                                     ggCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                   >  1:1443881/1‑148 (MQ=255)                                     ggCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                   <  2:1443881/148‑1 (MQ=255)                                        ggAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAg                                                                                                                                      <  1:1332831/125‑1 (MQ=255)                                        ggAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAg                                                                                                                                      >  2:1332831/2‑126 (MQ=255)                                            cTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTCGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTAcgcg                                                                                                          <  1:1567883/150‑1 (MQ=255)                                              gCTCGATCGGCATCAAGCCCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCggg                                                                                                        <  1:563467/150‑1 (MQ=255)                                              gCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGa                                                                                                                                              <  2:695340/112‑1 (MQ=255)                                              gCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGa                                                                                                                                              >  1:695340/1‑112 (MQ=255)                                               cTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGa                                                                                                       >  1:95971/1‑150 (MQ=255)                                                  gATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGatat                                                                                                    >  1:712797/1‑150 (MQ=255)                                                       gCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGc                                                                                               >  1:1083912/1‑150 (MQ=255)                                                        cATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGa                                                                                                                                                       >  1:290394/1‑93 (MQ=255)                                                        cATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGa                                                                                                                                                       <  2:290394/93‑1 (MQ=255)                                                           cAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTgg                                                                                           >  1:226468/1‑150 (MQ=255)                                                                     gggTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTc                                                                                 <  2:1083912/150‑1 (MQ=255)                                                                            aGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAAcc                                                                          <  1:971838/150‑1 (MQ=255)                                                                                    gCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGa                                                                                                                          <  2:524911/94‑1 (MQ=255)                                                                                    gCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGa                                                                                                                          >  1:524911/1‑94 (MQ=255)                                                                                    gCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAg                                                                  >  1:282045/1‑150 (MQ=255)                                                                                      aCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                   >  2:566197/1‑99 (MQ=255)                                                                                      aCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                   <  1:566197/99‑1 (MQ=255)                                                                                      aCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGc                                                                                               >  1:75405/1‑119 (MQ=255)                                                                                      aCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGc                                                                                               <  2:75405/119‑1 (MQ=255)                                                                                       cTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTc                                                               <  2:654179/150‑1 (MQ=255)                                                                                          cgcATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCgag                                                            >  1:1163703/1‑150 (MQ=255)                                                                                           gcATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCgaga                                                           >  1:586249/1‑150 (MQ=255)                                                                                           gcATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCgaga                                                           <  2:1601998/150‑1 (MQ=255)                                                                                           gcATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCgaga                                                           <  1:225499/150‑1 (MQ=255)                                                                                                 gTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTg                                                                                                                  >  2:659441/1‑89 (MQ=255)                                                                                                 gTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTg                                                                                                                  <  1:659441/89‑1 (MQ=255)                                                                                                   gACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATGTTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCg                                                   <  2:83388/150‑1 (MQ=255)                                                                                                      aaCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGCCAATGTCGATgc                                                                                                                                           >  1:1602164/1‑58 (MQ=255)                                                                                                      aaCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGCCAATGTCGATgc                                                                                                                                           <  2:777009/59‑2 (MQ=37)                                                                                                      aaCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGCCAATGTCGATgc                                                                                                                                           >  1:777009/1‑58 (MQ=255)                                                                                                      aaCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGCCAATGTCGATgc                                                                                                                                           <  2:1602164/59‑2 (MQ=37)                                                                                                       aCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAg                                               <  2:95971/150‑1 (MQ=255)                                                                                                        cGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAgg                                              <  2:282045/150‑1 (MQ=255)                                                                                                              cATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTgg                                                                                                                 >  2:698192/1‑77 (MQ=255)                                                                                                              cATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTgg                                                                                                                 <  1:698192/77‑1 (MQ=255)                                                                                                              cATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGAGTTCGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTa                                        >  2:1050454/1‑150 (MQ=255)                                                                                                                tGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAAc                                      <  2:586249/150‑1 (MQ=255)                                                                                                                   gcgcCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCaa                                   >  2:1498084/1‑150 (MQ=255)                                                                                                                       cAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTgg                                                                                                                 <  1:155322/68‑1 (MQ=255)                                                                                                                       cAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTgg                                                                                                                 >  2:155322/1‑68 (MQ=255)                                                                                                                               gCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGa                                                                                                       <  1:31335/70‑1 (MQ=255)                                                                                                                               gCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGa                                                                                                       >  2:31335/1‑70 (MQ=255)                                                                                                                                         tGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGCTACATTTGCGGACGtt             >  2:1109789/1‑150 (MQ=255)                                                                                                                                          gATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGc                    <  2:1228698/142‑1 (MQ=255)                                                                                                                                          gATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGc                    >  1:1228698/1‑142 (MQ=255)                                                                                                                                           aTCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTgg           <  1:97580/150‑1 (MQ=255)                                                                                                                                            tCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGc          <  2:712797/150‑1 (MQ=255)                                                                                                                                              gCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCag        >  1:369398/1‑150 (MQ=255)                                                                                                                                               cATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCaga       >  1:1029537/1‑150 (MQ=255)                                                                                                                                                aTCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCagag      >  2:986415/1‑150 (MQ=255)                                                                                                                                                    aCATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTATATTTGCGGACGTTGGCAGAGGGGa  <  1:1498084/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CCCGGTACGTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCAGAGGGGA  >  NC_002947/466297‑466593

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