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%	122.8 / NA	36	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 (15/21);  total (15/21)

TGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTG‑CCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGA  >  NC_000913/2449114‑2449266                                                                               |                                                                            tgtTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGc                                                                           <  2:435775/81‑1 (MQ=255) tgtTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGc                                                                           <  2:435774/81‑1 (MQ=255)      attgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAgggg                                                                      >  2:435811/1‑81 (MQ=255)      attgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAggg                                                                       <  2:435776/80‑1 (MQ=255)       ttgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGc                                                                     <  2:435777/81‑1 (MQ=255)           aGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAAt                                                                 <  2:435778/81‑1 (MQ=255)            gTCGAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATg                                                                >  2:435810/1‑81 (MQ=255)              cGAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGcc                                                              >  2:435813/1‑81 (MQ=255)               gAAGAGTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGccc                                                             >  1:435814/1‑81 (MQ=255)                 agagTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCcct                                                           >  1:435817/1‑81 (MQ=255)                 agagTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCcct                                                           <  1:435779/81‑1 (MQ=255)                   agTGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCcctgc                                                         >  2:435822/1‑81 (MQ=255)                     tGGCTGCCGCCGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCcctgcct                                                       <  1:435780/81‑1 (MQ=255)                          gccgccGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGgtg                                                   >  2:435823/1‑80 (MQ=255)                          gccgccGGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGgtg                                                   >  2:435829/1‑80 (MQ=255)                            cgccgGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTAt                                                <  1:435782/81‑1 (MQ=255)                            cgccgGAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTAt                                                <  2:435781/81‑1 (MQ=255)                              ccgGAGTTGCCCCTGAGGAGTTGAGCGCGTTACGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCt                                              <  2:435783/81‑1 (MQ=255)                                ggAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCt                                              >  1:435832/1‑79 (MQ=255)                                 gAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAg                                           <  2:435784/81‑1 (MQ=255)                                 gAGTTGCCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAg                                           <  2:435785/81‑1 (MQ=255)                                          cTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAAc                                  >  2:435812/1‑81 (MQ=255)                                            gaggagTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACt                                 >  2:435825/1‑80 (MQ=255)                                             aggagTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACt                                 <  2:435786/79‑1 (MQ=255)                                             aggagTTGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCa                               >  1:435855/1‑81 (MQ=255)                                                  ttGAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAg                          <  1:435787/81‑1 (MQ=255)                                                    gAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGt                         <  2:435788/80‑1 (MQ=255)                                                    gAGCGCGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTa                        >  1:435821/1‑81 (MQ=255)                                                      gcgcGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCTCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATg                      >  2:435820/1‑81 (MQ=255)                                                      gcgcGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATg                      <  2:435789/81‑1 (MQ=255)                                                      gcgcGTTGCGCCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATg                      >  2:435831/1‑81 (MQ=255)                                                              cgcCAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGaa              <  1:435790/81‑1 (MQ=255)                                                                ccAGGTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAAcc            <  1:435791/81‑1 (MQ=255)                                                                   ggTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGc         <  2:435792/81‑1 (MQ=255)                                                                    gTACGGGTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCa        <  2:435793/81‑1 (MQ=255)                                                                          gTAATAGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGa  <  1:435794/81‑1 (MQ=255)                                                                               |                                                                            TGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTG‑CCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGA  >  NC_000913/2449114‑2449266

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