Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	16,966	A→G	noncoding (15/59 nt)	sokC →	small regulatory RNA antitoxin SokC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	16,966	0	A	G	100.0%	68.1 / NA	25	noncoding (15/59 nt)	sokC	small regulatory RNA antitoxin SokC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/7);  total (18/7)

ATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACA  >  NC_000913/16726‑17194                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                       agggaacccccccccGCCGTTGCCCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGt                                                                                                                                                                                                                              <  1:171849/240‑1 (MQ=255)  gggtatccccccccGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGtt                                                                                                                                                                                                                             <  1:456250/241‑1 (MQ=255)         aCTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTggtagt                                                                                                                                                                                                                       >  1:15036/1‑245 (MQ=255)                           tACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTTCGGTCAGTAAGAGGCACTCTACCTGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCtt                                                                                                                                                                                                    >  1:484162/1‑251 (MQ=255)                                 gATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTTACCTTACGGTCCGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTc‑tct                                                                                                                                                                                             >  2:343961/1‑251 (MQ=255)                                                        aaCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGCTCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTg                                                                                                                                                                      >  2:237777/1‑251 (MQ=255)                                                                tCTGGCCAGTTCGGATGTGAACCTCACAGAGGACTTTTCTCGTTACCAGCGCTGCCACTACGGCGGTGATAAAGATGACGATCAGGGCGACAATCCTCGCCTTATGCTGCTTCCATGC‑TCTCTTCCCCTTGACCTTACGGTCAGGAAGAGGCACTCTCCATGGGTTCCGCATATAGGGGGCCTCGGGGTGATAGGAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGagac                                                                                                                                                                  >  1:554746/1‑247 (MQ=255)                                                                                      cTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGCTACAGCTGACGATAAGGGCGCCAATCATCGCCTTATGCTGATTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCCGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCAATAGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCtct                                                                                                                                        >  1:388717/1‑249 (MQ=255)                                                                                                     tctcGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGGAAATCATCGCCTTAAGCTGCTTCAATGC‑TCTCTTCTCCTTGACCTTACGGACAGCAAGAGGCACCCTACATGTGCTGAGCATATAGGGGGACTCGGGATGATAGGAACATATCACTCTGGGACTTTC‑TCTATCTGCCGTTCCGGTAATGCCTGAGACATATAGCCTCCAGCACCCGCCGCTATTATCTCgtttttttt                                                                                                                         >  1:608440/1‑243 (MQ=255)                                                                                                           tACCAGCGCCGCCACTACGGCGGTGATAAAGATGAAGAACAGGGAGACAAACAACGCCTTATTCTGCTTCAATGA‑TCTCTTCTCCTTTACCCTACGGTCAGTAAAAGGCACTCTACATGTGGTCCGCCTAAATGGGGCCTCGGGTTGGTTGTAAAATATCACTCAGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGAAAGACTCAATCACCCGCCGCTATTATATCTCTCCCCTTTAAccc                                                                                                                   >  2:445256/1‑251 (MQ=255)                                                                                                               ggctccgccACTACGGCGGTGATACAGATGACGATCGGGGCGACAATCATCGCCTTATGCTGCTGCATTGC‑TCTCTGCTCCGTGACCTTACTGTCAGTAATAGTCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCa                                                                                                                                                                                                              <  2:306943/153‑1 (MQ=255)                                                                                                               aGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCCCCTTATGCTGCTTCCTTGC‑TCTCTTCTCCTTGAACTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCa                                                                                                                                                                                                              >  1:306943/1‑157 (MQ=255)                                                                                                                                         gATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTCTCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGAc                                                                                    >  1:254263/1‑251 (MQ=255)                                                                                                                                          aTGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGA‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGATGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTTAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCCTTCTGTTTTATCGATTCTAATCCTgtcgacg                                                                                   >  2:476368/1‑245 (MQ=255)                                                                                                                                                 tCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGc                                                                             >  1:216223/1‑250 (MQ=255)                                                                                                                                                       cGACAAACACCGCTTTATGCTGCTTC‑TTGCCTCTTTTCCCCTTGCCCTTACGGTAAGAAAGAGGCACTCTACATGTGTTCAGCATAAGGGGGGCCGGGGTTTGAAAATAAAATATCCCCGGGGGCTTTTC‑TTTTTCTGCCGTTCAGCTAAGGCTGGAGACATATAGCCTCAACCCCCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTtg                                                                      <  1:263560/251‑1 (MQ=255)                                                                                                                                                        gACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTtgt                                                                     >  1:591731/1‑251 (MQ=255)                                                                                                                                                                       atcctgctTCTTTTC‑TCCCTTCCCCTTGCCCTTGGGGTAGGTAAGAGCCACTCTACTTGTTTTCAGCATAGAGGGGGCCTCGGGTGAATAAAAAAATATCAC‑CGGGGGTTTTCCTTTCCCTCCCGTTCATTTAATGCCTGAAACAGACAGCCTAATGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTaa                                                      <  2:431476/248‑1 (MQ=255)                                                                                                                                                                         gctgctTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATg                                                    >  1:205211/1‑251 (MQ=255)                                                                                                                                                                           gctttcttttc‑tctctTCTCCTTGTCTTTACGGTAATTAAGAGCCCCTCTACATTTGTCCAGCATATGGGGGCCCTGGGGTTAATAGTAAAATATCCCTCGGGTTTTTTC‑TCTATCTGCGGTTCACCTAATCCCGGAGACAGATACCCACAAGCATCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt                                                  <  1:13536/241‑1 (MQ=255)                                                                                                                                                                                                  tGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAAttt                             >  1:209526/1‑250 (MQ=255)                                                                                                                                                                                                    aCCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATCTAtc                          >  1:5222/1‑251 (MQ=255)                                                                                                                                                                                                           ggTCAGTAAGAGGCACTCTACTTGTGTTCGGAATAAAGGGGGCCGCGGTTTGAAAGTAAAATTTCCCGCGGGGTTTTTT‑TCTATCTGCCGTGCAGCTAAGGCCGGAAACAGAAACCCACAATCGTCCGCCGCTATTATATCGCTCTC‑TTTAACCCAGTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTAAATTATTTAACTTTAATTTATCTCtttat                   <  1:269246/251‑4 (MQ=255)                                                                                                                                                                                                                  aaGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAatt            >  1:291636/1‑251 (MQ=255)                                                                                                                                                                                                                             tACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACa  >  1:51532/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                       ATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGC‑TCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTC‑TTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACA  >  NC_000913/16726‑17194

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