Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	17,385	C→A	intergenic (+375/‑104)	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,385	0	C	A	100.0%	42.1 / NA	14	intergenic (+375/‑104)	sokC/nhaA	small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/4);  total (10/4)

TCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTA  >  NC_000913/17169‑17601                                                                                                                                                                                                                           |                                                                                                                                                                                                                          tctcTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGATATTCCTACACTATAATc                                                                                                                                                                                          >  1:193165/1‑251 (MQ=255)            cAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACg                                                                                                                                                                               <  1:182297/251‑1 (MQ=255)                      cgacAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTcc                                                                                                                                                                                                      >  1:15264/1‑218 (MQ=255)                              tGGATTATTTTTGGAATATTGTCCTAACAAATAACGACGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCACAAAA‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTACTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTa                                                                                                                                                              >  2:96204/1‑250 (MQ=255)                                                     ctaacaaataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCAc                                                                                                                                     >  2:226260/12‑251 (MQ=255)                                                                  ccGACTGACAAGACATTAATTATTAATTTTACTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATATGGACTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGG‑TACAAATAGTACAAACGATCTATTCACCTGAAAGCGAAAt                                                                                                                        >  1:190029/1‑251 (MQ=255)                                                                                                          cccTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGc                                                                                >  1:9626/1‑251 (MQ=255)                                                                                                          cccTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGACGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGc                                                                                >  2:113017/1‑251 (MQ=255)                                                                                                              cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTATGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCACCGATTCTTTACCAGTGATGCCTc                                                                             <  2:24472/246‑1 (MQ=255)                                                                                                               agaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCgg                                                                           >  2:186550/4‑251 (MQ=255)                                                                                                                                 cTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc                                                          >  2:74961/1‑250 (MQ=255)                                                                                                                                         aTCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCt                                                 >  2:75436/1‑251 (MQ=255)                                                                                                                                             gggCTTAGGGGCACTCGTTAGCGCTTACAGCCGTCA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTA‑AAATAGTAAAAACGTTCTATTCACCTGAAAGAGAAATAAGAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGGCTCGGGAGGCATTATTCTTATCATTGCCGCt                                                 <  1:20650/246‑1 (MQ=255)                                                                                                                                                                                         cATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGAGTTGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTa  <  1:151168/251‑1 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                          TCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCA‑AAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTA‑AAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTA  >  NC_000913/17169‑17601

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