Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,733,008	G→A	V270I (GTT→ATT)	xylG →	xylose ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,733,008	0	G	A	85.7%	21.7 / ‑0.3	14	V270I (GTT→ATT)	xylG	xylose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (4/8);  total (4/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

GATTTGCGTTATTCGCGACGGACAGCACATTGGTACGCGTGATGCTGCCGGAATGAGTGAAGACGATATTATCACCATGATGGTCGGGCGAGAGTTAACCGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGGTTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGTTGGTGCCGGACGTACCGAGACCATTCAGTGCCTGTTTGGTGTGTGGCCCGGACAATGGGAAGGAAAAATTTATATTGATGGCAAACAGGTAGATA  >  NC_000913/3732833‑3733182                                                                                                                                                                                |                                                                                                                                                                               gATTTGCGTTATTCGCGACGGACAGCACATTGGTACGCGTGATGCTGCCGGAATGAGTGAAGACGATATTATCACCATGATGGTCGGGCGAGAGTTAACCGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGATTAATCGTCATATTAAACGAGTTAatgat                                                                                                                                                   <  1:190209/205‑1 (MQ=255)                                               gcCGGAATGAGTGAAGACGATATTATCACCATGATGGTCGGGCGAGAGTTAACCGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGATTAATCGTCATATTAAACTAGTTAATGATGTCTCTTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGTTGGTGc                                                                                           <  1:293163/215‑1 (MQ=255)                                                  ggAATGAGTGAAGACGATATTATCACCATGATGGTCGGGCGAGAGTTAACCGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGGTTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGTTGGTGCCTGACGTACCGAGACCATTc                                                                       <  1:91268/232‑1 (MQ=255)                                                               aCGATATTATCACCAGGATGGTCGGGCGAGAGTTAACCGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGATTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGTTGTTGCCGGACGTACCGAGACCATTCAGTGCCTGTTTGGTGTGTGGCCCGGACAATgg                                       <  1:265585/251‑1 (MQ=255)                                                                                       ggCGAGAGTTAACCGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCAGGGCATCCGATTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGan                                                                                                     <  1:37294/165‑2 (MQ=255)                                                                                         cGAGAGTTAACCGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGATTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGTTGGTGCCGGACGTACCGAGAc                                                                            >  1:289210/1‑188 (MQ=255)                                                                                                   ccGCGCTTTANCCTAATGAACCACATACCACCGGAGATGAAANANTACGTATTGAACATCTGACGGCATGGCATCCGATTAATCGTCATATTANACGAGTTAATNATGNNTCGTNTNCNCTGAAACGTGGCGAAATATTGGGTATTGCNGGACTCGTTGGTGCCGGACGTACNNAGNNNATNNAGTGc                                                                  <  2:243090/188‑1 (MQ=255)                                                                                                   ccGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGATTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGTTGGTGCCGGACGTACCGAGACCATTCAGTGc                                                                  >  1:243090/1‑188 (MQ=255)                                                                                                        cTTTACCCTAATGAACCACATACCACNGGAGATNNACNATTACGTATTGAACATCTGNCGNNNTNNNANACGATNAANNGTCATATTANACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGc                                                                                                                           >  2:74287/1‑126 (MQ=255)                                                                                                        cTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGATTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGc                                                                                                                           <  1:74287/126‑1 (MQ=255)                                                                                                             cccTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGATTAAGCGTCATATTAAACGATTTAATGATCTCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGtt                                                                                                <  1:292577/148‑1 (MQ=255)                                                                                                              ccTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGGTTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTCGCGAAATATTGGGTATTGCCGGACTCGTTGGTGcc                                                                                          <  1:141178/153‑1 (MQ=255)                                                                                                                                                       tGAACNTCTGGNGNNATGGCATCCGATTAATCGTNATATTAAACGAGTTANTGATNGCTCGNTTNc                                                                                                                                        >  2:135190/1‑66 (MQ=255)                                                                                                                                                       tGAACATCTGCCGGCATGGCATCCGATTAATCGTCATATTAAACGATTTAATTATGGCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTAGTGCCGGACTCGGTTGTGTCTGACGTACCGAGACTATGTAGGGCCTGTTTAGTGTGTGGCCCGGACAACGGGAAGGAAAAACTTATATTGATGGCAAACAGGTAGata  <  1:135190/200‑1 (MQ=255)                                                                                                                                                                               ggTTAATCGTCATATTAAACGAGTCAATGATGTCTCGTTTGCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGTGGGTGCCGGACGTACCGAGACCATTCAGTGCCTGTTTGGTGTGTGGCCCGGACAATGGGAAGGAAAAATTTATATTGATGGc               <  1:164365/163‑1 (MQ=255)                                                                                                                                                                                |                                                                                                                                                                               GATTTGCGTTATTCGCGACGGACAGCACATTGGTACGCGTGATGCTGCCGGAATGAGTGAAGACGATATTATCACCATGATGGTCGGGCGAGAGTTAACCGCGCTTTACCCTAATGAACCACATACCACCGGAGATGAAATATTACGTATTGAACATCTGACGGCATGGCATCCGGTTAATCGTCATATTAAACGAGTTAATGATGTCTCGTTTTCCCTGAAACGTGGCGAAATATTGGGTATTGCCGGACTCGTTGGTGCCGGACGTACCGAGACCATTCAGTGCCTGTTTGGTGTGTGGCCCGGACAATGGGAAGGAAAAATTTATATTGATGGCAAACAGGTAGATA  >  NC_000913/3732833‑3733182

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