Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,449,191	G→A	intergenic (+34/+37)	smrB → / ← yfcO	putative endonuclease SmrB/DUF2544 domain‑containing protein YfcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,449,191	0	G	A	100.0%	81.9 / NA	24	intergenic (+34/+37)	smrB/yfcO	putative endonuclease SmrB/DUF2544 domain‑containing protein YfcO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/15);  total (9/15)

ATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAAT  >  NC_000913/2449119‑2449268                                                                         |                                                                              attgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAgggg                                                                        <  1:387571/80‑1 (MQ=255) attgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGc                                                                       <  1:387572/81‑1 (MQ=255)       gTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGc                                                                 <  2:387573/81‑1 (MQ=255)           aaGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGGGCGCGTTGCGCCAGGTACGGGTAATAGCAGGAGCGAATGCCcct                                                             <  2:387574/81‑1 (MQ=255)            agagTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCcctg                                                            <  1:387575/81‑1 (MQ=255)             gagTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCcctgc                                                           >  2:387612/1‑81 (MQ=255)               gTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCcctgcct                                                         <  1:387576/81‑1 (MQ=255)                tGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCcctgcctg                                                        <  2:387577/81‑1 (MQ=255)                 ggCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTgg                                                       >  1:387609/1‑81 (MQ=255)                      ccgccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTAt                                                  <  2:387578/81‑1 (MQ=255)                         ccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTa                                               <  1:387579/81‑1 (MQ=255)                           ggAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCTGGGGCGAATGCCCCTGCCTGGTGTATCTAAg                                             >  2:387611/1‑81 (MQ=255)                                      gaggagTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACt                                   >  2:387620/1‑80 (MQ=255)                                         gagTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCAtt                               >  1:387616/1‑81 (MQ=255)                                          agTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTc                              <  2:387580/81‑1 (MQ=255)                                             tGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAg                            >  1:387624/1‑80 (MQ=255)                                              gAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGt                           <  2:387581/80‑1 (MQ=255)                                                 cgcgTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGt                       <  1:387582/81‑1 (MQ=255)                                                     ttGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCTCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTGTGTACTg                   >  1:387614/1‑81 (MQ=255)                                                     ttGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTg                   >  1:387613/1‑81 (MQ=255)                                                        cgcCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGa                 <  1:387583/80‑1 (MQ=255)                                                                  gggTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTa      <  1:387584/81‑1 (MQ=255)                                                                     tAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGaa   <  2:387585/81‑1 (MQ=255)                                                                      aaTAGCAGTGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTATAAt  >  2:387617/1‑81 (MQ=255)                                                                         |                                                                              ATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAAT  >  NC_000913/2449119‑2449268

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