Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	147,213	2 bp→TT	coding (246‑247/903 nt)	rpnC →	recombination‑promoting nuclease RpnC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	147,213	0	C	T	96.3%	80.5 / ‑3.6	27	H82H (CAC→CAT)	rpnC	recombination‑promoting nuclease RpnC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (12/14);  total (12/15)
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.94e-01
*	NC_000913	147,214	0	C	T	96.2%	83.3 / ‑3.5	27	Q83* (CAA→TAA)	rpnC	recombination‑promoting nuclease RpnC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (12/13);  total (12/15)
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.95e-01

GGATGCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCAT  >  NC_000913/146970‑147446                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                         ggatgCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagca                                                                                                                                                                                                                                    <  1:161731/251‑1 (MQ=255)         cGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGAt                                                                                                                                                                                                                            >  1:141840/1‑251 (MQ=255)                  cacCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATCCATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACAAAGCCCGGACGTGCTCTATTACGTGCAAATGCCGGGCAATCCCGGTTATCTGCATGTTGTGGTTGAACATTAAAGCAAGCCTGATAAGaaaagg                                                                                                                                                                                                                   >  2:351570/1‑249 (MQ=255)                       cATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCtt                                                                                                                                                                                                              <  1:63474/251‑1 (MQ=255)                          gACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATGAAAGCAAGCCGGATAAGAAAATGGCCGTTc                                                                                                                                                                                                            <  2:158450/250‑1 (MQ=255)                            cgcgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTAGCAAGCGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGCTTCCTTGAAGAGGGCCTGAAAGGAAACCGCAAGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGcc                                                                                                                                                                                                         >  1:350792/1‑250 (MQ=255)                             gcgGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAATCAAGCCGGATAAGAAAATGGCCTTTCGCat                                                                                                                                                                                                        <  1:334313/251‑1 (MQ=255)                                                                 gCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCt                                                                                                                                                                    <  2:251331/251‑1 (MQ=255)                                                                    cgcgACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGa                                                                                                                                                                 <  1:305531/251‑1 (MQ=255)                                                                                  gaTACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGAAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTATATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGAt                                                                                                                                                    >  2:164183/1‑250 (MQ=255)                                                                                             cAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCtg                                                                                                                                        >  2:189465/1‑251 (MQ=255)                                                                                                         gcgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATa                                                                                                                            <  1:259849/251‑1 (MQ=255)                                                                                                                         cAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGc                                                                                                             <  1:189465/250‑1 (MQ=255)                                                                                                                               gCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAAGCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCca                                                                                                      <  2:46469/251‑1 (MQ=255)                                                                                                                                  tCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCcac                                                                                                     <  1:258362/249‑1 (MQ=255)                                                                                                                                   cATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACAcc                                                                                                  <  2:141840/251‑1 (MQ=255)                                                                                                                                    aTTTAGAGTCGATGAGTTTAATTGAGAAGAGCCTGAAAGGACACAGCCCGAGCGTGCTCTATTCCGTGCAAATGCCGGGCAATCCCGGATCTCTGCATGTTGTGATTTAAAATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCAACGGCCTATGTGGGCTGACCCCGATAAGCTGCCGCTGGTGATGCCGATACTGTTTTATCCAGGCGAAGCCACACCt                                                                                                 >  2:94465/1‑251 (MQ=255)                                                                                                                                                        aTTGAAGAGAGCCTGAAAGGACACAGCACGAACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAAAATTAAAGCACGCCGGATAAGAAAATGGCCTTTCGCCTGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATACGCTATCACtgtgatg                                                                             >  1:351232/1‑248 (MQ=255)                                                                                                                                                              gagagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCACTGTGCTGGTTTgg                                                                       >  2:72743/1‑250 (MQ=255)                                                                                                                                                                           gACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCg                                                          >  1:293540/1‑251 (MQ=255)                                                                                                                                                                              acaGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATCAGAAACTGGCCTTTCGCATGATGCGTTAATCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGCTACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCg                                                       >  2:110668/1‑251 (MQ=255)                                                                                                                                                                                 gCACGGACGTGCTCTATTCCGTGCAAATGAAGGGCAATCCCGGTTCTCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCgg                                                      >  2:242123/1‑249 (MQ=255)                                                                                                                                                                                      gACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGgc                                               <  1:227613/251‑1 (MQ=255)                                                                                                                                                                                                          aaTGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTAGCAGGGCGAGGCCACACCTTATCCGCTATCATTGTGCTGGTGTGATAGGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTc                           <  1:205051/251‑1 (MQ=255)                                                                                                                                                                                                             gCAGGGCAATCCCGGTTGTCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTATATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCtt                        <  2:194397/251‑1 (MQ=255)                                                                                                                                                                                                              cAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCttt                       <  1:176549/251‑1 (MQ=255)                                                                                                                                                                                                                                   cATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCAt  >  1:84550/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                         GGATGCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCAT  >  NC_000913/146970‑147446

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