Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	17,039	C→T	intergenic (+29/‑450)	sokC → / → nhaA	small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	17,039	0	C	T	100.0%	23.0 / NA	11	intergenic (+29/‑450)	sokC/nhaA	small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/1);  total (10/1)

ACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGA‑TGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTT  >  NC_000913/16833‑17262                                                                                                                                                                                                                |                                                                                                                                                                                                                                aCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTAAGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGGTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAACACCCGCCGCTATTATATCGCTCTCTTTAAcccct                                                                                                                                                                                     >  1:6746/1‑249 (MQ=255)               cgcgggggtaaatcggATGACCATCAGGGCTACATTCACCTCCTTCTGTTGCTTCATTCCCCTCTTCCCCTTCATCTCCCGGTCGGAAAGGGCCCCTCTACATGTTTTCAGTAGATGGGGGGCCTCGG‑TTGAATAGAAAAATATCCCCCGGGTTTTTTTTCTATCTCCCGTTCAGATAATCCCTAAGACATATAGCCAAAACCTCCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGAt                                                                                                                                                                       <  1:522452/236‑1 (MQ=255)                  ggcggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCt                                                                                                                                                                    >  1:204489/1‑251 (MQ=255)                                              aCAATCATCGCCTTATGCTGCTTCATTGATATCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGACAACGCCTCGCAtttttttg                                                                                                                                        >  2:395651/1‑248 (MQ=255)                                                                     cATTGCTCTCTTCTACGTGAACTTACGGTCAGTAAGAGGCACTCTACCTGTGTTCAGCATATAGAGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTACTCCTGAAGACGCATCGCATATCTGAGGCGTAATTTTTTGATGATTTAat                                                                                                                 >  2:470459/1‑251 (MQ=255)                                                                            ctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTa                                                                                                           >  1:300693/1‑250 (MQ=255)                                                                                    cTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGctctct                                                                                                                                                                                              >  1:62723/1‑159 (MQ=255)                                                                                                        gagGCACTCTTCATGTGGTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCATGAGCCAGATAGCCTCAAGCACCCGCCGCTATTATCTCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGCCGCCGCATCGAATTCTTGTGGCGAATTTTTTTACTGCTTTAATTATTTAACTTTATTTTATATCtttatagtaatttt                                                                              >  2:150418/1‑239 (MQ=255)                                                                                                                                                ttGA‑TAGTAAAAAATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTGTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCATAGCATTTTTGTGGCGTACTTTTTTCATGGTTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGaaaacttgtcttaac                                       >  2:249995/1‑238 (MQ=255)                                                                                                                                                                                    cTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGAACCCGCCGCTATTATATCGCTCTCTTTAAACCATTCTGTTTTATCGATTCTATTCCTGAAGACGCATCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACACCAAGCTGGCTTATTTTTGAAATATTGTCCTAACAAATACCGACGACTGACAAGAAATTAATTATTACtt   >  2:8731/1‑251 (MQ=255)                                                                                                                                                                                     tGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACttt  >  1:161050/1‑251 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                                ACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGA‑TGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTT  >  NC_000913/16833‑17262

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