Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,605,791	A→G	intergenic (‑128/‑13)	fhuF ← / → yjjZ	ferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,605,791	0	A	G	100.0%	80.4 / NA	25	intergenic (‑128/‑13)	fhuF/yjjZ	ferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/10);  total (15/10)

CGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGG  >  NC_000913/4605554‑4606035                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     cgctGTTTTGCGATCGTCGCACACCCGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGCGACGCGTTCGCCCGATGACCTCTTCATAGAGCGGTGCGGAACGATAGGCCCTAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATc                                                                                                                                                                                                                                          >  2:127972/4‑250 (MQ=255)         tGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCa                                                                                                                                                                                                                                  >  2:202733/1‑250 (MQ=255)                       gCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGc                                                                                                                                                                                                                     >  1:37626/1‑249 (MQ=255)                           ggCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGAGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGc                                                                                                                                                                                                               >  1:124559/1‑251 (MQ=255)                                  gcAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAACTCATGTTGCAACGTACGCTGGGCAGTGGCTgggg                                                                                                                                                                                                          >  1:57200/1‑249 (MQ=255)                                                    ggctggcaataagtttgCCAGATGACATCTTCATAGAGCGGTGCGGACCGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTACTCCTAGCGCTAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCACATGATACAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGc                                                                                                                                                                                      >  2:243126/17‑251 (MQ=255)                                                                                               cGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTg                                                                                                                                           <  1:148515/251‑1 (MQ=255)                                                                                                 gAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGTTCAGTGGc                                                                                                                                         >  1:127094/1‑251 (MQ=255)                                                                                                         ggCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTAt                                                                                                                                                                  >  1:253948/1‑218 (MQ=255)                                                                                                         ggCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAAGTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTAt                                                                                                                                                                  <  2:253948/218‑1 (MQ=255)                                                                                                                                    gataatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTAATGACGcc                                                                                                       <  2:194828/250‑1 (MQ=255)                                                                                                                                           atTGCTAATCCTAGCGATAGGTTTACCCGAGAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTAATTCTGCTGGGATTAGTAATGACGCCGATGATGc                                                                                               >  1:258969/1‑251 (MQ=255)                                                                                                                                                                     ccGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGTGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTAGCGGCATTAc                                                                     <  1:110932/251‑1 (MQ=255)                                                                                                                                                                                  aTTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGCGTGTTGCTGCCGGGATTGCTGCTTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACCGTTGCGGCATTACATTTTGCTACCac                                                        >  2:89121/1‑250 (MQ=255)                                                                                                                                                                                          ggCTTGCAAATGATAAAAATTATCATATGCTATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTTACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTCACATCGTGCCTg                                                >  2:112855/1‑251 (MQ=255)                                                                                                                                                                                                cAAATGATAAAACTTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGAAGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACtt                                          >  2:189165/1‑251 (MQ=255)                                                                                                                                                                                                     gATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGc                                     <  2:57200/251‑1 (MQ=255)                                                                                                                                                                                                           aaTTATCATATGATATTGATTATCATTATCACTGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGAAGAGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCGTTCTGATGGGTTTAGTATTGAAGCCGATGGTGCTGTATCACAAACCATTGCGTCATTTCCTTTTGCTGCCATCCTGACTGGCACTTATTGCTGGcat                               >  2:27251/1‑251 (MQ=255)                                                                                                                                                                                                                 cATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTTGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCAGTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGc                          <  1:108565/250‑1 (MQ=255)                                                                                                                                                                                                                      gATATTGGTTATCATTATCAATGAGAGAGATGAACTCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTCACATCGTGCCTGGCACTTATTGCTGGCCTCATGCtgatg                     >  2:33728/1‑250 (MQ=255)                                                                                                                                                                                                                                  cattatcaATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGTATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTTCATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCAGCATGCTGATGATAATGAATTTg         <  2:108234/250‑1 (MQ=255)                                                                                                                                                                                                                                  cattatcaATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCTGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTg         <  1:243126/250‑1 (MQ=255)                                                                                                                                                                                                                                  cattatcaATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTg         <  1:89965/250‑1 (MQ=255)                                                                                                                                                                                                                                   agtagcaATGAGAGAGATGAAATCAGGGTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATGGCTGGCATCATGCTGATGATAATGAATTTGaa       <  1:127972/246‑1 (MQ=255)                                                                                                                                                                                                                                        caATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGACTCAgg  >  2:62291/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     CGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGG  >  NC_000913/4605554‑4606035

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