Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,053,455	A→G	intergenic (+37/+192)	hemN → / ← yshB	coproporphyrinogen III dehydrogenase/uncharacterized protein YshB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,053,455	0	A	G	100.0%	90.1 / NA	26	intergenic (+37/+192)	hemN/yshB	coproporphyrinogen III dehydrogenase/uncharacterized protein YshB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/21);  total (5/21)

CCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTT  >  NC_000913/4053379‑4053529                                                                             |                                                                           ccAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAgc                                                                        <  2:798671/81‑1 (MQ=255)    gAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGcc                                                                     <  1:798672/81‑1 (MQ=255)       aGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATc                                                                  <  1:798673/81‑1 (MQ=255)        gcgcGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATc                                                                  >  2:798712/1‑80 (MQ=255)        gcgcGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCa                                                                 <  2:798674/81‑1 (MQ=255)           cgGATGCAGGAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAAtg                                                              <  2:798675/81‑1 (MQ=255)             gATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAAtgtg                                                            <  2:798676/81‑1 (MQ=255)                  agcagTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTgag                                                       <  1:798677/81‑1 (MQ=255)                       ttCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACaa                                                  <  2:798678/81‑1 (MQ=255)                              gggTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGAtt                                              <  2:851478/78‑1 (MQ=255)                              gggTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTg                                             <  1:798679/79‑1 (MQ=255)                                      tAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATgcgc                                   >  2:798716/1‑81 (MQ=255)                                        aaTAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCt                                  <  2:798680/80‑1 (MQ=255)                                          tAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTggg                               >  2:798719/1‑81 (MQ=255)                                          tAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTggg                               <  1:798681/81‑1 (MQ=255)                                           aGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTggg                               >  2:798710/1‑80 (MQ=255)                                            gTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGc                              <  2:798682/80‑1 (MQ=255)                                            gTGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGc                              <  2:798683/80‑1 (MQ=255)                                             tGAATGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGcaa                            <  2:798684/81‑1 (MQ=255)                                                 tGGCGCTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACaaa                        <  1:798685/81‑1 (MQ=255)                                                      cTTCGTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTg                   <  2:798686/81‑1 (MQ=255)                                                          gTTTACAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTg                <  2:798687/80‑1 (MQ=255)                                                               cAAAGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAatt           <  1:798688/80‑1 (MQ=255)                                                                aaaGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAattat         <  1:798689/81‑1 (MQ=255)                                                                 aaGTAGAAAACGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAattatt        >  1:798727/1‑81 (MQ=255)                                                                       aaaaCGGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCtttt  <  1:798690/81‑1 (MQ=255)                                                                             |                                                                           CCAGAAAGCGCGGATGCAGCAGTTCTCTCGGGTGATTTAAATAGTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTT  >  NC_000913/4053379‑4053529

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