Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,118	1	.	C	73.1%	43.1 / 14.7	26	coding (1254/1455 nt)	xylB	xylulokinase
Reads supporting (aligned to +/- strand):  ref base . (7/0);  new base C (10/9);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 5.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATATATAGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCC‑G‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTAC  >  NC_000913/3727868‑3728362                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       gaTATATAGATGTGAATTATCCCCCACCCGGTCCGGCAGGGGATAACGTTTACGCCATTAATGGCAGCCGTTGCTGATAGAGGCGCCGGCACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGCCTGTTCTAACGGTAGTTGAGGCAACAATTCCATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGCGCTGGCCCC‑CCA‑Ccccc‑t‑ccc                                                                                                                                                                                                                                                         >  2:106957/1‑251 (MQ=255)      gtaGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTAGCCGCGATCTGCGCCAGACTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑Tcccc‑c‑ccccgaac                                                                                                                                                                                                                                                    >  2:10219/2‑246 (MQ=255)      ataGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTAAGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACCATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑CCA‑Tcccc‑c‑ccccgtcc                                                                                                                                                                                                                                                    >  1:130167/1‑246 (MQ=255)       taGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTACTGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCCGTGCTGGCCCC‑ACA‑Taccc‑c‑ccccccacc                                                                                                                                                                                                                                                   >  1:284205/1‑245 (MQ=255)             tgAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATc                                                                                                                                                                                                                                           >  1:7802/1‑251 (MQ=255)                 ttATCCCCCACCCGGACAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCCGCCTTGCTGCGCCCCGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCC‑CCTACG‑GTaaccacgc                                                                                                                                                                                                                                       >  2:117240/1‑245 (MQ=255)                    tCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCCAC‑CCC‑C‑TACG‑GCCACCGAgctgc                                                                                                                                                                                                                                    >  1:217720/1‑251 (MQ=255)                         cACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGAcc                                                                                                                                                                                                                               <  1:44654/251‑1 (MQ=255)                              ggTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑CACG‑GTCAACGAGCCTCTGACCGccgc                                                                                                                                                                                                                          >  1:5100/1‑248 (MQ=255)                                                   tACGCCATTACTGGCAGACGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAAAGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCC‑C‑CACCCGGAATCCAGCTGCTCTCCGCTGATCTCCGCCAGCATCTGACGcc                                                                                                                                                                                                     >  2:140690/1‑251 (MQ=255)                                                        cATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTg                                                                                                                                                                                                          >  1:234573/1‑241 (MQ=255)                                                        cATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGGCGGGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTg                                                                                                                                                                                                          <  2:234573/241‑1 (MQ=255)                                                         aTTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACCATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACC‑TCCCC‑CACCC‑C‑TAAA‑GTACTCGCGCTGCTTACCGATGATAGCTCCCAGCCTCTGCCGCCAgaacc                                                                                                                                                                                               >  2:11546/1‑247 (MQ=255)                                                         aTTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATAGCCGCCAGCATCTGACGCCAGTACt                                                                                                                                                                                               <  1:117240/251‑1 (MQ=255)                                                                    aaGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGGGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCACAGTGCTGGCCCC‑ACAACCCCC‑C‑CCC‑G‑CACG‑GTAATCGATCTGCTGACCGCTGCTATCCGCCAGCACCAGACGCCAggcatccctactagcc                                                                                                                                                                                    >  2:228533/1‑236 (MQ=255)                                                                                                tctcGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGc                                                                                                                                                         <  1:312771/250‑1 (MQ=255)                                                                                                      cgtGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TCCG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCACCATCTGACGCCAGTACTAACTACGCCCCCCCCCCCCAATCaaaccaacaccttgcggataa                                                                                                                                                   >  1:76073/1‑231 (MQ=255)                                                                                                                    ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑GCCG‑GGAATCGGGCTGCGGACCGCCGATACCCGCCAGCCTCTGACGCCAGTAATCCTTACGCGCCCCGCCCCCAATCCACGTAACACTATACGGTTTAATACCGCAGGCacgc                                                                                                                                    >  2:167869/1‑248 (MQ=255)                                                                                                                    ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTATCTCGATTCGCAGCGATATGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCACG‑CACCC‑C‑CCC‑G‑TAAG‑GTGATCACGCCGCTGGCCGTTGCTGTTCGCCAGCATCTGAATCCAGTACTCCCTACGCACCCCGCCCCCCATCAACGTAAGACttgcaggaataaaaccccagacaacc                                                                                                                                    >  2:302723/1‑227 (MQ=255)                                                                                                                    ggCATAACGCTGAGCATCTGGTAGATGCGACTGTTCTAACGGTCGTTGCGGCAACAATTCCATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGATGCGGACAGCTGATATCCGCAAGCACCTGACGCCAGAACGCACCACGCGCCCCGCCCCCAATCCACGCAACACCTTGCGGCTTAATACCGCAGGccagc                                                                                                                                    >  2:25327/1‑247 (MQ=255)                                                                                                                           cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAGCAACGTAACACTTGGCGGTTTAATACCGCAGGCATGCACGACAt                                                                                                                             <  1:119188/251‑1 (MQ=255)                                                                                                                                          aGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCACTGCGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCATTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTCATCGAGCTGCTGCCCGCTGATATCCGCAAGCATCTGCCGCCAGTACTCACTACGCGCCCCGCCCCCACTCACCGTAACAATTTCCGCTTTAATACCGCAGGCATGCACGCCACCCATGccccctgcc                                                                                                               >  1:203662/1‑243 (MQ=255)                                                                                                                                            aggCGACTGTTCTAACGGTAGTTGCGGCAACATTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCTCCCCCAATCAACGTAACACTTGGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc                                                                                                            <  1:136334/249‑1 (MQ=255)                                                                                                                                                                              tCAATGAGCGATTTCTCTGGATTCGCCGCGATCCGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACAACCCCC‑C‑CCC‑C‑TACC‑GGAATCGAGATACTGCCCGCTGATACCCGCCAGCAACAGAACCCCGTACCCACTACGCGCCCCCCCCCCCATCAACGCAACACTTTGCGGGTTAATACCGCAGACCTCCACGACATACATGCCATCTCCCCACGCACAACCAACCCCTTCCAGCACTGCTCGCGccc                                                                          >  1:179912/1‑250 (MQ=255)                                                                                                                                                                                   gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTc                                                                      <  1:11546/250‑1 (MQ=255)                                                                                                                                                                                                                       ccTTGCTGGCCCCAGGGCTGGCCCC‑CCA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGGGCTGCTGGCCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGTTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAcc                                 <  1:2516/251‑1 (MQ=255)                                                                                                                                                                                                                              gcgcCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACCTCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCCGTTCATTGGGGCCATGTTGATGGGGCAAACCCACGACAACCCCCTTcc                          >  2:253627/1‑250 (MQ=255)                                                                                                                                                                                                                                   cAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTgg                     <  1:167869/251‑1 (MQ=255)                                                                                                                                                                                                                                         tGGCCCC‑ACA‑TCCCC‑C‑CCCCG‑TACC‑GGACTCGTGCTGCTGAACGCTGATATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACCTTGCGGTTTAATACCGCACGCATACACGCCATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGATCATTGGGGCCATGTTGATGAGTCAAAACAAAGAAAAACCCCTTCGCCAGGGGatt                >  1:109787/1‑250 (MQ=255)                                                                                                                                                                                                                                              cc‑ACA‑TCCCC‑C‑CCCCG‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATAAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtg          >  1:174333/1‑251 (MQ=255)                                                                                                                                                                                                                                                       ccc‑c‑ccc‑G‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTAc  <  1:25327/251‑1 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       GATATATAGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCC‑ACA‑TCCCC‑C‑CCC‑G‑TACG‑GTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTAC  >  NC_000913/3727868‑3728362

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