Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,948,216	C→G	21.4%	R494R (CGG→CGC)	PP_1749 ←	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,948,216	0	C	G	21.4%	57.0 / 12.3	28	R494R (CGG→CGC)	PP_1749	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base G (3/3);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

GCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGA  >  NC_002947/1948068‑1948358                                                                                                                                                     |                                                                                                                                               gCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGcgcc                                                                                                                                               >  1:182550/1‑150 (MQ=255)  cTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGcgccg                                                                                                                                              <  1:386830/150‑1 (MQ=255)       cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTCCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACgg                                                                                                                                         >  2:4188/1‑150 (MQ=255)           aTCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAg                                                                                                                                     <  1:543046/150‑1 (MQ=255)            tCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGGCGCACGGCCaac                                                                                                                                    <  1:462211/150‑3 (MQ=255)            tCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAgg                                                                                                                                    >  2:85128/1‑150 (MQ=255)                gCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAgcgcct                                                                                                                                >  1:551751/1‑145 (MQ=255)                gCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCt                                                                                                                                <  1:424471/150‑1 (MQ=255)                    cATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGGCGCACGGCCaac                                                                                                                                    >  2:495103/1‑140 (MQ=255)                    cATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGGCGCACGGCCaac                                                                                                                                    <  1:495103/142‑3 (MQ=255)                     aTGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGGCGCACGGCCa                                                                                                                                      <  1:181939/139‑1 (MQ=255)                     aTGAAGTACAGCACCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGGCGCACGGCCa                                                                                                                                      >  2:181939/1‑139 (MQ=255)                             cAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCaa                                                                                                                   >  1:228188/1‑150 (MQ=255)                              aGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAac                                                                                                                  <  1:438089/150‑1 (MQ=255)                                  ccaccaCCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgt                                                                                                              <  1:85128/150‑1 (MQ=255)                                   caccacTGGAATTTCCAGTGCCCGTGCCGTGCGCACGGCAGCATCGGCCAGCACCGGGTGCTGGCGTTCGGTGACGCCCTCCCGGGTGCCGCCGGTGTGCAGGTTGGCCGCGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgtc                                                                                                             <  2:369046/150‑1 (MQ=255)                                    accaccGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgtcg                                                                                                            >  2:41382/1‑150 (MQ=255)                                     ccaccGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGCCCGGCAACAcgtcgt                                                                                                           >  2:99129/1‑150 (MQ=255)                                              tttCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCaa                                                                                                  >  2:353070/1‑150 (MQ=255)                                                  cAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCaa                                                                                                                   <  2:402269/129‑1 (MQ=255)                                                  cAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCaa                                                                                                                   >  1:402269/1‑129 (MQ=255)                                                               gcgcgcACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGTAAg                                                                                 <  2:426487/150‑1 (MQ=255)                                                                 gcgcACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGt                                                                               <  1:424849/150‑1 (MQ=255)                                                                           gcaTCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGtt                                                                     <  1:90884/150‑1 (MQ=255)                                                                                         cGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGAACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCCCGGTTTCGTCGTCCAGCgg                                                       >  2:284441/1‑150 (MQ=255)                                                                                                              gTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGt                                  >  2:505369/1‑150 (MQ=255)                                                                                                                         tGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGc                       <  2:374057/150‑1 (MQ=255)                                                                                                                                      cAGGTTGGCCGTGCGCCGCACGGCCAGGCGTTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTg          >  1:456792/1‑150 (MQ=255)                                                                                                                                        gcctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTgcg        >  2:237382/4‑150 (MQ=255)                                                                                                                                              ccGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGACACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGa  >  1:305259/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                               GCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGA  >  NC_002947/1948068‑1948358

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