Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,728,118	(C)8→9	coding (1254/1455 nt)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,110	1	.	C	81.5%	67.5 / 11.3	27	coding (1262/1455 nt)	xylB	xylulokinase
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (7/15);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑T‑C‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAC  >  NC_000913/3727978‑3728234                                                                                                                                       |                                                                                                                                tAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑t                                                                                                                                  >  1:435471/1‑133 (MQ=255)      ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑T‑c‑c‑ccc                                                                                                                          >  1:330451/1‑133 (MQ=255)       gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑T‑c‑c‑cccc                                                                                                                         >  1:45958/1‑133 (MQ=255)          tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑T‑C‑C‑CCCCC‑Cg                                                                                                                     <  1:17404/133‑1 (MQ=255)                        tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACCGTCATCTAg                                                                                                        >  2:499389/1‑133 (MQ=255)                          gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAgct                                                                                                      >  1:31178/1‑133 (MQ=255)                             gATGCGACCGTTCTAACGATAGCTGCGGCACCAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAgctgct                                                                                                   <  1:187177/133‑1 (MQ=255)                              aTGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGctgctg                                                                                                  <  2:306125/133‑1 (MQ=255)                               tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGa                                                                                                 <  2:8382/133‑1 (MQ=255)                               tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑C‑C‑C‑CCCCCACGTGCGATAATCGTGCTGCTGa                                                                                                 >  1:132581/1‑133 (MQ=255)                                 cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑T‑CTC‑CCCCC‑CGTACAGTCATCGAGGTGCTGAcc                                                                                               >  1:145394/1‑133 (MQ=255)                                   aCTGTTCTAATGGTAATTGCGGCATCAACTCACTGAGCGATATCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGc                                                                                             <  2:132581/133‑1 (MQ=255)                                    cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑T‑C‑CACCCGC‑CGTACGGTAATTGAGCAGCTGACCGCt                                                                                            >  2:187177/1‑133 (MQ=255)                                    cTGTTCTAACGGTAGATGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGCACCGTAATCGTGCTGCTGACCGCt                                                                                            >  2:514359/1‑133 (MQ=255)                                                 aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGc                                                                               >  2:109736/1‑133 (MQ=255)                                                 aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGc                                                                               >  1:26959/1‑133 (MQ=255)                                                     gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                           <  2:404322/133‑1 (MQ=255)                                                         caacaaTTCACTGAGCGATTTCTCCGGAGTCGACGCGATCTCCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                       <  2:145394/133‑1 (MQ=255)                                                         caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGGCAGCATCTGACg                                                                       <  2:105602/133‑1 (MQ=255)                                                         caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                       <  1:386348/133‑1 (MQ=255)                                                             aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAg                                                                   <  2:435471/133‑1 (MQ=255)                                                                 cAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGGTGATATCGGCCAGCGTCTGACTCCAGTACt                                                               >  1:135175/1‑133 (MQ=255)                                                                  aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTc                                                              <  2:483589/133‑1 (MQ=255)                                                                        cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATT‑C‑C‑CCCCC‑CGTACGGTAATCGACCCTCTGTCCGCTGATTTTCGCCAGCATCTCACGCCAGTAATCACAAcg                                                        >  2:454350/1‑133 (MQ=255)                                                                                       cgccgcGATCTGCGCCAGCCTTGCTGCGCTCAGTTCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTCCGGTAATCGAGCTGCTGTCCGCTGATATCCGCCAGCATCTGACGCCAGTGCTCCCTACGCGCCCCGCCCCCAAt                                         <  2:415095/133‑1 (MQ=255)                                                                                       cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                         <  1:514359/133‑1 (MQ=255)                                                                                       cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                         <  1:315168/133‑1 (MQ=255)                                                                                                      cAGCCTTGCTGCGCCCAGTGCTGGCCCGACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACGCTTTg                          <  1:454350/133‑1 (MQ=255)                                                                                                                     cAGTGCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa           >  2:330989/1‑133 (MQ=255)                                                                                                                         gCTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa       <  2:330451/133‑1 (MQ=255)                                                                                                                         gCTGGCCCCACA‑T‑C‑CACCCCC‑CGTCCGGCAATCGAGCTGCTTACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa       <  2:135175/133‑1 (MQ=255)                                                                                                                          cTGGCCCCACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt      <  2:26959/133‑1 (MQ=255)                                                                                                                              ccccACA‑TCC‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc  <  1:330989/133‑1 (MQ=255)                                                                                                                                       |                                                                                                                                TAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACA‑T‑C‑C‑CCCCC‑CGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAC  >  NC_000913/3727978‑3728234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑