Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	64,713	T→C	A356A (GCA→GCG)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	64,713	0	T	C	96.7%	91.6 / ‑6.1	30	A356A (GCA→GCG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (12/17);  minor base A (0/1);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTC  >  NC_000913/64494‑64923                                                                                                                                                                                                                            |                                                                                                                                                                                                                    gCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCa                                                                                                                                                                                       <  1:7183/250‑1 (MQ=255)           taCGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCTAGGAAAAATGGCATGATTTTAGTTTTATggt                                                                                                                                                                             >  2:35623/2‑249 (MQ=255)                          ttATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCa                                                                                                                                                              <  1:254165/250‑1 (MQ=255)                             taGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGGAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTc                                                                                                                                                          >  2:90426/1‑251 (MQ=255)                                 ccagcaccagtaCCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCCCTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCa                                                                                                                                                       <  2:212021/250‑1 (MQ=255)                                  cagcacaagtacCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCAGACAGGGCTGACGTGTGGCGGCACTTCGCCGAGCTTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGCTGACCAAACGCCGCCACCGATCACCCGTGGCGGTATACCGGCAGGCATGTCACAGTCGCTCGTTCGCGTTAATATGGCATGTTTTTAGTTTTTTGGAAAATCCACGTCCCCAGcgcacaa                                                                                                                                                      >  1:249418/1‑244 (MQ=255)                                     caccagtacCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCtt                                                                                                                                                  <  1:39194/251‑1 (MQ=255)                                        cagtaccgAATCATAAAGCCCTGGCCGTGAATCCAGCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCa                                                                                                                                               <  1:26429/244‑1 (MQ=255)                                             ccgAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCTCCAGGGCTGTCGTGTGTCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGATAATACCAAACGCCGCCACTGATCCCCCTTGGCTGTCTACCGTCAGGCCGTTCACCGTCGCTCTTTCGAGTAAAAATGGCATGTTTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTATCAGGTta                                                                                                                                          <  1:172065/249‑1 (MQ=255)                                                  tCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCAACAGCTCGCAATCTTTCAGGTTATAAGt                                                                                                                                     >  2:231923/1‑251 (MQ=255)                                                   aagaaagaccgggCCGTGAATCCGTCACGTAGCCGCCAGAGCTGGCGTGTGGAGGCACTTCGCCGAGATTAGGCGCGACATAGCCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCAACGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCAACGCCGCTCGTTCGAGTAAAAAGGGCATGATTGCAGTTTTATGGAAAAGCTGCGTCACCAGCTCGCAATCCTTCAGGTTATAAGt                                                                                                                                     <  2:73763/239‑1 (MQ=255)                                                        aGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAg                                                                                                                               <  1:129805/251‑1 (MQ=255)                                                                                   ccgccAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGATCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAAcggcg                                                                                                    >  2:120790/1‑251 (MQ=255)                                                                                            atGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCCGTCACCGTCGCGCGTTCGAGGAAAAAAGGCATGATTTCAGGTTTATGGGAAATCTGCTTCACCAGCTCGCAATCTTTCAGGTTCTAAGTCGCCAGCGCCGGTGTATCTTCGGAGAAACGGCGGTACAAttt                                                                                           >  2:208897/2‑251 (MQ=255)                                                                                                      gcggcACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGa                                                                                <  1:231923/251‑1 (MQ=255)                                                                                                                     gaTTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCg                                                                 <  1:203293/251‑1 (MQ=255)                                                                                                                         aGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAg                                                             <  1:246336/251‑1 (MQ=255)                                                                                                                             gcgACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAAATTCGGCCATTCGATCCCACGGGGTTTCGATaacatt                                                         >  2:220305/1‑246 (MQ=255)                                                                                                                               gACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTc                                                       <  2:52046/251‑1 (MQ=255)                                                                                                                                aCATAACCAGCGCGGTGCATACGAGGGAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCGCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTGTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTGCGTCCATTCGATCCCACGGGTTATCGATAGATTTTcc                                                      <  2:94827/251‑1 (MQ=255)                                                                                                                                aCATAACAAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGACAATCTGCGTCACCAGCTCGCACTCTTTCAGGTTATCAGTAGCCCGCGCAGGTTTATCTTCGGCGAAACGGCGGG‑CCATTTCGTCCATTTTATCCCACGGGGTATCTATAGATTTTcc                                                      >  2:53938/1‑251 (MQ=255)                                                                                                                                   tAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCCTTCGATCCCACGGGTTATCGATAGATTTTCCTTc                                                   >  1:151641/1‑251 (MQ=255)                                                                                                                                          cgcgGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAAt                                            <  1:214615/251‑1 (MQ=255)                                                                                                                                          cgcgGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGGTATCGATAGATTTTCCTTCTCCTAAt                                            >  1:82465/1‑251 (MQ=255)                                                                                                                                                      gaggAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑TAATTTCGTCCATTCGATCCCACGGGTTATCGATAGAATTTCCTTCGCCTACTAACTCCTGAGCg                                >  2:90344/3‑251 (MQ=255)                                                                                                                                                             aTAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGAcgggtt                         >  1:226213/1‑246 (MQ=255)                                                                                                                                                               agagATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTc                        <  1:148857/250‑1 (MQ=255)                                                                                                                                                                gagaTGACCAAACGCCGACACCGATTCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTAGATCGGTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCGTCGGCGAAACGGCGGT‑CGATTGCGTCCATTCGATCCCACGTGTGATCGATAGATTTTCCTTCGCCTAATAACGCCTGAGCgcagggttcta                      <  2:178985/251‑11 (MQ=255)                                                                                                                                                                                   aCCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAAtt   <  1:166707/251‑1 (MQ=255)                                                                                                                                                                                    ccgatcccccgTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTc  <  1:70926/243‑1 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                    GCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTC  >  NC_000913/64494‑64923

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