Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	658,165	C→T	intergenic (‑227/‑782)	crcB ← / → tatE	F(‑) channel protein/twin arginine protein translocation system ‑ TatE protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	658,165	0	C	T	100.0%	82.9 / NA	25	intergenic (‑227/‑782)	crcB/tatE	F(‑) channel protein/twin arginine protein translocation system ‑ TatE protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/12);  total (13/12)

TGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCGCCCGTTAT  >  NC_000913/657939‑658353                                                                                                                                                                                                                                   |                                                                                                                                                                                             tGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCg                                                                                                                                                                      >  1:25038/1‑251 (MQ=255)     aaaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGTCTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTTCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTa                                                                                                                                                                  <  1:47796/251‑1 (MQ=255)                          aTGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAg                                                                                                                                             <  1:175103/251‑1 (MQ=255)                                    aGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCt                                                                                                                                   <  1:99106/251‑1 (MQ=255)                                    aGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCt                                                                                                                                   <  1:172030/251‑1 (MQ=255)                                    aGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCGGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGTTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTccg                                                                                                                                   <  1:172031/251‑2 (MQ=255)                                             ccGAATAATAGTCGTTGGCGACTTTTACGACTATGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGTTCTATCGGTTAAATCAGCACAGCTGCTGAAAGGCGAATTCCTCGGCCttt                                                                                                                           >  2:7686/1‑250 (MQ=255)                                                    atCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCgag                                                                                                                   <  1:142569/251‑1 (MQ=255)                                                           ttGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTa                                                                                                             >  1:2045/1‑250 (MQ=255)                                                                 aaTTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTa                                                                                                       >  1:124716/1‑250 (MQ=255)                                                                     ttACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACatg                                                                                                  >  2:87408/1‑251 (MQ=255)                                                                                    gaggTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGGTACATCTGTGTGTGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGTACTGTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTg                                                                                   <  2:124716/251‑1 (MQ=255)                                                                                             aTGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCAt                                                                          <  1:248486/251‑1 (MQ=255)                                                                                                     tGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATGCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCtt                                                                   <  2:65689/250‑1 (MQ=255)                                                                                                      gCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTAAGAATTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACCACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCtta                                                                  >  1:91905/1‑249 (MQ=255)                                                                                                             gACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCgg                                                          >  1:257821/1‑251 (MQ=255)                                                                                                                aGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAAAGCCggggg                                                       >  2:39417/1‑250 (MQ=255)                                                                                                                                tgtgtgGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGACGCAGGAAGCGGGAAACGACACATCGCTGTTTGCCCTGCCGAAAGCAGTGCTGCGGCGTGATGCACATGATGCAGCCCTATCGGTTAAATCCGCACAGCTGCTGTAAGGCGAATTCCTCGGAATTTACGGCGAGACAGAAACCGTAACATGATGTCTACAATTTAGACGATTCATGTTCCTTCAACGCCGTGGCGCGCATGGAATATGCt                                       >  1:14647/1‑251 (MQ=255)                                                                                                                                                  tGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTcaggcagg                     >  1:82856/1‑251 (MQ=255)                                                                                                                                                         tgtgCCTCGTTGCTGGCGCAGGCGGCGGAACACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAACTCAGCACAGCTGCTGGACGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTACCATGATGCCATCAATTCTGACGCTTCATGTTCCTTCaa                                                                >  1:88769/1‑201 (MQ=255)                                                                                                                                                                cGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCGcc        <  2:25038/250‑1 (MQ=255)                                                                                                                                                                    gATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCGCCCGtt    <  1:156774/250‑1 (MQ=255)                                                                                                                                                                     aTGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCCCACCTGCTGGAAGGCGAATTCCTCGGACTTTCCGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGCCGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGCATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCGCCCGTTAt  >  1:178960/1‑251 (MQ=255)                                                                                                                                                                                               aTCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACAtcgtcg          >  1:160109/1‑217 (MQ=255)                                                                                                                                                                                               aTCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGTGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACAtcgtcg          <  2:160109/217‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                             TGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCGCCCGTTAT  >  NC_000913/657939‑658353

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