Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,292,161	+G	coding (650/1659 nt)	nrfE →	putative cytochrome c‑type biogenesis protein NrfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,292,159	1	.	G	100.0%	81.9 / NA	24	L216L (CTC→CTG)	nrfE	putative cytochrome c‑type biogenesis protein NrfE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (4/20);  total (4/20)

GAGTCTACTGCGCGGCGAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTC‑GGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTC  >  NC_000913/4292036‑4292284                                                                                                                             |                                                                                                                              gAGTCTACTGCGCGGCGAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCtgg                                                                                                                <  1:1034456/140‑1 (MQ=255) gAGTCTACTGCGCGGCGAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCtgg                                                                                                                <  1:677796/140‑1 (MQ=255) gAGTCTACTGCGCGGCGAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCtgg                                                                                                                <  2:362688/140‑1 (MQ=255)    tCTACTGCGCGGCGAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCtggtgg                                                                                                             <  1:435684/140‑1 (MQ=255)      tACTGCGCGGCGAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGtt                                                                                                           <  1:931463/140‑1 (MQ=255)         tGCGCGGCGAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTg                                                                                                        <  1:406388/140‑1 (MQ=255)             cggcgAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGCGCATCCTCCGTCCCTGGGGGGGGCACTGCGCACACGGGGTGGGGGGGGTGGGGGGTGTggggg                                                                                                    >  2:313869/1‑138 (MQ=255)                        aTGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGaaaa                                                                                         <  1:433913/140‑1 (MQ=255)                           gTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGc                                                                                      <  2:404721/140‑1 (MQ=255)                           gTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGc                                                                                      <  2:395536/140‑1 (MQ=255)                           gTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGc                                                                                      <  1:623864/140‑1 (MQ=255)                           gTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGc                                                                                      <  2:799231/140‑1 (MQ=255)                             gCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCct                                                                                    <  1:591787/140‑1 (MQ=255)                             gCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCct                                                                                    <  2:13408/140‑1 (MQ=255)                               cTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCctct                                                                                  <  1:196266/140‑1 (MQ=255)                                                        acagtACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTg                                                         <  1:549320/136‑1 (MQ=255)                                                                  gctggAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGca                                               <  1:2624/140‑1 (MQ=255)                                                                  gctggAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGca                                               <  1:590467/140‑1 (MQ=255)                                                                      gAGTGCATTAACGGCGGGGGGCATCCTCCGTTCCTGGGGGGGCTACCGCGAACACGGGCTGGGGGGGGGGGGGGGTTTGGGGTCCGGGGGGGAACGCCCCTTTTTTTTCCTCGCCTTCTTCCCCCCCGCCGCTGcccacg                                           >  2:871776/1‑138 (MQ=255)                                                                                             tcCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCgggg                    <  1:56617/140‑1 (MQ=255)                                                                                                tCGGTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGAt                 >  1:705281/1‑140 (MQ=255)                                                                                                  ggTTCCTGGTGGGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGAttt               >  1:214795/1‑140 (MQ=255)                                                                                                        tggtggGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGcc         >  1:360067/1‑140 (MQ=255)                                                                                                        tggtggGCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGcc         >  1:1071658/1‑140 (MQ=255)                                                                                                            gggCCTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTg     <  1:468134/140‑1 (MQ=255)                                                                                                               ccTACTGCGAACTCGGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAGAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTc  <  1:150764/140‑1 (MQ=255)                                                                                                                             |                                                                                                                              GAGTCTACTGCGCGGCGAGTTTGATGGTGCCTGCGCCCGAATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTC‑GGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTC  >  NC_000913/4292036‑4292284

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