Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,730,560	C→T	W69* (TGG→TAG)	xylA ←	xylose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,730,560	0	C	T	96.4%	79.7 / ‑4.4	28	W69* (TGG→TAG)	xylA	xylose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (12/15);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGG  >  NC_000913/3730316‑3730796                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             gTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCGCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGAc                                                                                                                                                                                                                                        <  1:165470/251‑1 (MQ=255)         ggTTTGTAAATCAGTTGGCCGTTCCCCACAGCAGCTGCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCAGTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAAGGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAAc                                                                                                                                                                                                                                <  1:287018/251‑1 (MQ=255)                aacgCAGTTGGCCGTTCCCCACATCAGCTGCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCGCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCca                                                                                                                                                                                                                         <  2:55793/248‑1 (MQ=255)                   gCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCacac                                                                                                                                                                                                                      <  1:60831/251‑1 (MQ=255)                       ttGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGCACATCCACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCCCATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTCCGCTTCGCCAACGCAAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCccaccac                                                                                                                                                                                                                   >  2:53578/1‑249 (MQ=255)                            cGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACatat                                                                                                                                                                                                             <  1:290804/251‑1 (MQ=255)                                         agcTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTACACGCCCCCACACCAAACATATCCGCCCCGCTccc                                                                                                                                                                                                >  2:321334/1‑250 (MQ=255)                                                     ctctctTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGtgtg                                                                                                                                                                                    >  1:371328/1‑251 (MQ=255)                                                                   gcctgccAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCGTGTAACGACGCGCCCTCAGGGGAAATATCCACATCGTGGAAGCAATAAAATGGCACATGTATCTTGTGTAAAAGCTCAAATGCTACATCTTCTTTACACTTCGCCAACGCCAGTGCCTCACCATGCTTCTGCTACGGACGATTAAACGCCCCCTCACCAAACATTGCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAggcg                                                                                                                                                                       <  2:331089/250‑1 (MQ=255)                                                                        ccAGGACATCCACCATTTGCGCAAAATTATTAATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCATTAAAATGGAACATGTACCTTATGGAAAAACTCAAATGCGAAATCTGCTTTACGCTTCGCAAACTCCATTGCCTCACCAGGATGCTGCTACAGACGATTAAACGCCCCCACAAAAAACATATCCGCCCCGTAACAGCAGAAGGTGTGCCAGTAGCAAGCGACaaa                                                                                                                                                                  >  2:154500/1‑250 (MQ=255)                                                                                                aaTTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCCATAAAATGGCACATGTAACTTGTGGAACAAATCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCCCACCAAACATATCCGACCAGTTCCAGCCGTAGGAGTGCCAGTAGCCGGCGGCACAACGCAAGAGCTCTTCCATAAGCtt                                                                                                                                          >  2:281117/1‑250 (MQ=255)                                                                                                        aTGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAacac                                                                                                                                 <  1:41617/251‑1 (MQ=255)                                                                                                                       gacgCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGAtt                                                                                                                  <  1:321334/251‑1 (MQ=255)                                                                                                                            gcCCTCAGGGGAAACATCAACACCGTAGAAGCAATAAAATAGCACATGTAACTTGTGGAAAAACTCAAATGCGAAATCTGATTTACGCTTCCCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCACACACCAAACCTATCCGAACCGTTCCAGAAGAAGTTGTACCCGAAACCGGCGGCAAAACGCAAGGGCTCTTCCATACGCTTACCCACCACCCGTTAGTCGGGAtttaagg                                                                                                             >  2:238621/1‑246 (MQ=255)                                                                                                                                  aGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCaa                                                                                                                                                            >  1:176319/1‑198 (MQ=255)                                                                                                                                   ggggAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTAAAATGCGACATCTGCTTTACGCTTAGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATGGTAGTGACGGa                                                                                                       >  1:51131/1‑250 (MQ=255)                                                                                                                                                caTCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGttt                                                                                         <  1:382465/251‑1 (MQ=255)                                                                                                                                                   cGTGTAAGCAATAATAGGGTACAGCTAACTTGTGGAAAAACTTATAGGCGACTTCGTATTTACGGTTCGCCTACGCCAGTGCCTCACCAGGCTGCTTCCACGGTCGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAg                                                                                      <  2:380602/251‑1 (MQ=255)                                                                                                                                                         aGCAATAAAATGTCAAATGTAACTGGTGGAAAAACATAAATGTGACATATGCTGTACGTTACGCAAACGCCATGGCATCACCAGGCTGCTGCTACGGACGTTTAAACGCCCCCACACCAAACTTAGCCGCCCCGGTCCAGCAGAAGGTGTTCCAGTAGCAGGCGGCAAAACGCCAGTGCCCGTCCATACGCTTACACAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGAtttt                                                                                <  2:151527/251‑1 (MQ=255)                                                                                                                                                                        caTGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGATGCTGCTACGGACGATTAAACGCCCCCACACCCAACATATCCGCCCCGTTCCAGCAGAAGTTGTGCCCGTAGCAGGCAGAAAAAAGCCAGTGCTCTTCCATACCCTTACCCAACACCAGTTCGTCGAGCTTTTAGTGACGGTATGCTAACGGGTTTGGGGATTTTGATCCTTccttacgc                                                                 >  2:342602/1‑244 (MQ=255)                                                                                                                                                                         aTGTAACTGGTGGAAAAACTCAAATGCGACACCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATaacgaa                                                                <  1:268689/251‑1 (MQ=255)                                                                                                                                                                                      aaaaaCTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTg                                                   <  1:152645/251‑1 (MQ=255)                                                                                                                                                                                         aaCTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTc                                                <  1:83284/251‑1 (MQ=255)                                                                                                                                                                                                              ttACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGa                            >  1:212604/1‑250 (MQ=255)                                                                                                                                                                                                                  gCTTCGTCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTcc                       <  1:146504/251‑1 (MQ=255)                                                                                                                                                                                                                       gCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAACCGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAAt                  <  1:37867/251‑1 (MQ=255)                                                                                                                                                                                                                                aGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTAGACGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGAATATTGAACTCCATAATCAGGTagtg         >  1:72467/1‑248 (MQ=255)                                                                                                                                                                                                                                       cACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCACTAATCAGGAAATGCCGCggg  >  2:190428/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             GTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGG  >  NC_000913/3730316‑3730796

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