Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,510,238	T→G	intergenic (+445/+452)	insO → / ← fecE	IS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,510,238	0	T	G	100.0%	120.5 / NA	35	intergenic (+445/+452)	insO/fecE	IS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/20);  total (15/20)

CCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGT  >  NC_000913/4510005‑4510483                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      ccgctGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCAGTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGGAGCAAACTCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAGTCCGTTCCAGGGTTGCTTTTGATTTAc                                                                                                                                                                                                                                      <  2:31146/248‑1 (MQ=255)      ggctaaaaaGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCGTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATACCAATCGATCCATCAACGGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACTCACACTGTGTAGAAAAAGCGATGGTCCGAAACTGGATGAGCAATTCCGTTCCAGGGCTGCTTTTGTTTTACGtt                                                                                                                                                                                                                                   <  2:385385/245‑1 (MQ=255)        cGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCg                                                                                                                                                                                                                               <  1:37478/251‑1 (MQ=255)            aaaGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTg                                                                                                                                                                                                                           <  2:118138/251‑1 (MQ=255)             aaGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACTCCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGa                                                                                                                                                                                                                          <  1:113416/251‑1 (MQ=255)               gggAGAAGGACGGGCCCGTCTGCCTTCACGTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGATGCTTTTGATTTACGTTGGCGTCTGATc                                                                                                                                                                                                                        <  2:255409/251‑1 (MQ=255)                    aaGGGCGGGCCTCTCTTCCTTCACTTTATGCCAGGTACGGTAAGCTGGGTGAGACCCCGTCTCATGAAGATATCAGTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGTCTGCCGGTGATTAAGAAAGCGGCACCTACCTAAGTAGGAGCAAACGCGCACTTTTTAGAAAAATCTATGGTCAGAAACTGGATTAGCAATTCCGTGCCAGGGTTGCGTTTGAGTTACGTTGGCGTCTGATCATTGa                                                                                                                                                                                                                   <  2:280357/251‑1 (MQ=255)                           ggCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGCTTTATcc                                                                                                                                                                                                            >  1:280570/1‑251 (MQ=255)                                         cACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAAc                                                                                                                                                                                               <  1:270226/250‑1 (MQ=255)                                                    aGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGGATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAAt                                                                                                                                                                                   <  2:344758/251‑1 (MQ=255)                                                                   ggTGAGACACCGACTCATGAAGAGATCATTGATAACAATCGATCCATCAACTGGCCGGTATAACGTGGCTGCCGGTTATTAAGAAAGCTGCACCTACCGAAGGAGTAGCAAACGCGCACGTTTTAGAAAAATCGATGCTCAGAACCGTGATTAGCAATTCCGTGCCAGGGTGGCTGTTGATGTTCGTTGGCGGCTGATCATTGATTTAACGCGAAGAGTTCACGCTCTTTGGTAATGAACCAGCTCCGTGa                                                                                                                                                                    <  2:277089/251‑1 (MQ=255)                                                                                              aTTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTACGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGaa                                                                                                                                                                   >  1:199881/1‑225 (MQ=255)                                                                                              aTTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTGTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGGTGCTTTTGATTTACGTTGGCGTCTGAGCATTGATTGATCCTCAAAAGCCCAGCCTCATTGGTAATGAACCAGCTCCGTGaa                                                                                                                                                                   <  2:199881/225‑1 (MQ=255)                                                                                                      cAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTACGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTCTTTGAtc                                                                                                                                 >  1:352524/1‑251 (MQ=255)                                                                                                       aaTCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGAtct                                                                                                                                <  1:151400/251‑1 (MQ=255)                                                                                                                  cAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGaa                                                                                                                      >  2:295565/1‑250 (MQ=255)                                                                                                                         ccTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAc                                                                                                              <  2:436452/251‑1 (MQ=255)                                                                                                                                   gTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGCTTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGt                                                                                                    >  2:145054/1‑251 (MQ=255)                                                                                                                                   gTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATTTCCGCTCCGGCCCAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGt                                                                                                    >  2:106012/1‑251 (MQ=255)                                                                                                                                      gCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATGCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCAc                                                                                                 <  1:392303/251‑1 (MQ=255)                                                                                                                                               gATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTTAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTgg                                                                                         <  2:35584/250‑1 (MQ=255)                                                                                                                                                     gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATAGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGc                                                                                                                                                            >  1:185641/1‑177 (MQ=255)                                                                                                                                                     gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTGTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGAGTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAAGGAACCAGCGCCGTGAAGGTCCGc                                                                                                                                                            <  2:185641/177‑1 (MQ=255)                                                                                                                                                       aaGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTaa                                                                                 >  2:305006/1‑250 (MQ=255)                                                                                                                                                        aGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTaaa                                                                                >  1:2569/1‑250 (MQ=255)                                                                                                                                                               cctacctaAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAACATTTTAAc                                                                        >  2:238418/1‑251 (MQ=255)                                                                                                                                                                ctacctaAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACt                                                                       >  1:170824/1‑251 (MQ=255)                                                                                                                                                                             agCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCg                                                           <  1:332370/250‑1 (MQ=255)                                                                                                                                                                                aaaCGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAGTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTTGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGc                                                       <  1:305006/251‑1 (MQ=255)                                                                                                                                                                                     cacaCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTa                                                   >  1:275576/1‑250 (MQ=255)                                                                                                                                                                                             ttAGAAAAATCGATGGTCAGAGACTGGATTAGCAATGCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTtat                                           <  2:82954/250‑1 (MQ=255)                                                                                                                                                                                               aGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATAc                                         >  2:328104/1‑250 (MQ=255)                                                                                                                                                                                                                          tAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCaaa              >  1:119629/1‑250 (MQ=255)                                                                                                                                                                                                                                    gTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCg   <  2:2569/251‑1 (MQ=255)                                                                                                                                                                                                                                     ttCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGt  >  1:169573/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      CCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGT  >  NC_000913/4510005‑4510483

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