Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	612,670	A→G	intergenic (‑176/‑67)	fepA ← / → fes	ferric enterobactin outer membrane transporter/enterochelin esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	612,670	0	A	G	100.0%	97.4 / NA	27	intergenic (‑176/‑67)	fepA/fes	ferric enterobactin outer membrane transporter/enterochelin esterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/8);  total (19/8)

AGACACGCAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAACATGGGGTTATGGTGTGACGGCGT  >  NC_000913/612592‑612746                                                                               |                                                                             aGACACGCAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTc                                                                            <  2:122346/81‑1 (MQ=255)  gACACGCAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCa                                                                           <  2:122347/81‑1 (MQ=255)   acacGCAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCaa                                                                          <  1:122348/81‑1 (MQ=255)   acacGCAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCaa                                                                          >  2:122357/1‑81 (MQ=255)    cacGCAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCaat                                                                         <  1:122349/81‑1 (MQ=255)            gcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAAtattat                                                                 >  2:122360/1‑81 (MQ=255)             cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAAtattat                                                                 >  1:122371/1‑80 (MQ=255)                ttGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCaa                                                              >  1:122365/1‑80 (MQ=255)                 tGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAata                                                            >  1:122373/1‑81 (MQ=255)                               ccGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGcaatc                                               >  2:122368/1‑80 (MQ=255)                               ccGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGcaatc                                               >  1:122378/1‑80 (MQ=255)                                     ttCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGaa                                         >  2:122363/1‑80 (MQ=255)                                      tCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGaaa                                        >  2:122364/1‑80 (MQ=255)                                        ccACGGCGCGCCATTACACTATTGCAAATGCAAATAGTTGTCAATAATATTATCAAGATATTTCTGCAATCAATGAAAAAt                                     >  2:122385/1‑81 (MQ=255)                                          aCGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTg                                   >  1:122376/1‑81 (MQ=255)                                          aCGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTg                                   >  1:122359/1‑81 (MQ=255)                                           cGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGc                                  <  1:122350/81‑1 (MQ=255)                                            ggCGCGCCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGca                                 >  1:122374/1‑81 (MQ=255)                                                 gcCATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTa                            <  1:122351/81‑1 (MQ=255)                                                   cATTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTaaa                          >  1:122370/1‑81 (MQ=255)                                                    aTTACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAAc                         <  2:122352/81‑1 (MQ=255)                                                     ttACGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAACa                        <  1:122353/81‑1 (MQ=255)                                                        cGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAACATgg                     >  1:122377/1‑81 (MQ=255)                                                        cGCTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAACATgg                     >  2:122372/1‑81 (MQ=255)                                                          cTATTGCAAATGCAAATAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAACATgggg                   >  1:122379/1‑81 (MQ=255)                                                                caaatgcaaatAGTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAACATGGGGTTATgg             >  1:122369/1‑81 (MQ=255)                                                                            gTTGTCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAACATGGGGTTATGGTGTGACGGCGt  >  2:122381/1‑80 (MQ=255)                                                                               |                                                                             AGACACGCAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCTGCAATCAATGAAAAATTGCACAGTAAACATGGGGTTATGGTGTGACGGCGT  >  NC_000913/612592‑612746

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