Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,717,873	G→A	85.1%	R131H (CGC→CAC)	fabB →	3‑oxoacyl‑ACP synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,717,873	0	G	A	85.1%	130.2 / 10.3	55	R131H (CGC→CAC)	fabB	3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (22/24);  total (27/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 = 9.94e-01

AAGGACGCTGGCCTGACCGAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGT  >  NC_002947/4717731‑4718010                                                                                                                                                |                                                                                                                                          aaGGACGCTGGCCTGACCTAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTAcc                                                                                                                                    <  1:511024/150‑1 (MQ=255) aaGGACGCTGGCCTGACCGAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTAcc                                                                                                                                    <  1:510995/150‑1 (MQ=255)    gACGCTGGCCTGACCGAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGc                                                                                                                                 >  2:287687/1‑150 (MQ=255)            ccTGACCGAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATggg                                                                                                                         >  2:27127/1‑150 (MQ=255)                 ccGAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGgcagca                                                                                                                    >  1:155304/1‑150 (MQ=255)                  cGAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCAc                                                                                                                   >  1:1109513/1‑150 (MQ=255)                  cGAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCAc                                                                                                                   >  1:704648/1‑150 (MQ=255)                   gAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACTCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCAcc                                                                                                                  <  2:504473/150‑1 (MQ=255)                             tttCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCAt                                                                                                        <  2:50497/150‑1 (MQ=255)                                ccAACCCGCGTACCGGCCTGGTGGCGGGTTCAGGCGGTCCACCCGCTGCCAATCAGATGGAAGCGCTGGACCCAACGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGcc                                                                                                     <  1:287687/150‑1 (MQ=255)                                   aCCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGGGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTgg                                                                                                  >  2:777004/1‑150 (MQ=255)                                   aCCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTgg                                                                                                  <  2:695514/150‑1 (MQ=255)                                    cccGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGGGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGc                                                                                                 <  1:777004/150‑1 (MQ=255)                                           aCCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCg                                                                                          <  1:27127/150‑1 (MQ=255)                                             cGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGACACCCCGtt                                                                                        <  2:155304/150‑1 (MQ=255)                                                  tggtggCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGa                                                                                   <  2:847549/150‑1 (MQ=255)                                                       gCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAg                                                                              <  2:704648/150‑1 (MQ=255)                                                               cggcggCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCCTACCACGTTACCCGCACCATGGGCAGCCCCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCaa                                                                      <  2:1111776/150‑1 (MQ=255)                                                               cggcggCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCaa                                                                      <  2:333356/150‑1 (MQ=255)                                                               cggcggCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCaa                                                                      >  1:1214072/1‑150 (MQ=255)                                                                        cTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGCAAGGCGT‑GCGCGCCGTCGGCCCATAGCTCGGTACCCTCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGAt                                                             >  1:798528/1‑150 (MQ=255)                                                                         tCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGAtc                                                            >  1:1234397/1‑150 (MQ=255)                                                                         tCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATc                                                                        >  2:82733/1‑138 (MQ=255)                                                                         tCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATc                                                                        <  1:82733/138‑1 (MQ=255)                                                                            aCCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAactact                                                                 >  1:74373/1‑142 (MQ=255)                                                                            aCCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAactact                                                                 <  2:74373/142‑1 (MQ=255)                                                                              ccTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAac                                                                     <  2:409602/136‑1 (MQ=255)                                                                              ccTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGAACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAac                                                                     >  1:409602/1‑136 (MQ=255)                                                                               cTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCtcgtcg                                                      <  2:798528/150‑1 (MQ=255)                                                                               cTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCtcgtcg                                                      <  2:791698/150‑1 (MQ=255)                                                                                tGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCgg                                                     >  1:950865/1‑150 (MQ=255)                                                                                    ccAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTg                                                 >  2:761184/1‑150 (MQ=255)                                                                                     cAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTgc                                                <  2:1214072/150‑1 (MQ=255)                                                                                     cAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTgc                                                >  1:598174/1‑150 (MQ=255)                                                                                       gATGGAAGCGCTGGACACCCTGCGTGAGAACGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTgcgc                                              <  2:291695/150‑1 (MQ=255)                                                                                            aaGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCt                                         >  1:333770/1‑150 (MQ=255)                                                                                                gcTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCaa                                                                                     <  1:888322/102‑1 (MQ=255)                                                                                                gcTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCaa                                                                                     >  2:888322/1‑102 (MQ=255)                                                                                                    gACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCAcc                                                                                             >  2:3641/1‑90 (MQ=255)                                                                                                    gACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCAcc                                                                                             <  1:3641/90‑1 (MQ=255)                                                                                                    gACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCAc                                 >  2:205797/1‑150 (MQ=255)                                                                                                         ccTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGAt                                                             >  2:1028210/1‑117 (MQ=255)                                                                                                         ccTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGAt                                                             <  1:1028210/117‑1 (MQ=255)                                                                                                           tGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGgc                                                                                                                        >  2:898746/1‑56 (MQ=255)                                                                                                           tGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGgc                                                                                                                        <  1:898746/56‑1 (MQ=255)                                                                                                                gagaAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATg                                                                                                                           <  2:776376/48‑1 (MQ=255)                                                                                                                gagaAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATg                                                                                                                           >  1:776376/1‑48 (MQ=255)                                                                                                                  gaAAGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGc                   >  1:888512/1‑150 (MQ=255)                                                                                                                   aaaGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCa                  >  2:6113/1‑150 (MQ=255)                                                                                                                     aGGCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACt                <  1:1150186/150‑1 (MQ=255)                                                                                                                       gCGTCAAGCG‑CGTCGGCCCATACCACGTTACCCACACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTgg              >  1:911803/1‑150 (MQ=255)                                                                                                                        cGTCAAGCG‑CGTCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCACCACCTCGGCTCACTGCATCGGTACCGCACTGGa             >  1:1097548/1‑150 (MQ=255)                                                                                                                                    tCGGCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGt  <  2:598174/150‑1 (MQ=255)                                                                                                                                       gCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAg            <  2:712978/137‑1 (MQ=255)                                                                                                                                       gCCCATACCACGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAg            >  1:712978/1‑137 (MQ=255)                                                                                                                                                |                                                                                                                                          AAGGACGCTGGCCTGACCGAAGAGCAGGTTTCCAACCCGCGTACCGGCCTGGTGGCGGGTTCCGGCGGCGCCTCGACCCTGAACCAGATGGAAGCGCTGGACACCCTGCGTGAGAAAGGCGTCAAGCG‑CGTCGGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGT  >  NC_002947/4717731‑4718010

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