Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,534,321	C→T	P100S (CCA→TCA)	nhoA →	arylamine N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,534,321	0	C	T	96.3%	69.0 / ‑5.0	27	P100S (CCA→TCA)	nhoA	arylamine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (6/20);  total (6/21)
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.33e-01

CGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCC‑GTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGTACCGTTTTGATCTCTGCGAGC  >  NC_000913/1534094‑1534558                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               cGCTTCGTGCATTGCACATTAAACACAATTGCACCATTCC‑GTTTGAAACCCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGGTTACTGTTTGGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCAt                                                                                                                                                                                                                         <  2:288698/251‑1 (MQ=255)                       acacAATTGCACCATTCC‑GTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTgg                                                                                                                                                                                                   <  1:195467/251‑1 (MQ=255)                         acaagttcACCATTCCTGTTT‑AAATCCTCGACTGTTTGCTGCCGAGGGAAATACATCTTGATAATCAAGCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTGACGCGAGCTGGGGGTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTTCTGTTGGTGGAACTGGaa                                                                                                                                                                                                 <  2:370127/244‑1 (MQ=255)                                              gAAAACCTCGCCGGTTTGCTGCCGCGTGAAATACCACTTGAGAATCAAGCGCCGGAAGAGAATCTGTTGATAGCCCGTCGTGTCGGTGACTGTTTTGATCGGAATGGCGTGTTTGAGCGGGTGGGACGAGAGCTGGTGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTGAGCATTACCGCCGCGCATCCATCGGTTGCTGTTGGTGTAACTGGAAGATGAAAAATGGAGTGCTg                                                                                                                                                                              <  2:395454/250‑1 (MQ=255)                                                  aCCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑c                                                                                                                                                                        <  1:386316/251‑1 (MQ=255)                                                       gACGTTTGGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGttt                                                                                                                                                                   <  1:361419/251‑1 (MQ=255)                                                                                        aatcaatcGCCGGAAGAGAACCTGGTTATAGCCCGTCGTGGCGGTTTCTTTTTTGAGCAGAAGGGCGTTTTTTAGCGGGTGTGACTCGAGCTTGGGTTTAACGTTCGCAGCTTGTTAGGGCGTGTGGTGTTATCAAATCCGTCAGCAGTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGttt                                                                                                                                  <  2:456800/251‑1 (MQ=255)                                                                                                  cGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTc                                                                                                                                      <  1:461172/237‑1 (MQ=255)                                                                                                          aaaCTGGTGATTGCCCGGCGTGCCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGATCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTGCCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATTGATTGCTGATGT‑CGGTTTCGGGGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCa                                                                                                                <  2:520154/251‑1 (MQ=255)                                                                                                                         cgtcgtGGCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACgg                                                                                                 <  1:68433/251‑1 (MQ=255)                                                                                                                                       cTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGgttgt                                                                                   >  1:93391/1‑251 (MQ=255)                                                                                                                                              gAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGcc                                                                                                                                           <  1:511997/188‑1 (MQ=255)                                                                                                                                                       ggCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGtgatg                                                                    <  1:17521/250‑1 (MQ=255)                                                                                                                                                                            cgctggCTGGGGTTTAACGGTCGCGGCTTGTTAGGGCGTGTTGTGTTCTCAAGTCCGTCAGCATTCCCGCCGCGCACCCTGCGTTTGCGGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCATACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGTAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGCGTGTTGCAGTGTAatca                                              <  2:518807/246‑1 (MQ=255)                                                                                                                                                                            cgcgAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAatc                                               >  1:59072/1‑250 (MQ=255)                                                                                                                                                                              ggAGCGGGGGTTTCACGTTCGCAGCTTTTTAGGTGGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTGCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCATGAGGGTGAGGATTGGGTGTTGCATTTTAatcatc                                            <  2:516668/250‑1 (MQ=255)                                                                                                                                                                                   tGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAAtcatcatca                                        <  2:222297/250‑1 (MQ=255)                                                                                                                                                                                   tGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAAtcatcatca                                        >  1:222297/1‑250 (MQ=255)                                                                                                                                                                                   tGGGGTTTAACGTTCGCAGATTGTTAGGGCGCGTAGTGTTATAAAATCAGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAGTGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGATAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGAAGCTTGGGTGTTGCAGTTTCatcttgatc                                         >  2:248409/1‑243 (MQ=255)                                                                                                                                                                                     gggTTTAACGTTGGCGGCTTGTTAGGGCGGGTAGTGTTAACAAATCCGTCAGCATTGCCGCCGCGCACCCATCGTTTGCTGGTGGTGGAACTGGAAGAGGAAAAATGGATT‑CGGGGGTTCGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTGTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTGAATCATGAgcagca                                     <  2:436742/251‑1 (MQ=255)                                                                                                                                                                                       gTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATGACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCAtt                                   <  1:248409/251‑1 (MQ=255)                                                                                                                                                                                             cGTTCGCCGCTTTTTAGGGCGCGTATTGTTATCAAATCCGTCAGCATTACCGTCGCGCACCCATCGTTTGCTGTGTGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTGCGGTGGTCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACTGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGt                             <  2:572164/251‑1 (MQ=255)                                                                                                                                                                                                    aGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGt                       >  1:469239/1‑250 (MQ=255)                                                                                                                                                                                                                  cgTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGTACCGTTTTGATctct        <  1:497885/251‑1 (MQ=255)                                                                                                                                                                                                                   gTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTg                                                         <  2:175537/201‑1 (MQ=255)                                                                                                                                                                                                                   gTAGTGTTATCAAATCCGTCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTg                                                         >  1:175537/1‑201 (MQ=255)                                                                                                                                                                                                                         gtattaaaTCCGCCCGCATAAACTCCGAGCACCCATGGTGTGCTGTGGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGG‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGTACCGTTTTGATCTCTGCGagc  <  2:281860/245‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               CGCTTCGTGCATTGCACCTGAAACACAATTGCACCATTCC‑GTTTGAAAACCTCGACGTTTTGCTGCCGAGGGAAATACAGCTTGATAATCAATCGCCGGAAGAGAAACTGGTGATAGCCCGTCGTGGCGGTTACTGTTTTGAGCAGAATGGCGTGTTTGAGCGGGTGTTACGCGAGCTGGGGTTTAACGTTCGCAGCTTGTTAGGGCGCGTAGTGTTATCAAATCCGCCAGCATTACCGCCGCGCACCCATCGTTTGCTGTTGGTGGAACTGGAAGAGGAAAAATGGATTGCTGATGT‑CGGTTTCGGTGGGCAGACGCTAACCGCGCCGATTCGTTTAGTTTCCGATCTCGTGCAGACCACGCCACACGGAGAGTATCGGTTGTTGCAGGAGGGTGATGATTGGGTGTTGCAGTTTAATCATCATCAGCATTGGCAGTCGATGTACCGTTTTGATCTCTGCGAGC  >  NC_000913/1534094‑1534558

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