Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,084,827	G→A	intergenic (+180/‑165)	efeB → / → phoH	deferrrochelatase, periplasmic/ATP‑binding protein, putative PhoH family P‑loop ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,084,827	0	G	A	100.0%	149.3 / NA	48	intergenic (+180/‑165)	efeB/phoH	deferrrochelatase, periplasmic/ATP‑binding protein, putative PhoH family P‑loop ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/26);  total (22/26)

ACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTGTCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATGCGAGATCTGTTGTACTTATTAAGTAGCAGCGGAAGTTCCCGG  >  NC_000913/1084692‑1084960                                                                                                                                        |                                                                                                                                      aCGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTAtca                                                                                                                                    >  2:327449/1‑139 (MQ=255)   gCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTatcatc                                                                                                                                  <  2:671471/139‑1 (MQ=255)    cTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTAtcatca                                                                                                                                 <  1:737067/139‑1 (MQ=255)    cTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTAtcatca                                                                                                                                 <  1:36620/139‑1 (MQ=255)        cgcgTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCActct                                                                                                                             <  2:105099/139‑1 (MQ=255)             cTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCAt                                                                                                                        <  1:530051/139‑1 (MQ=255)                  cATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTc                                                                                                                   <  1:327449/139‑1 (MQ=255)                          cAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGaa                                                                                                           >  2:421824/1‑139 (MQ=255)                               aGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTaa                                                                                                      >  1:107273/1‑139 (MQ=255)                               aGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTaa                                                                                                      <  2:393417/139‑1 (MQ=255)                                 tGCGGGTTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATg                                                                                                    >  1:229600/1‑139 (MQ=255)                                 tGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATg                                                                                                    >  1:766413/1‑139 (MQ=255)                                  gCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGt                                                                                                   >  1:730900/1‑139 (MQ=255)                                          ggTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGaaa                                                                                           <  1:192675/139‑1 (MQ=255)                                           gTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAAt                                                                                          >  1:733289/1‑139 (MQ=255)                                                cTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGgtgt                                                                                     >  1:362060/1‑139 (MQ=255)                                                cTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGgtgt                                                                                     >  2:430149/1‑139 (MQ=255)                                                   gATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGtttt                                                                                  >  1:81670/1‑139 (MQ=255)                                                    aTAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTta                                                                                 <  1:430149/139‑1 (MQ=255)                                                    aTAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTta                                                                                 <  1:598898/139‑1 (MQ=255)                                                    aTAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTta                                                                                 <  2:630439/139‑1 (MQ=255)                                                     tAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTta                                                                                 >  2:556720/1‑138 (MQ=255)                                                     tAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTta                                                                                 <  1:556720/138‑1 (MQ=255)                                                      aaGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTtata                                                                               >  1:522077/1‑139 (MQ=255)                                                        ggCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGt                                                                             <  2:733289/139‑1 (MQ=255)                                                          cGTTCACGTCGCATCCGGCGATCGTGCACTAATGCCTGATGCAAATCGTACTGAAAGGACATAGCTTTGTTCATCGCTAGCATTACTTTTTCTTCTTTCCATTAGAAACTAATGTAACGTAAACGGTGTTTTATCGTTa                                                                           >  2:542314/1‑139 (MQ=255)                                                              cACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAata                                                                       <  1:421824/139‑1 (MQ=255)                                                                                     aCTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATa                                                <  2:766413/139‑1 (MQ=255)                                                                                         atgcctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATg                                                                                         <  2:578758/94‑1 (MQ=255)                                                                                         atgcctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATg                                                                                         >  1:578758/1‑94 (MQ=255)                                                                                         atgcctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAatat                                                                      <  2:112663/113‑1 (MQ=255)                                                                                         atgcctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAatat                                                                      >  1:112663/1‑113 (MQ=255)                                                                                         atgcctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATg                                            >  1:107517/1‑139 (MQ=255)                                                                                         atgcctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATg                                            >  2:403687/1‑139 (MQ=255)                                                                                         atgcctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAA‑‑ATATTGTTGCAATAAATGCg                                          >  1:695544/1‑139 (MQ=255)                                                                                         atccctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATg                                            >  1:358787/4‑139 (MQ=255)                                                                                           gcctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAA‑‑ATATTGTTGCAATAAATGCgag                                        <  2:695544/139‑1 (MQ=255)                                                                                             ctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATGCgag                                        <  2:362060/139‑1 (MQ=255)                                                                                             ctgatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATGCgag                                        <  1:403687/139‑1 (MQ=255)                                                                                               gatgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATGCGAGAt                                      <  2:358787/139‑1 (MQ=255)                                                                                                 tgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCa                                                   >  2:394312/1‑124 (MQ=255)                                                                                                 tgcAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCa                                                   <  1:394312/124‑1 (MQ=255)                                                                                                     aaTCCTGCTGAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATGCGAGATCtgttg                                <  2:644091/139‑1 (MQ=255)                                                                                                           gctgAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAg                                                                              <  1:684104/87‑1 (MQ=255)                                                                                                           gctgAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAg                                                                              >  2:684104/1‑87 (MQ=255)                                                                                                           gctgAAAGCACACAGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATGCGAGATCTGTTGTACtta                          <  2:107273/139‑1 (MQ=255)                                                                                                                      acaGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATg                                                                                         >  1:19129/1‑65 (MQ=255)                                                                                                                      acaGCTTTTTTCATCACTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATg                                                                                         <  2:19129/65‑1 (MQ=255)                                                                                                                                   tcaCTATCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATGCGAGATCTGTTGTACTTATTAAGTAGCAGCGGAAGTTCCCgg  <  2:229600/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      ACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACACAGCTTTTTTCATCACTGTCATCACTCTGTCATCTTTCCAGTAGAAACTAATGTCACTGAAATGGTGTTTTATAGTTAAATATAAGTAAATATATTGTTGCAATAAATGCGAGATCTGTTGTACTTATTAAGTAGCAGCGGAAGTTCCCGG  >  NC_000913/1084692‑1084960

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