Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	949,985	G→A	W106* (TGG→TGA)	ycaK →	putative oxidoreductase YcaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	949,985	0	G	A	100.0%	60.8 / NA	18	W106* (TGG→TGA)	ycaK	putative oxidoreductase YcaK
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (3/15);  total (3/15)

CAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCA  >  NC_000913/949749‑950215                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                       cAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGggctggc                                                                                                                                                                                                                          <  1:157225/251‑1 (MQ=255)      ggCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATg                                                                                                                                                                                                                     <  2:157468/251‑1 (MQ=255)                                                 cGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTgg                                                                                                                                                                           <  1:22092/250‑1 (MQ=255)                                                            cAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTggag                                                                                                                                                               <  2:68950/251‑1 (MQ=255)                                                               cccAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGgaggag                                                                                                                                                            <  2:213678/251‑1 (MQ=255)                                                                   gTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCGGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCGATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAAGTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACaa                                                                                                                                                        <  1:387358/251‑1 (MQ=255)                                                                                           gACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATTGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTg                                                                                                                                <  1:32522/251‑1 (MQ=255)                                                                                                         ggATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGTGCTACTTGAACATGACACGTTAGTGGTGGTTTTTCATCTCTGGTGGTACAGCTTCCCCGCACTGCTAACAGGATATCTTGACAGAGTTTGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAATGTTCGATGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTACAGATGGGCTGGGAACAAAAtata                                                                                                                   >  1:49348/1‑250 (MQ=255)                                                                                                                cGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAg                                                                                                                                                      <  1:67189/208‑1 (MQ=255)                                                                                                                cGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAg                                                                                                                                                      >  2:67189/1‑208 (MQ=255)                                                                                                                 gTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGattat                                                                                                          <  1:77358/251‑1 (MQ=255)                                                                                                                             agGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTTTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTTGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATAtg                                                                                              <  2:53841/251‑1 (MQ=255)                                                                                                                                          ttagtCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGTTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCt                                                                                   <  2:259229/245‑1 (MQ=255)                                                                                                                                                                gACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGAGTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCa                                                            <  2:388967/250‑1 (MQ=255)                                                                                                                                                                                tttttCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATAc                                            >  1:242640/1‑250 (MQ=255)                                                                                                                                                                                        ctctGGTGGTACATCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGATATGGGCACAAATGACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTa                                     <  1:402467/249‑1 (MQ=255)                                                                                                                                                                                                    aGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGATATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGTGCTGGGAAAAAAATATAAGCGATGATTTAAAAAATATGTGCAGTTATCTTGTTATTGAAGATGCCGATGTCACTTTCTTGTGTAAGACAGTTGTAGTCGAGGCGgaaga                       <  2:265719/251‑1 (MQ=255)                                                                                                                                                                                                                         ggATATATTGACAGAGTATGAAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCa  <  2:156323/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                       CAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCA  >  NC_000913/949749‑950215

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