Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,310,770	A→G	H32H (CAT→CAC)	yciI ←	protein YciI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,310,770	0	A	G	100.0%	66.8 / NA	20	H32H (CAT→CAC)	yciI	protein YciI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/9);  total (11/9)

GTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATC  >  NC_000913/1310544‑1311009                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                gTGCCCTTAAATATTGAATCTCTATTCAGAACTCTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTaaa                                                                                                                                                                                                                         <  2:323998/251‑1 (MQ=255)            tatTGAAACTCTCTTCAGAACACTTTCTTAAATGGCTTCACTGAAACGTGTTCATCGACTCCTGCCGCTACGTACGGGTCCGCATCGGCCCCGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGATAGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTCAATGTGCCGGACg                                                                                                                                                                                                              >  2:223220/1‑251 (MQ=255)                gAATCGCTATTCAGAACACTTTCTTAAATGGTTCCACTGAAACGTGGTCAGAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAAc                                                                                                                                                                                                           <  1:361501/250‑1 (MQ=255)                            aGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCAAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGtttt                                                                                                                                                                                              >  1:66364/1‑251 (MQ=255)                             gAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGttttt                                                                                                                                                                                             >  1:16170/1‑251 (MQ=255)                                             ggTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCtt                                                                                                                                                                              <  1:52677/250‑1 (MQ=255)                                                         aCGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGt                                                                                                                                                                  >  1:53232/1‑250 (MQ=255)                                                                                 acgtacgGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGGTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGTCACaa                                                                                                                                          <  2:19670/250‑1 (MQ=255)                                                                                         gTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGtt                                                                                                                                 <  1:268101/251‑1 (MQ=255)                                                                                                                   ctTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTAttt                                                                                                       >  1:23127/1‑251 (MQ=255)                                                                                                                       cAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCgg                                                                                                                              >  1:23494/1‑224 (MQ=255)                                                                                                                       cAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCgg                                                                                                                              <  2:23494/224‑1 (MQ=255)                                                                                                                             ttcaactacAGCAATAAGGGTTGCGCCAGGCAATCCAGCAGCCCCTGGATCGTTACTGTCTACAGCTGGAATTGGTCCAGCTGTCAACAAACGACCTTCATCGCGCAGTAACTGTAAACGTGCTAAATGTGCCGGGCGAACGGAAAGGCGTATTTCGAGGCTATCAACTTTGTTTTGAGCGAAAATAACTTCCACCACGGGCACAACTCCTTGTTCGGGGTAGTTGTAAATTACGTTATGTAAAAGGGCg                                                                                              >  2:152814/9‑250 (MQ=255)                                                                                                                                        cAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGt                                                                                                                        >  1:28912/1‑213 (MQ=255)                                                                                                                                        cAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAAc                                                                                  >  1:213091/1‑251 (MQ=255)                                                                                                                                        cAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTGTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGt                                                                                                                        <  2:28912/213‑1 (MQ=255)                                                                                                                                               ggTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCATCATCGTGCAGTAACTGTAAACGAGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGaa                                                                              >  1:164162/1‑248 (MQ=255)                                                                                                                                                                           ggATCGTTACTGTCTACCGCTGGCCTTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAACTGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCCACGGAAAGATGATGTCTTTGTTAAGGCCATGCATaa                                                >  2:364832/1‑250 (MQ=255)                                                                                                                                                                                                ggCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCaaa                           <  2:66364/250‑1 (MQ=255)                                                                                                                                                                                                                        cGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCActgatcctgatc  <  1:224014/251‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                GTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATC  >  NC_000913/1310544‑1311009

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