Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	892,896	G→C	10.4%	S274R (AGC→AGG)	pta ←	phosphate acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	892,896	0	G	C	10.4%	140.8 / 7.8	48	S274R (AGC→AGG)	pta	phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (23/20);  new base C (2/3);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

TTGGCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTG  >  NC_002947/892747‑893045                                                                                                                                                      |                                                                                                                                                      ttGGCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAgg                                                                                                                                                       <  2:1098170/150‑1 (MQ=255)   ggCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCt                                                                                                                                                     <  1:587110/150‑1 (MQ=255)   ggCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCt                                                                                                                                                     <  1:488811/150‑1 (MQ=255)     cGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGAGACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGa                                                                                                                                                   <  1:1385141/150‑1 (MQ=255)      ggTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGAc                                                                                                                                                  >  2:1123898/1‑150 (MQ=255)           cGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCCGCATCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGa                                                                                                                                             <  1:1105480/150‑1 (MQ=255)           cGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGTCAGCCTGGCAGCGa                                                                                                                                             <  1:773832/150‑1 (MQ=255)           cGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGa                                                                                                                                             <  2:669123/150‑1 (MQ=255)               cGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgat                                                                                                                                         <  1:1123898/150‑1 (MQ=255)               agcgCCGGTGCTGACCGAAAGGATCGGCAGGGCGGCCTGCAGGGCACCGCGGCAGAGGCCAAGGATACGCCCGTCAGGCTTGCTGTTGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCGAGGCTGACAGCGAGgat                                                                                                                                         <  2:469309/147‑1 (MQ=255)                gAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgatg                                                                                                                                        <  1:443962/150‑1 (MQ=255)                  gCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgatgat                                                                                                                                      <  2:549953/150‑1 (MQ=255)                  gCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCAGTGATCGCTGCCAGGCTGACAGCGAGgatgat                                                                                                                                      <  1:1402999/150‑1 (MQ=255)                              cGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCTTCGATTGGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTc                                                                                                                          >  1:1156510/1‑150 (MQ=255)                                aaaGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGc                                                                                                                        >  2:491172/1‑150 (MQ=255)                                     aTCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCggg                                                                                                                   <  1:1145334/150‑1 (MQ=255)                                      cctgCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCgggg                                                                                                                  >  2:635280/4‑150 (MQ=255)                                             gccggccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTaccac                                                                                                           >  1:452556/1‑150 (MQ=255)                                               cggccTGCAGGGCACCGCGGCACAGGCCAAGCATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGACAGGGGGACGCCATTGATCGCTGCCAGGCTCACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTAccacca                                                                                                         >  2:646355/1‑150 (MQ=255)                                                    tGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTg                                                                                                    >  2:1104316/1‑150 (MQ=255)                                                         ggCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGa                                                                                               <  2:720556/150‑1 (MQ=255)                                                              cgcgGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGAAAGCGTGGTACGGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCa                                                                                          <  2:1107251/150‑1 (MQ=255)                                                                           aaGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCgt                                                                             >  2:201919/1‑150 (MQ=255)                                                                                   gCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCa                                                                     >  1:1149306/1‑150 (MQ=255)                                                                                       gTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgt                                                                 >  2:1165095/1‑150 (MQ=255)                                                                                        tCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgtg                                                                >  1:214829/1‑150 (MQ=255)                                                                                         cAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgtgc                                                               <  2:1119460/150‑1 (MQ=255)                                                                                           ctctggCTGTCGCGGGTCAGCAG‑‑‑GCGGGGCAGTGGGACGCCATTGATCGCTGCGAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCGGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGca                                                          <  2:82747/145‑1 (MQ=255)                                                                                              ttGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGca                                                          >  1:682831/1‑150 (MQ=255)                                                                                              ttGCTGTCGCTGGT‑‑‑CAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAg                                                       >  1:82747/1‑150 (MQ=255)                                                                                                   gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgat                                                     >  1:264031/1‑150 (MQ=255)                                                                                                   gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgat                                                     >  2:40225/1‑150 (MQ=255)                                                                                                     cGCCGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatga                                                   >  1:49607/1‑150 (MQ=255)                                                                                                            cagcagCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCtt                                            <  1:478849/150‑1 (MQ=255)                                                                                                             agcagcAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTg                                           >  2:832764/1‑150 (MQ=255)                                                                                                              gcagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGc                                          >  1:527975/1‑150 (MQ=255)                                                                                                               cagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAg                                                                                                        >  2:1148351/1‑87 (MQ=255)                                                                                                               cagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAg                                                                                                        <  1:1148351/87‑1 (MQ=255)                                                                                                                           aaGGGGGGCGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGGCGTCGCGGTCGCCCGAGGTTACCACCAGGGAGCCCGAGGCCAGTAGCGGGCCCGTATTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCgctg                             >  2:89487/2‑150 (MQ=255)                                                                                                                                    gCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATc                    >  2:476873/1‑150 (MQ=255)                                                                                                                                     ccATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCg                   >  2:681061/1‑150 (MQ=255)                                                                                                                                       aTTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTAACACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGcc                 <  1:491172/150‑1 (MQ=255)                                                                                                                                           aTCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCa             <  2:527975/150‑1 (MQ=255)                                                                                                                                            tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatgatgat                                               <  2:33396/115‑1 (MQ=255)                                                                                                                                            tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatgatgat                                               >  1:33396/1‑115 (MQ=255)                                                                                                                                            tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCAGGACCGTGtt                                                                          <  2:794453/88‑1 (MQ=255)                                                                                                                                            tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCAGGACCGTGtt                                                                          >  1:794453/1‑88 (MQ=255)                                                                                                                                            tCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAt            >  2:843643/1‑150 (MQ=255)                                                                                                                                             cGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAtt           <  1:832764/150‑1 (MQ=255)                                                                                                                                                     ggCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAtt           >  2:1213666/1‑142 (MQ=255)                                                                                                                                                     ggCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCAtt           <  1:1213666/142‑1 (MQ=255)                                                                                                                                                      gCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCGGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTg  >  2:1262634/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TTGGCGGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTG  >  NC_002947/892747‑893045

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