Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,949,591	A→T	7.3%	K125I (AAA→ATA)	PP_2583 →	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,949,591	0	A	T	7.3%	171.0 / 4.3	55	K125I (AAA→ATA)	PP_2583	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (23/28);  new base T (2/2);  total (25/30)
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.97e-01

TTTCATCAACAGTGGCTGCGCCCACGAGGTCGTGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCG  >  NC_002947/2949460‑2949730                                                                                                                                    |                                                                                                                                            tttCATCAACAGTGGCTGCGCCCACGAGGTCGTGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGtt                                                                                                                           <  1:813181/150‑1 (MQ=255)        aaCAGTGGCTGCGCCCACGAGGTCGTGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAAc                                                                                                                    >  2:1189267/1‑150 (MQ=255)             tGGCTGCGCCCACGAGGTCGTGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCaaa                                                                                                                                <  1:1096478/133‑1 (MQ=255)             tGGCTGCGCCCACGAGGTCGTGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCaaa                                                                                                                                >  2:1096478/1‑133 (MQ=255)                  gcgcCCACGAGGTCGTGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCg                                                                                                          <  2:597674/150‑1 (MQ=255)                                 tGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTtgat                                                                                                                        >  2:1268693/1‑121 (MQ=255)                                 tGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTtgat                                                                                                                        <  1:1268693/121‑1 (MQ=255)                                    tGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCAGCAGGGCAAATGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGcggcgg                                                                                        <  1:596480/150‑1 (MQ=255)                                    tGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGcggcgg                                                                                        <  2:341790/150‑1 (MQ=255)                                    tGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGcggcgg                                                                                        <  1:1128940/150‑1 (MQ=255)                                      tCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCATAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGc                                                                                      >  1:992231/1‑150 (MQ=255)                                         gacgaGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGaa                                                                                                                 >  2:775366/1‑120 (MQ=255)                                         gacgaGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGaa                                                                                                                 <  1:775366/120‑1 (MQ=255)                                         gacgaGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCggg                                                                                       <  2:803275/146‑1 (MQ=255)                                         gacgaGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCggg                                                                                       >  1:803275/1‑146 (MQ=255)                                               gCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATg                                                                                                  >  2:993132/1‑129 (MQ=255)                                               gCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATg                                                                                                  <  1:993132/129‑1 (MQ=255)                                               gCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGa                                                                               <  2:247214/148‑1 (MQ=255)                                               gCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGa                                                                               >  1:247214/1‑148 (MQ=255)                                                 tttCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGAGCATGCACGCGGCGGGCAAGCAGACCAt                                                                           <  1:1189267/150‑1 (MQ=255)                                                 tttCAAAGAGCTGGTCGGCAATGAACCGGCGGCTGCGAACGTGGGGATCATGTCGGAGGTTGTAGACCGGGACCACCAGGGCAAAGGCTAGGCCAAGGTGTGGATGAACGAATTCAGCGGGCGCATGCACGCGGCGGGCAAGCAGACCAt                                                                           >  2:134082/1‑150 (MQ=255)                                                   tCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTc                                                                         >  2:290226/1‑150 (MQ=255)                                                         agCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCATAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACATGa                                                                   <  2:992231/150‑1 (MQ=255)                                                         agCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGa                                                                   <  1:769697/150‑1 (MQ=255)                                                            tGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCAGCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCAcc                                                                      >  1:252649/1‑144 (MQ=255)                                                            tGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCAcc                                                                      <  2:252649/144‑1 (MQ=255)                                                            tGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGAtgt                                                                >  2:677846/1‑150 (MQ=255)                                                                   ccATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCATAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCaa                                                            <  1:228569/147‑1 (MQ=255)                                                                   ccATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCATAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCaa                                                            >  2:228569/1‑147 (MQ=255)                                                                    cATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGAt                                                                  <  1:189952/140‑1 (MQ=255)                                                                    cATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGAt                                                                  >  2:189952/1‑140 (MQ=255)                                                                      tGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGagca                                                      <  1:345017/150‑1 (MQ=255)                                                                        aCCCGGCGGCCGCGAACGTCGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTAATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGagcagc                                                    <  2:843015/150‑1 (MQ=255)                                                                                gcCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGcc                                            >  2:492475/1‑150 (MQ=255)                                                                                  cgcgAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCt                                          <  2:497503/150‑1 (MQ=255)                                                                                     gAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTa                                       >  2:789674/1‑150 (MQ=255)                                                                                     gAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTa                                       <  1:789674/150‑1 (MQ=255)                                                                                        cGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCa                                    <  1:290226/150‑1 (MQ=255)                                                                                             tgATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCAt                                                                                                   <  1:1086063/82‑1 (MQ=255)                                                                                             tgATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCAt                                                                                                   >  2:1086063/1‑82 (MQ=255)                                                                                               aTCATGTCGGTGGTGGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTg                             <  1:501257/150‑1 (MQ=255)                                                                                                 cATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTggg                           <  2:631814/150‑1 (MQ=255)                                                                                                     tCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTAc                       <  1:256803/150‑1 (MQ=255)                                                                                                        gtggtTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTACCGa                    >  2:1046618/1‑150 (MQ=255)                                                                                                        gtggtTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTACCGa                    >  1:512936/1‑150 (MQ=255)                                                                                                           gtTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTACCGATAc                 >  1:241228/1‑150 (MQ=255)                                                                                                             tGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTa                                       <  1:436909/126‑1 (MQ=255)                                                                                                             tGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTa                                       >  2:436909/1‑126 (MQ=255)                                                                                                                    ccGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCggg                                                                                       >  2:1224251/1‑71 (MQ=255)                                                                                                                    ccGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTACCGATACGt               >  1:905392/1‑143 (MQ=255)                                                                                                                    ccGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTACCGATACGt               <  2:905392/143‑1 (MQ=255)                                                                                                                    ccGGCCCACCAGGGCAAAGGCTACGCCAAAATGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCggg                                                                                       <  1:1224251/71‑1 (MQ=255)                                                                                                                          caccaGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTACCGATACGTGAAAccgtcgccg  >  2:428447/1‑150 (MQ=255)                                                                                                                                  cAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGAtgt                                                                >  1:923515/1‑80 (MQ=255)                                                                                                                                  cAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGAtgt                                                                <  2:923515/80‑1 (MQ=255)                                                                                                                                    |                                                                                                                                            TTTCATCAACAGTGGCTGCGCCCACGAGGTCGTGATGTCCGACGAGGCTTTCAAAGAGCTGGTCGGCCATGACCCGGCGGCCGCGAACGTGGTGATCATGTCGGTGGTTGTAGACCCGGCCCACCAGGGCAAAGGCTACGCCAAAGTGTTGATGAACGAATTCATCGGGCGCATGCACGCGGCGGGCAAGCAGACCATTCACCTGATGTGCAAGGAGCAGCATGTAGCCCTGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCG  >  NC_002947/2949460‑2949730

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