Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,083,432	Δ1 bp	coding (57/1272 nt)	efeB →	heme‑containing peroxidase/deferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,083,430	0	A	.	100.0%	102.4 / NA	25	coding (55/1272 nt)	efeB	heme‑containing peroxidase/deferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (16/9);  total (16/9)

GTGCTGGGACTGGATTAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTGAAAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGCACAAAAAACGCAAAGTGCGCC  >  NC_000913/1083353‑1083500                                                                              |                                                                       gtgCTGGGACTGGATTAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTGaagt                                                                     >  2:219647/1‑79 (MQ=255)      gggACTGGATTAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAg                                                               >  1:219638/1‑81 (MQ=255)       ggACTGGATTAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAgg                                                              <  1:219610/81‑1 (MQ=255)           aggATTAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCa                                                          >  2:219636/2‑81 (MQ=255)            ggATTAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCAc                                                         >  2:219635/1‑81 (MQ=255)                tAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCActggc                                                     >  2:219660/1‑81 (MQ=255)                tAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCActggc                                                     >  2:219651/1‑81 (MQ=255)                   gCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCActggcgct                                                  >  2:219641/1‑81 (MQ=255)                        aTGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCgg                                             >  2:219645/1‑81 (MQ=255)                        aTGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCgg                                             >  2:219640/1‑81 (MQ=255)                                  aaGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTcc                                   <  2:219613/81‑1 (MQ=255)                                  aaGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTcc                                   <  1:219611/81‑1 (MQ=255)                                  aaGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTcc                                   <  1:219612/81‑1 (MQ=255)                                  aaGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTc                                    <  1:219614/80‑1 (MQ=255)                                     aTGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGt                                <  2:219615/81‑1 (MQ=255)                                       gAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCg                              >  2:219644/1‑81 (MQ=255)                                          aaCGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCt                            >  1:219646/1‑80 (MQ=255)                                          aaCGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTc                           <  1:219616/81‑1 (MQ=255)                                           aCGGCGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCa                          >  1:219652/1‑81 (MQ=255)                                               cGTGAATGAACCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGca                      >  2:219659/1‑81 (MQ=255)                                                    atgaaCCGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGCACaaa                  >  2:219654/1‑80 (MQ=255)                                                         ccGTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGCACAAAAAACGc            <  2:219617/81‑1 (MQ=255)                                                           gTCACGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGCACAAAAAACGCaa          <  2:219618/81‑1 (MQ=255)                                                               cGCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGCACAAAAAACGCAAAGTg      >  2:219655/1‑81 (MQ=255)                                                                gCCGACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGCACAAAAAACGCAAAGTgc     >  2:219657/1‑81 (MQ=255)                                                                   gACGTTTACTG‑AAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGCACAAAAAACGCAAAGTGCGcc  >  2:219648/1‑81 (MQ=255)                                                                              |                                                                       GTGCTGGGACTGGATTAAGCGTTATGCAGTATAAAGATGAAAACGGCGTGAATGAACCGTCACGCCGACGTTTACTGAAAGTGATAGGTGCACTGGCGCTGGCGGGAAGTTGTCCGGTCGCTCATGCACAAAAAACGCAAAGTGCGCC  >  NC_000913/1083353‑1083500

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