Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	540,118	T→C	15.4%	D831D (GAT→GAC)	rpoB →	DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	540,118	0	T	C	15.4%	55.3 / 4.3	26	D831D (GAT→GAC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (12/10);  new base C (2/2);  total (14/12)
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.05e-01

GCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGCACTGAACAAGCTGGACGAAGCCGGTATCGTTTACGTAG  >  NC_002947/539974‑540260                                                                                                                                                 |                                                                                                                                               gCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGAt                                                                                                                                                <  2:311853/145‑1 (MQ=255) gCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGAt                                                                                                                                                >  1:311853/1‑145 (MQ=255)                     aCGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATcc                                                                                                                                <  2:356650/141‑1 (MQ=255)                     aCGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATcc                                                                                                                                >  1:356650/1‑141 (MQ=255)                     aCGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAgg                                                                                                                       <  1:28468/150‑1 (MQ=255)                         tccGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACt                                                                                                                   >  1:94495/2‑150 (MQ=255)                         cccGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACt                                                                                                                   >  1:59276/1‑150 (MQ=255)                                   gACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCtgtgtg                                                                                                         <  1:264069/150‑1 (MQ=255)                                        gggTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGcgcg                                                                                                    <  1:337235/150‑1 (MQ=255)                                                         gTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGcc                                                                                   >  2:284/1‑150 (MQ=255)                                                          tCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGccc                                                                                  >  1:280410/1‑150 (MQ=255)                                                           cAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCa                                                                                 <  2:280410/150‑1 (MQ=255)                                                                            gTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACTCCAAGCTCGGccc                                                                                  <  2:262834/132‑1 (MQ=255)                                                                            gTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGccc                                                                                  >  1:262834/1‑132 (MQ=255)                                                                                 tGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAgg                                                                             <  1:95183/132‑1 (MQ=255)                                                                                 tGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAgg                                                                             >  2:95183/1‑132 (MQ=255)                                                                                     gTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGaa                                                                           <  1:255341/130‑1 (MQ=255)                                                                                     gTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGaa                                                                           >  2:255341/1‑130 (MQ=255)                                                                                                acAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAg                                            <  2:94495/149‑1 (MQ=255)                                                                                                 cAACTTCGAAGACTCCATCTGCCTGTCCGATCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGc                                           <  1:284/150‑1 (MQ=255)                                                                                                    cTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGc                                        >  2:237982/1‑150 (MQ=255)                                                                                                              tCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCGCGGCACAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGCACTGAACAAg                              >  1:140366/1‑150 (MQ=255)                                                                                                                 aTCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCg                                                                  >  1:334668/1‑111 (MQ=255)                                                                                                                 aTCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCg                                                                  <  2:334668/111‑1 (MQ=255)                                                                                                                      ccTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGAACTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATACCGAACGTGGGTGAAGCCGCACTGAACAAGCTGGACGa                      >  2:190577/1‑150 (MQ=255)                                                                                                                                       ttCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCaga                                                                               <  2:266606/76‑1 (MQ=255)                                                                                                                                       ttCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCaga                                                                               >  1:266606/1‑76 (MQ=255)                                                                                                                                          aGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGCACTGAACAAGCTGGACGAAGCCGGTAATGACTACGTAg  >  2:95545/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                               GCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGCACTGAACAAGCTGGACGAAGCCGGTATCGTTTACGTAG  >  NC_002947/539974‑540260

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