Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,291,288	G→A	G271G (GGG→GGA)	nrfD →	putative menaquinol‑cytochrome c reductase subunit NrfD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,291,288	0	G	A	100.0%	144.3 / NA	40	G271G (GGG→GGA)	nrfD	putative menaquinol‑cytochrome c reductase subunit NrfD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/19);  total (21/19)

TAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGT  >  NC_000913/4291210‑4291360                                                                               |                                                                         tAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCgggctgggact                                                                        <  1:859315/81‑1 (MQ=255)  aaaGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCgggctgggact                                                                        >  2:859391/1‑80 (MQ=255)    aGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGAt                                                                     <  2:859316/81‑1 (MQ=255)     gTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGAtt                                                                    <  1:859317/81‑1 (MQ=255)        cGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATtgtg                                                                 >  2:859371/1‑81 (MQ=255)         gTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGc                                                                >  2:859393/1‑81 (MQ=255)          tGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGcc                                                               <  1:859318/81‑1 (MQ=255)           gcgcTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGcc                                                               >  2:859384/1‑80 (MQ=255)             gcTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAAt                                                            >  1:859401/1‑81 (MQ=255)              cTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAAtg                                                           <  1:859319/81‑1 (MQ=255)               tggtggCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCAGGGCTGGGACTGATTGTGCCAAtgc                                                          >  2:859381/1‑81 (MQ=255)                 gtggCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAAtgct                                                         <  1:859321/80‑1 (MQ=255)                 gtggCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAAtgct                                                         <  2:859320/80‑1 (MQ=255)                   ggcggcATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATgctgct                                                      >  2:859397/1‑81 (MQ=255)                   ggcggcATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATgctgct                                                      >  1:859362/1‑81 (MQ=255)                     cggcATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCa                                                    <  1:859322/81‑1 (MQ=255)                       gcATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCaaa                                                  >  2:859402/1‑81 (MQ=255)                             ggTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTgg                                            >  1:859418/1‑81 (MQ=255)                              gTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTggg                                           >  2:859378/1‑81 (MQ=255)                              gTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTggg                                           <  2:859323/81‑1 (MQ=255)                               tGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGt                                          >  2:859388/1‑81 (MQ=255)                               tGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGt                                          <  2:859324/81‑1 (MQ=255)                                ggCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGt                                          >  2:859392/1‑80 (MQ=255)                                ggCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGt                                          >  2:859389/1‑80 (MQ=255)                                 gCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTc                                         <  2:859325/80‑1 (MQ=255)                                     tttCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCg                                    <  2:859326/81‑1 (MQ=255)                                        cTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCa                                  <  1:859327/80‑1 (MQ=255)                                          ggACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGtt                               >  1:859394/1‑81 (MQ=255)                                              gtggtgGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCgg                           >  2:859412/1‑81 (MQ=255)                                              gtggtgGTTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCgg                           <  2:859328/81‑1 (MQ=255)                                                 gtggtTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCAt                        >  2:859395/1‑81 (MQ=255)                                                  tggtTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCAtt                       <  1:859330/81‑1 (MQ=255)                                                  tggtTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCAtt                       <  2:859329/81‑1 (MQ=255)                                                   ggtTCTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTc                      >  2:859387/1‑81 (MQ=255)                                                       cTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTg                   <  2:859331/80‑1 (MQ=255)                                                       cTGGCTTGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTg                   <  1:859332/80‑1 (MQ=255)                                                             tGGTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCt            >  1:859423/1‑81 (MQ=255)                                                              ggTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCt            <  2:859333/80‑1 (MQ=255)                                                              ggTGTCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTg           >  2:859396/1‑81 (MQ=255)                                                               gtgtCGCCGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTgg          <  1:859334/81‑1 (MQ=255)                                                                      cGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCgtg   >  1:859399/1‑81 (MQ=255)                                                                      cGGGCTGGGACTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCgtg   >  1:859428/1‑81 (MQ=255)                                                                       gggctgggactgATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCgtgt  >  1:859409/1‑81 (MQ=255)                                                                               |                                                                         TAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGT  >  NC_000913/4291210‑4291360

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: ‑