Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,852,273	G→C	43.7%	R129R (CGG→CGC) ‡	pncA →	nicotinamidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,852,273	0	G	C	43.7%	‑1.0 / 24.5	22	R129R (CGG→CGC) ‡	pncA	nicotinamidase
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base C (6/3);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

TCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGAT  >  NC_000913/1852141‑1852407                                                                                                                                     |                                                                                                                                       tCACTGTGTGCAGAAGAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCCTTTTTGATAACGGCCGCCGGCAGAAAACct                                                                                                                              >  2:271217/1‑143 (MQ=255)       tgtgCAGAACAGTGAAGGCGCACAATCGTTGTTTTAAGAGCTCCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTGGATAACGGCCGTTGCCAGAAACCCTCTCTCg                                                                                                                        >  2:177077/1‑143 (MQ=255)         cggaGAACAGTGAGGGCGCGCAATTAACTCCGTTACTTCACCAAAAATGGATCGCAGCGGTGTGCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGAt                                                                                                                      <  2:255043/140‑1 (MQ=255)           cAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGa                                                                                                                    <  1:61191/143‑1 (MQ=255)           cAGAACAGAGAAGGCGCACAATTACAGCCGTGACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGa                                                                                                                    <  1:33823/143‑1 (MQ=255)             gAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGGCGGCAGAAAACATCGCGCGATGACt                                                                                                                  >  2:189550/1‑143 (MQ=255)                     aaGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGAGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTAATAACGACCATTGACAGAAAACCTCTCTCGATTACTGATTAcaa                                                                                                          >  2:105968/1‑141 (MQ=255)                              aaTTACATCCGTTACTGCACC‑AAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGATACAGTGCCTTTTTTGATAACGGCCGTCGGCAGACAACCTCCATCGAAGACAGACTACGCTATCAtcaga                                                                                                >  2:34628/1‑139 (MQ=255)                               aTTACATCCTTTACTACACCAAAAATCGATCGCAGAGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGaaa                                                                                                <  1:138083/143‑1 (MQ=255)                                               cACCAAAAAGCGATCGCAGCGGGGTTCCATAAAGGCGAAAATCCTGTAGTTGACAGTTACAGTGCCTGGTTTGATAACGGCCGTTGCCAGAAAACCTCGCTCGAGGACTGGGTACGGGAGCATGAGAGCGAGGGATAGATCga                                                                                >  1:47888/1‑142 (MQ=255)                                                             cGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTa                                                                  <  1:221879/143‑1 (MQ=255)                                                                 tcGGTGGTCCATAAAGCCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGa                                                              <  2:10088/142‑1 (MQ=255)                                                                   ggTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACt                                                            <  1:57478/143‑1 (MQ=255)                                                                       ttCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTAGTg                                                        >  1:76747/1‑143 (MQ=255)                                                                         ccATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCg                                                      <  1:166373/143‑1 (MQ=255)                                                                                 gCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTAGGCGTTCATGAAGTTGGGGTAAGGATCGTTATGGGCCTAGATACTGGATATTGCGTGGAGttt                                              >  1:59697/1‑143 (MQ=255)                                                                                 gCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTCCGTGAAGttt                                              >  1:225886/1‑143 (MQ=255)                                                                                 gCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGttt                                              >  2:25400/1‑143 (MQ=255)                                                                                           tttAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACGTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGATATGAGCCTGGCTGCTGACTATTGCGTGAAGTTTACCGTACtag                                    >  2:61586/1‑141 (MQ=255)                                                                                                ttGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGCATATTGCGTGAAGTTTACCGTGCTGGAcgcg                               >  2:111127/1‑143 (MQ=255)                                                                                                     aGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAGATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTAcc                          >  1:276324/1‑142 (MQ=255)                                                                                                                             ggCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTACACGTgct  >  1:205389/1‑141 (MQ=255)                                                                                                                                     |                                                                                                                                       TCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGAT  >  NC_000913/1852141‑1852407

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