Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,729,955	2 bp→AA	coding (810‑811/1323 nt)	xylA ←	xylose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,729,955	0	G	A	100.0%	75.7 / NA	26	H271Y (CAC→TAC)	xylA	xylose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/16);  total (10/16)
*	NC_000913	3,729,956	0	G	A	100.0%	68.7 / NA	24	N270N (AAC→AAT)	xylA	xylose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/16);  total (8/16)

TTATCAGTACTTTGACGACGTACTTTGGCATCGAAGTTCAGACCACCGGTGGTGAAACCGCCTGCTTTGAGAATTTCATACATCACCAGCGCATTCTCTTCCACACTGTTCGGGAACTGGTCGGTGTCCCAGCCCAGTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTGGTTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCCTGGAAACCGATTTTATGTTTATGCTCAACCACCATCTGCATAAAGCGGCCCAGTTGTTCACGCTCCTGACGCAAGTCGGTATTTAACAGCGTTTCGTAACCTTCAC  >  NC_000913/3729720‑3730194                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                               ttATCAGTACTTTGACGACGTACTTTGGCATCGAAGTTCAGACCACCGGTGGTGAAACCGCCTGCTTTGAGAATTTCATACATCGCCAGCGCATTCTCTTCCACACTGTTCGGGAACTGGTCGGTGTCCCAGCCCAGTTGCGCATCGCCACGGGTGGCGTCGACATAACCGAACAGGCCAAGCGCAATGGCGTGGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGt                                                                                                                                                                                                                                   <  1:203482/250‑1 (MQ=255)                   cgtgcTTTTGCATCGAAGTTCAGACCACCGGGGGTGAAACCGCCTGCTTTGAGAATTTCATACATCACCAGCGCATTCTCTTCCACACTGTTCGGGAACTGGTCGGTGGCCCAGCCTAGTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGTCCAAGCGCAATGGCGGTGGCTATTTCAGGATGTAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTc                                                                                                                                                                                                                <  1:314810/247‑1 (MQ=255)                                      tCAGACCTCCGGTGGTGAAACCGCCTGCTTTGAGAATTTCATACATCACCAGCGCATTCTCTTCCACACTGTTCGGGAACTGGTCGGTGTCCCAGCCCAGTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCAGGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAg                                                                                                                                                                                              <  1:65172/250‑1 (MQ=255)                                                   ggtgAAACCGCCTGCTTTGAGAATTTCATACATCACCAGCGCATTCTCTTCCACACTGTTCGGGAACTGGTCGGTGTCCCAGCCCAGTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTGTTTCCAGACCAAACTGTTTCAGGAATCCATAGACNg                                                                                                                                                                                <  1:348185/251‑1 (MQ=255)                                                                                                    tCCACACTGTTCGGGAACTGGTCGGTGTCCCAGCCCAGTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTc                                                                                                                                 <  1:231709/249‑1 (MQ=255)                                                                                                       acacAGTTCGGGAACTGGTCGGTGTCCCAGCCCAGTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTc                                                                                                                            <  1:382410/251‑1 (MQ=255)                                                                                                               cGGGAACTGGTCGGTGTCCCAGCCCAGTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCATACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGc                                                                                                                     <  1:215140/250‑1 (MQ=255)                                                                                                                             tggcccagcccagTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTGCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACGGTGTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTc                                                                                                                                 <  1:406001/221‑1 (MQ=255)                                                                                                                                cccagcccagTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGAGGGAAAGAGGGACCTGCCAGCGTCGCGTAATTAGCTTCAATTTTCAGTTTAATCTCTTGTTCCATACCAAACTGTTTCAGGAAGCCATAGACCGTCTCGGCATCTTAAGCATATTGAGGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCCTGGAAAcc                                                                                                    <  1:288518/250‑1 (MQ=255)                                                                                                                                  cagcccagTGGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGTGCAATGGCGGTGGCGATTTCATGATGGAAAGAGGGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAGCTGTTTCGGGACGCCATATACCGGCGCGGCATCGTAATCAGATGGAGGTTTGTTCGGTTATTGCGGTGTCGGTTCGATAAGCAACGTGCCCTGGAAACCGAt                                                                                                 <  1:381157/251‑1 (MQ=255)                                                                                                                                            gcgcATCGCCACGGTTGGCGTCGACAGAACCCAACAGGCCAAGCGCAATGGCGGTGCCTCTTTCATGATGGCAAGAGTGACATGCCAGCGTCGCGAAATTAGCTTCAATGTTCAGTTTAATCTCttt                                                                                                                                                                                                                   >  1:135401/1‑127 (MQ=255)                                                                                                                                               cATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCACCGTCGCGCAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAg                                                                                                                      >  1:421208/1‑217 (MQ=255)                                                                                                                                                                  gACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTGCAGTTTAATCTCTTTGTCCAGACCAAGCTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTgg                                                                                                                                                  <  1:220008/170‑1 (MQ=255)                                                                                                                                                                        aCCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGAGAAGCAACGTGCCCTGGAAACCGATTTTATGTTTATGCTCAACCACCATCTGCATAAAGCGGc                                                           <  1:399109/251‑1 (MQ=255)                                                                                                                                                                           naaCAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGAGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTAAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCCTCGTAATCAAATTGATGTTTGTTCGGGTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCCTGGAAACCGATTTTATTTTTATGCGCAACCACCATCTGCATAAAGCGGccc                                                         >  1:281384/2‑250 (MQ=255)                                                                                                                                                                                                 gTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCCTGGAAAcc                                                                                                    >  1:177541/1‑185 (MQ=255)                                                                                                                                                                                                       aTTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCACTGTTCAGTTTAATCTGTTTTTCCAGACAAAACTGTTTCAGGAAGACATAGAACGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTAGGTTCGCTAAGCAACGTGCCCTGGAACCCGATTTTATGTTTATGCTCaa                                                                                >  1:231656/1‑199 (MQ=255)                                                                                                                                                                                                           cATGANGGAACNANNCACCTGCCAGCGTCGCGTANTTAGCTTCAATGTTCNGTTTNATCTCntantn                                                                                                                                                                                                                >  2:170941/1‑63 (MQ=255)                                                                                                                                                                                                                 ggAAAGAGTGACCTGCCAGCGTCGCGTAATTAGCTTCAAGGTGCAGTTTAATCTCTTTTTCCAGACCAAACTGGTTCAGGAAGCCATAGCCCGGCGCGGGATCGTAATCATATTGATGTGTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCATGGAAAc                                                                                                     <  1:209719/168‑1 (MQ=255)                                                                                                                                                                                                                     agagTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCAGAGAcc                                                                                                                                                                                 <  1:342962/88‑1 (MQ=255)                                                                                                                                                                                                                        gTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCCGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGCAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCCTGGAAACCGATTTTATGTTTATGCTCAACCACCATCTGCATAAAGCGGCCCAGTTGTTCACGCTCCTggc                                      >  1:288886/1‑222 (MQ=255)                                                                                                                                                                                                                        gTGACCTGCCAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCtttttn                                                                                                                                                                                                                <  1:128745/54‑2 (MQ=255)                                                                                                                                                                                                                        gTGACCTGCCAGCATCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCCTGGAAACCGATTTTATGTTTATGCTCAACCACCATCTGCATAAAGCGGCCCAGTTGCTCACGCTCCTGACGCAAGTCGGTATTTAACAGCGTTTCGt           >  1:361233/1‑251 (MQ=255)                                                                                                                                                                                                                           aCCTGCCAGCGTCGCGTAATTAGCTTNAATGTTNNGTNTAATCTCTTTTTCCAGACCNAANNNTNNNANGAAGCNATNNACCGTCGCGNCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCNACGTGcc                                                                                                             >  2:94073/1‑150 (MQ=255)                                                                                                                                                                                                                                 cAGCGTCGCGTAATTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCCTGGAAACCGATTTTATGTTTATGCTCAACCACCATCTGCATAAAGCGGCCCAGTTGTTCACGCTCCTGACGCAAGTCGGTATTTAACAGCGTTTCGTAACCTTCAc  <  1:425146/251‑1 (MQ=255)                                                                                                                                                                                                                                        gcgTANTTAGCNTNNATGTTCAGTTTAATCTCTTNTTCCAGACCAAACTGNTTCANGAAGCNaa                                                                                                                                                                                      >  2:168116/1‑63 (MQ=255)                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                               TTATCAGTACTTTGACGACGTACTTTGGCATCGAAGTTCAGACCACCGGTGGTGAAACCGCCTGCTTTGAGAATTTCATACATCACCAGCGCATTCTCTTCCACACTGTTCGGGAACTGGTCGGTGTCCCAGCCCAGTTGCGCATCGCCACGGTTGGCGTCGACAGAACCGAACAGGCCAAGCGCAATGGCGGTGGCTATTTCATGATGGAAAGAGTGACCTGCCAGCGTCGCGTGGTTAGCTTCAATGTTCAGTTTAATCTCTTTTTCCAGACCAAACTGTTTCAGGAAGCCATAGACCGTCGCGGCATCGTAATCATATTGATGTTTGGTCGGTTCTTGCGGTTTCGGTTCGATAAGCAACGTGCCCTGGAAACCGATTTTATGTTTATGCTCAACCACCATCTGCATAAAGCGGCCCAGTTGTTCACGCTCCTGACGCAAGTCGGTATTTAACAGCGTTTCGTAACCTTCAC  >  NC_000913/3729720‑3730194

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