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.8%	88.0 / NA	33	W69* (TGG→TAG)	xylA	xylose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (18/12);  minor base G (0/1);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTGAT  >  NC_000913/3730314‑3730800                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 ccGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACgg                                                                                                                                                                                                                                              >  1:103534/1‑251 (MQ=255)                               gTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTGGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATc                                                                                                                                                                                                                <  2:175807/251‑1 (MQ=255)                                  ccccACAGCAGATTCACGCCGCTCTCTTCTTACTTGCCTGCCAGGCCATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGc                                                                                                                                                                                                             >  1:102628/1‑251 (MQ=255)                                    ccACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCAAATCGTGGAAGCAATAAAATGGCACATGTAACTTGCGGAAAAACTCAAATGCGACATAGGCTTTCCGCTTCGCCAACGCCAGTGCCTAACCAGGCTGCTGCTCCGGACGCTTAAACGCCCCCACACCAAACATATCCGccc                                                                                                                                                                                                           >  2:24500/1‑251 (MQ=255)                                          gagcTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGAAGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAAAGCCCCCACACCACAAATATCCGCACCGtacc                                                                                                                                                                                                     >  2:245571/2‑248 (MQ=255)                                           ggcTTCACGCCGCACGCTTCTTGCTTGCCTGCCAGGACATCACCCATTTGCGCAATATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATAGTGGAAGCAATAAAAAGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTAAGCTTCGCAAACGCCAGTGCCTAACACGGCTGCTGCTACGGACGATTAAACGCACCCACACCAAACTGATCCGCCCAGTTca                                                                                                                                                                                                     >  2:10959/2‑249 (MQ=255)                                                         ctctTCTGGCTTGCCTGCCAGGACATCAACCAGTTGCGCCAAAGTATTGATGTACTCTTTTAACTACGCTGCCTCAGTTGAAACATCCTCATCGTGTAAGCAATAAAATGGCACATGTAACTTGTGTAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCTCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGtgtg                                                                                                                                                                                        <  2:102628/249‑1 (MQ=255)                                                            ttcttGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCacac                                                                                                                                                                                                                          >  2:28766/1‑212 (MQ=255)                                                            ttcttGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCacac                                                                                                                                                                                                                          <  1:28766/212‑1 (MQ=255)                                                                 gcttgcCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGc                                                                                                                                                                                >  1:152625/1‑249 (MQ=255)                                                                 gcttgcCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATAGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCAAATGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGAAGc                                                                                                                                                                                >  2:61125/1‑249 (MQ=255)                                                                              gACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTAACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGCTTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGc                                                                                                                                                                  >  1:184995/1‑250 (MQ=255)                                                                                   cAACCATTTGCGCAAAATTAATGATGTACTCATGTCAAGACGCGACCTCAGGGGAAACAACCACAGAGTGGACGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCACAGGCGACATCTGCATTTCGCTTAGCCAACGCCAGAGCCTCACCAGGCGGCTGCAACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGCGCCAGGAGCAGGCGGCAAAACGCAAGTGc                                                                                                                                                            <  2:11357/251‑1 (MQ=255)                                                                                             cgcAAAATTATTGATGTACTCTTGTAACGACGCGCCCTCAGGGGAAACATCCTCTACGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTTCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTTCTACGGACGATTAAACGCCCTCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATAc                                                                                                                                                  <  1:272225/251‑1 (MQ=255)                                                                                                     tattGATGTACTCATTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTAcc                                                                                                                                           <  2:103534/250‑1 (MQ=255)                                                                                                          aTGTACTCTTTTAACGACGCGCCCGCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAacac                                                                                                                                     <  1:275438/251‑1 (MQ=255)                                                                                                               ctctTTTAACGAAGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGtt                                                                                                                                >  1:64095/1‑251 (MQ=255)                                                                                                                    ttAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTtcgtcg                                                                                                                           >  1:153470/1‑251 (MQ=255)                                                                                                                      aaCGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCggg                                                                                                                         <  1:61516/251‑1 (MQ=255)                                                                                                                           cgcgCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTCCCCAACACCAGTTCGTCGGGATTGt                                                                                                                    >  1:138901/1‑251 (MQ=255)                                                                                                                              gcCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGt                                                                                                                 >  1:130356/1‑251 (MQ=255)                                                                                                                                     ggggAACCATCCACGTCGTGGAAGCACTAAAATGGCAAATGTAACTTGTGGAAAAACTCAAATGCGACATCTGATTTACGCTTCGACAACGCCAGAGCCTAACCAGACTGCTGCTACGGACGATTAAACGCCCCCACACCAAACCAATCCGCCCCGTTCCAGCAGAACGTGTGCCAGTAGCAGGCGGCAAAACGCAACGGCTCTCCCATACGCTGACCCAACACCAGATCGTAGGAATTGtggggaaggc                                                                                                           >  2:33946/1‑241 (MQ=255)                                                                                                                                                   atagTGGAATCAATAAAATGGACCATGTAACATGTGGAAAAAATCAAATGCGACATCTGCTTTACGTTTCGCCAACGACAGTGCCTCACCAGGCTGCGGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACTCCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGCTTg                                                                                            <  2:170870/248‑1 (MQ=255)                                                                                                                                                                       gCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGCGCCTTCATa                                                                         >  1:286478/1‑250 (MQ=255)                                                                                                                                                                              tAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAAcgc                                                                 >  1:90466/1‑251 (MQ=255)                                                                                                                                                                                  ttGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATc                                                             >  1:73229/1‑251 (MQ=255)                                                                                                                                                                                         aaaaCTCAAATGCGACATCTGCTTTACTCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGGTCCAGCAGAAGGTGTGCCAGTAGCAGGAGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTgg                                                      <  1:16348/251‑1 (MQ=255)                                                                                                                                                                                         aaaaCTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTgg                                                      <  1:194639/251‑1 (MQ=255)                                                                                                                                                                                                             gccttctgCTTCTCCAACGCCAGTGCCCCCCTAGGCTTCTCCGCCGGACGATTAATCGCCCCCACTCCACACAAATCCGCCCCGTTCCAGCTGAAGGTGTGCCAGTAGCAAGCGGCAAAACGCAAGTTCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTATGATGATTTTGAGCCTTCATAACAAACGCGATCGAGCTGGTTAAAATAGGCTTTCatat                                   <  2:43067/243‑1 (MQ=255)                                                                                                                                                                                                                 tACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCAAGCTGGTCAAAATAGGCTTGCATATTGAAc                              >  2:103172/1‑251 (MQ=255)                                                                                                                                                                                                                              cGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGAGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAggg                 >  2:264512/1‑250 (MQ=255)                                                                                                                                                                                                                                     gCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCcg          <  1:139779/251‑1 (MQ=255)                                                                                                                                                                                                                                             aGGCTGCTTCTACGGCCGAGTAAACGCCCCCACACCAAACAGAACCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGGTCGTCGGGATTGTAGTGACGGAATGCTAACTTGGTTGAGGATTTGGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTgat  <  2:286478/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 CCGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTGAT  >  NC_000913/3730314‑3730800

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