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%	90.5 / ‑6.0	28	W69* (TGG→TAG)	xylA	xylose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/13);  minor base A (0/1);  total (14/14)
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.97e-01

TTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGC‑CTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTGATGG  >  NC_000913/3730338‑3730802                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   ttGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGC‑CTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCaa                                                                                                                                                                                                                         >  2:91310/1‑250 (MQ=255)  tGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGC‑CTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCaaa                                                                                                                                                                                                                        <  1:55730/250‑1 (MQ=255)      cGTTCCACTCAGCAGCTCCACGCCGCTCTAGTCTTGGTTGCACT‑CCATGACATCAAGCATCGGCGCAAAATTATTGATGTAGTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTAGAAGCAATAAACTTGCTGATGGAACTTGTGTAAAAACCGAAAGGTGACCTGTGCTTGACGCTTCGCCATCGCTAGTTCATCACCAGGCTGCTGCAACGGACGATTaaccgtcctcatatcgatcatag                                                                                                                                                                                                                   <  2:157647/251‑22 (MQ=255)             cacaGCAGCTTCACGCCGCTCTCTTCTTGCTTGC‑CTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGcccc                                                                                                                                                                                                            <  1:142761/251‑1 (MQ=255)                  cagcTTCACGCCGCTCTCTTCTTGCTTGC‑CTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTcc                                                                                                                                                                                                       <  1:236545/251‑1 (MQ=255)                                                  gccAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCaa                                                                                                                                                                         >  1:157177/1‑250 (MQ=255)                                                  gccAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCa                                                                                                                                                                          <  2:48186/249‑1 (MQ=255)                                                       gccagcAACCATTTGCGCAAAATGATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTGTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCACCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCa                                                                                                                                                                   <  2:183917/246‑1 (MQ=255)                                                               ccATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACACTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGctct                                                                                                                                                           <  2:166741/251‑1 (MQ=255)                                                                        cAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGc                                                                                                                                                  >  2:238165/1‑251 (MQ=255)                                                                         aaaaTTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACTCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATTCGCt                                                                                                                                                 <  1:132990/251‑1 (MQ=255)                                                                                      tACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACAcc                                                                                                                                      >  1:172655/1‑249 (MQ=255)                                                                                        cgtgtttAACGAAGCGCCGTCAGGGGAAACGTCCACATTGTGGAAGCAATAAAAGGGCCCATGTAACTGGTGGAAAAACTCAAATGCGACATCTGCTGTACGATTCGTCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCAGACGCTTACCCAACACCAGtt                                                                                                                                  <  2:232180/247‑1 (MQ=255)                                                                                            tttAACGACGCGCCCTCAGGGGCAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTtcgtc                                                                                                                              >  2:155054/1‑251 (MQ=255)                                                                                               aacgaagagacctccggcaaaaCATCCACATCGTCGAACAAATAAAATGGCACCTGTAACTTGTGGAAACACTCAAATGCGACCTCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGCGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCggg                                                                                                                           >  2:97118/20‑251 (MQ=255)                                                                                                             aGGGGAACCATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACgg                                                                                                              >  2:154965/1‑250 (MQ=255)                                                                                                                ggAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCAAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGAGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATg                                                                                                          >  2:209504/1‑251 (MQ=255)                                                                                                                  aaaCATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCt                                                                                                        >  2:53480/1‑251 (MQ=255)                                                                                                                   aaCATCCACATCGTGGAAGCAAGAAAATGGCCCATGTAACGTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAACCGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCATTACGTCGGTATTTTAGTGACGGAATGCTa                                                                                                       <  2:240816/251‑1 (MQ=255)                                                                                                                        ccACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACgg                                                                                                   <  1:187001/250‑1 (MQ=255)                                                                                                                                      cAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGCTTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTg                                                                                     >  2:80343/1‑250 (MQ=255)                                                                                                                                                         aCTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAAcgcg                                                                  >  2:173903/1‑250 (MQ=255)                                                                                                                                                                  aaaactcaaccccGACATCTGCTTTCCGCTTCGACAACGCCAGTGCCTCACCAGGCTGCTGCTACGGCAGATTCAACGCCCCCACACCACACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCTGGCGGCAAAACGCCAGTGCTCTTCCATACGCTTACCCAACACCAGGTCGTCGAGATTGTAGGGCCGGAATGCTAACGGCGTTGAGGCTTTTGAGCCTTCGTAACGAACGCGAACCAGCTgg                                                        >  2:39147/13‑251 (MQ=255)                                                                                                                                                                               gACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCtt                                           <  1:224219/251‑1 (MQ=255)                                                                                                                                                                                                   ccAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATc                       >  1:101511/1‑251 (MQ=255)                                                                                                                                                                                                     aaCGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAg                     <  1:45360/251‑1 (MQ=255)                                                                                                                                                                                                               ccTCACCAGGCGGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCcgc           <  2:64143/251‑1 (MQ=255)                                                                                                                                                                                                                         ctgctgCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCACCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTgatgg  >  1:153173/1‑250 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   TTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGC‑CTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTGATGG  >  NC_000913/3730338‑3730802

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