Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,310,841	C→T	D9N (GAC→AAC)	yciI ←	protein YciI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,310,841	0	C	T	100.0%	68.8 / NA	20	D9N (GAC→AAC)	yciI	protein YciI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/10);  total (10/10)

ACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTGCCTTATTGAATATGATTGCTATTTGCATTTAAAATCGAGACCTGGTTTTTCTACTGAAATGATTATG  >  NC_000913/1310600‑1311082                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  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                                                                                                                                                                                                       <  2:23494/224‑1 (MQ=255)                                                               cAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCgg                                                                                                                                                                                                       >  1:23494/1‑224 (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)                                                                                                                                                                                  aaCTGTAAACGTGCTAAATGGGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGGTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGATCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGGCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTGCCTTATTGAat                                                          <  1:286266/250‑1 (MQ=255)                                                                                                                                                                                      gTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAAGGACATTTTCACTGATCCTGATCGTCGTGCCTTATTGAATATGa                                                      <  2:53232/250‑1 (MQ=255)                                                                                                                                                                                           cGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTGCCTTATTGAATATGATTGCt                                                 >  1:257324/1‑250 (MQ=255)                                                                                                                                                                                                                  ggAAAGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTGCCTTATTGAATATGATTGCTATTTGCATTTAAAATCGAGACCTg                         <  1:319054/251‑1 (MQ=255)                                                                                                                                                                                                                      aGGCGTTTTTCGAGGCTATCAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTGCCTTATTGAATATGATTGCTATTTGCATTTAAAATCGAGACCTGGttt                     >  2:190348/1‑251 (MQ=255)                                                                                                                                                                                                                                         cAGCTTTGTTTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTGCCTTATTGAATATGATTGCTATTTGCATTTAAAATCGAGACCTGGTTTTTCTACTGAAATGATTATg  <  2:213091/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  ACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTGCCTTATTGAATATGATTGCTATTTGCATTTAAAATCGAGACCTGGTTTTTCTACTGAAATGATTATG  >  NC_000913/1310600‑1311082

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