Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,852,271	C→T	R129C (CGG→TGC)	pncA →	nicotinamidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,852,271	0	C	T	97.3%	86.9 / NA	37	R129C (CGG→TGC)	pncA	nicotinamidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/16);  minor base G (0/1);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGA  >  NC_000913/1852142‑1852406                                                                                                                                  |                                                                                                                                        cACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTGCATAAAGGCGACAATCCTTTAGTTGACAGTTGCAGTGCCTTTTTTGATAACGGCCGGTGGCAGAAAACttc                                                                                                                            >  2:694245/1‑143 (MQ=255) cACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACctc                                                                                                                            >  1:475087/1‑143 (MQ=255)        tgCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGAt                                                                                                                     <  1:657595/143‑1 (MQ=255)           aGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGAc                                                                                                                  >  1:611057/1‑143 (MQ=255)            gAACGGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGTAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACt                                                                                                                 <  2:514392/143‑1 (MQ=255)            gAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACt                                                                                                                 >  1:609027/1‑143 (MQ=255)            gAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCGTTCTTGATAACGGCCGTTGCCGGAAAACCTCTCTCGATGAcg                                                                                                                 >  2:496221/1‑142 (MQ=255)             aaCAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCCCTCGATGAATg                                                                                                                >  2:123130/1‑143 (MQ=255)                  tGAAGGTGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTAc                                                                                                           <  1:371122/143‑1 (MQ=255)                       gcgcACAATTACTTCTGTTACTGCACCAAAAAGCGATCGCAGCGGTGTGCCATAAAGTCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGAt                                                                                                      <  2:270059/143‑1 (MQ=255)                        ctcacgATGACATCCGGTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATc                                                                                                     <  2:497588/141‑1 (MQ=255)                         gcACTATTACATCCCTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAACCCCAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCa                                                                                                    <  1:148761/143‑1 (MQ=255)                           ataaTTACGTCCGTTGCTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATg                                                                                                  <  1:554546/141‑1 (MQ=255)                                   aTCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGAt                                                                                          >  1:187321/1‑143 (MQ=255)                                     ccGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGAGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGAGTGGTTACGCGATAATGAAATAGATGa                                                                                        >  1:474418/1‑143 (MQ=255)                                      cGTTACTGCACCAAAACGCTATCGCAGCGGTGTTTCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGaa                                                                                       <  2:367521/143‑1 (MQ=255)                                            tGCACCAAGAAGCCATCGCAGCGGTGTGCCATTAATGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATc                                                                                 <  2:2324/143‑1 (MQ=255)                                            tGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATc                                                                                 >  1:772902/1‑143 (MQ=255)                                            tGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGGTACAGTGCCTTTTTTGTGAGCTGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATc                                                                                 >  1:184316/1‑143 (MQ=255)                                               aCCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAACATCCTTTAGTTGACAGTGACAGTGCCGTTTTTGATACCGGCCGTTGCCAGAAAACCGCGCTCGATGACTGGTTGCGCGATAATGGGATCGCTGAATTGGTCgag                                                                              >  2:373007/1‑141 (MQ=255)                                                ccAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTa                                                                             >  1:179685/1‑143 (MQ=255)                                                 cAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGCAAACCTCTCTCGATGACTGGTTACGCGATCATACAATCGATGAATTGATCGCTac                                                                            >  2:372152/1‑142 (MQ=255)                                                  aaaaaGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATg                                                                           >  1:54209/1‑143 (MQ=255)                                                  aaaaaGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAAACGATGAATTGACCGTCATg                                                                           >  2:668460/1‑143 (MQ=255)                                                          aTCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGc                                                                   <  1:723709/143‑1 (MQ=255)                                                                gCGGTGTTCGATGAAGGCGCAAGTCCGTTAGTGGTCAGTTACAGTGCCGTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGa                                                             <  2:199331/143‑1 (MQ=255)                                                                                   aataTCCTTCATTTGTCATTTACTGTGCCTTTTTGTATAACGGCCGTTGCCATAATACCTCTCTCGATGACTTGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTAcc                                          <  2:33432/140‑1 (MQ=255)                                                                                   aaaaTCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGGCTGGCTACTGTATATTTCGTGGAACTTAcc                                          >  2:416153/1‑143 (MQ=255)                                                                                    aaaTCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCg                                         <  1:17602/143‑1 (MQ=255)                                                                                                   cAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGAATATTGCGTGAAGTTTACCGTGCTGGACGCGTTAc                          <  1:115750/143‑1 (MQ=255)                                                                                                    aGTTACAGTGTCGTTTTTGATAGCGGCCGTGGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAGATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACa                         <  2:609027/143‑1 (MQ=255)                                                                                                    aGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACa                         <  1:158288/143‑1 (MQ=255)                                                                                                      ttACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGt                       >  1:766892/1‑143 (MQ=255)                                                                                                                        tAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGAGCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGGTTACCGTGCTGGACGCGCTACAGTTAGGTTAGAAGGtcgacg     >  1:239616/1‑138 (MQ=255)                                                                                                                         aaCGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTGCTAAGGTACACGt    >  1:283747/1‑143 (MQ=255)                                                                                                                           cGGCCTTTGCCAGAAAACTTCGATCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGa  <  2:533778/143‑1 (MQ=255)                                                                                                                           cGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGa  >  1:737944/1‑143 (MQ=255)                                                                                                                                  |                                                                                                                                        CACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGA  >  NC_000913/1852142‑1852406

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