Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,411,028	C→T	T159I (ACC→ATC)	yjfJ →	PspA family protein YjfJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,411,028	0	C	T	100.0%	56.4 / NA	20	T159I (ACC→ATC)	yjfJ	PspA family protein YjfJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/12);  total (8/12)

AAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACTGCAACGCCAACAGGGCGAGTAAT  >  NC_000913/4410788‑4411252                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 aaGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGACGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAACTGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGgcgc                                                                                                                                                                                                                         >  1:179172/1‑250 (MQ=255)                  gTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGAc                                                                                                                                                                                                       >  2:70346/1‑251 (MQ=255)                                   tAACGAAGTTGCTGAATAAATCGCGCGCCTTGAGAATCTCATTACCGGTGAATAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTa                                                                                                                                                                                       <  2:11764/250‑1 (MQ=255)                                             gCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCa                                                                                                                                                                            <  2:3576/251‑1 (MQ=255)                                                  agaaATCGCGCGCCTTGCGAATCTCATAACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCATGATGGCGTGGAAAACGCAGTTAAAGCGACAGCGCAGCGTATCGCTCAGTTTGCGCAGCAACTGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTACCACCCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGACTCCTTAACACGCCTGCACAcgc                                                                                                                                                                       >  1:52933/1‑251 (MQ=255)                                                                         tcATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTGACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGt                                                                                                                                                         <  2:228762/242‑1 (MQ=255)                                                                         tcATTACCGCTGAAGAGCAACTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGt                                                                                                                                                         >  1:228762/1‑242 (MQ=255)                                                                            ttACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATAGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATTTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTg                                                                                                                                             >  1:124088/1‑251 (MQ=255)                                                                                          cAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCTCTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATTTTTCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCATACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGcaca                                                                                                                               <  2:124088/251‑1 (MQ=255)                                                                                               gttgtCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTgg                                                                                                                          <  2:135379/251‑1 (MQ=255)                                                                                                                            gCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGAc                                                                                             <  1:4453/251‑1 (MQ=255)                                                                                                                             cGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACg                                                                                            <  1:30623/251‑1 (MQ=255)                                                                                                                                      aGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAGATGAAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAAAAACCTCTATCGTTGGCGCATCTCCCAGCGTTTCGACAGCGGCAGACTCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACCGTTGGAGCGAGTCGCAGACGGTCACGCCCTTGaacaatagctgg                                                                                   >  2:3445/1‑240 (MQ=255)                                                                                                                                                    cagcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCTCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAAtt                                                                     <  1:184628/251‑1 (MQ=255)                                                                                                                                                     agcgcagcgTATCGCTCATTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCGTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTg                                                                    <  2:207690/251‑1 (MQ=255)                                                                                                                                                                  gCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCCCTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGCATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCGGGAGCAATa                                                        >  1:162001/1‑250 (MQ=255)                                                                                                                                                                                            cGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCa                              <  2:182201/250‑1 (MQ=255)                                                                                                                                                                                                 aaGCCACTGAAGCCATGCAGCGTGCACAACCGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACTg                        >  1:84167/1‑251 (MQ=255)                                                                                                                                                                                                               aTGCAGCGTGCACAACAGGCTGTAACAACCTCTATCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACTGCAACGCCAACAggg          <  2:98076/251‑1 (MQ=255)                                                                                                                                                                                                                        gCACAACTGGCTGTATCAACCTCTATCGTTTGCGCATCTTCCAGCGTTTCGACAGCGGCATAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCAATACTGCAACGCCAACAGGGCGAGTAAt  <  2:3505/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 AAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGGCTCGCCTGGATGCTGCCGCACAGTTGGAGAAAGTCGCAGACGGTCGCGACCTTGACGAAAAGCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACTGCAACGCCAACAGGGCGAGTAAT  >  NC_000913/4410788‑4411252

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