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%	95.2 / NA	33	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 (14/19);  total (14/19)

GCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCG  >  NC_000913/4509992‑4510482                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     gCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGGAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAGCTGGCTGGTATAACGTGTCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGGAGCAAACGCGCACGTTTTAGAAAAAGCGATTGTCAGAAACTGGATTAGCAATTCCGTGCCAGGGt                                                                                                                                                                                                                                                   <  2:204289/250‑1 (MQ=255) gCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGt                                                                                                                                                                                                                                                   <  1:348565/250‑1 (MQ=255)  cTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGAAAGATGGGAGAGACACCGACTCATGAAGATATCATTGATAACAATAGCTCCATCAACTGGCCTGTATACCGTGGCTGCCGGTGCTTAAGAAAGCTGCACCTACCTAAGTCTTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTg                                                                                                                                                                                                                                                 >  1:292978/1‑251 (MQ=255)        gcTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCGTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTGTTTAGAAAAATCGATGTTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTg                                                                                                                                                                                                                                           <  2:470454/251‑1 (MQ=255)             cccTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGAttt                                                                                                                                                                                                                                       <  1:196802/250‑1 (MQ=255)              ccTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTAc                                                                                                                                                                                                                                     <  1:44182/251‑1 (MQ=255)                    gCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGTCTCATGAAGATATCATTGATAACAATCGAGCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTTCCTAAGTAGTAGCAAACGCACTCTTTTTATAAAAATCGATGGTCAGAAAGTCGATTAGCAAGTCCTTTCCAGGGGTGCTGTTGATGTACGTTGGc                                                                                                                                                                                                                               <  1:426795/251‑1 (MQ=255)                      gAAAAAGGGAGCAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGt                                                                                                                                                                                                                             >  1:374697/1‑251 (MQ=255)                                   ggACGGGCCCGGCTGCCTTCACTTTATGCCAGCTGCGGTAAGCTTGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTGAGAAACTGGATTAGCAATTCCGGTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAtt                                                                                                                                                                                                                <  1:447124/251‑1 (MQ=255)                                      cggacccGTCTGCCTTCACTTTATGCAAGCTACGGTAAGCTGGGTGAGACACCGAATCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTCGCTGCCGGTGATTAAGAAAGCTACACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTAt                                                                                                                                                                                                             >  2:198785/5‑251 (MQ=255)                                             gTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCCTTGATTTATCCTCaa                                                                                                                                                                                                       >  1:162550/1‑250 (MQ=255)                                               cTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTTAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCCTCACCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGAACCTACCTAGGTAGTAGCAAACGCGCACTTTTTAGAAAAATCGATGGTCAGAACTTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGTCGTCTGATCTTTGATTTATCCTaaaa                                                                                                                                                                                                      >  2:269504/1‑249 (MQ=255)                                                                          aaGCTGGGGGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGAAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCAGTGGTAATGAACCAGCt                                                                                                                                                                         <  1:354477/251‑1 (MQ=255)                                                                               gggTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTg                                                                                                                                                                    <  1:199489/251‑1 (MQ=255)                                                                                      acacCGTCTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAAGGAACCAGCTCCGTGAATGTcc                                                                                                                                                             <  2:480313/251‑1 (MQ=255)                                                                                                       tatCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGACGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGGTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCAAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCg                                                                                                                                             >  1:65650/1‑250 (MQ=255)                                                                                                                           ccATCAACTGGCCTGTATAACGTGGCTGTCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCTAACGCACTCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCGCTGGCACAGAGCGAAATTTTTTGATCTCCCACCt                                                                                                                        <  2:316211/251‑1 (MQ=255)                                                                                                                                                   gcttccGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGTCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCAc                                                                                                <  2:374697/247‑1 (MQ=255)                                                                                                                                                       ccGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGttt                                                                                            >  2:408547/1‑251 (MQ=255)                                                                                                                                                         ggTGATTCAGACAGATGCCCCTACCTAAGAAGTCGCAAGCGCCCACTCTTTAGATAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGCTCATTGTTTTATCCTCAACAGCCCAACCTCATTGGTAATCAACCAGCTCCGTGAATGTCCGCTCTGGCACCGAGCGAAATTTTTTGCTCTCCCTCCTGAAATCTAAACTTAGTCATGTCCCGttttt                                                                                          >  2:80297/1‑251 (MQ=255)                                                                                                                                                                                 ctaAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTAAGTTGGCGTCTGATCATTGCTTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGACACGTTTTTGGGTTTCTAAAATTTTAAATTcgc                                                                  >  2:416104/1‑251 (MQ=255)                                                                                                                                                                                        gtagCAAACGCACTCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTATCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGAGCATTTATTTATCCTCAAAAGCTCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGATTTTc                                                           <  2:363938/251‑1 (MQ=255)                                                                                                                                                                                            cAAACGCACACTTTTTAGAAAAATCGATGGTCCGAAACTGGATTCGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGACATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTCAAATTTTAACTTCGCGTTTTTCGTTg                                                       >  2:300458/1‑251 (MQ=255)                                                                                                                                                                                                      cTTTTTAGAAAAATCGATGGTCAGCAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCATCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACCGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAACCTTAGTCATGTCAAGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGt                                             >  2:502865/1‑251 (MQ=255)                                                                                                                                                                                                        ttttAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCATCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGtt                                            <  1:300458/250‑1 (MQ=255)                                                                                                                                                                                                            aGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAATCCCAACCTCAGTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATAc                                        <  2:119766/250‑1 (MQ=255)                                                                                                                                                                                                                  aTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATg                                 <  1:272584/251‑1 (MQ=255)                                                                                                                                                                                                                            aGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAACTTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGc                       >  1:10519/1‑251 (MQ=255)                                                                                                                                                                                                                            aGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGc                       <  1:291680/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTGGCCGTAAGGGTTATACAGAAATGGCCGTTAAGCAGAGTTCAAAatt          >  1:270812/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACGTCGCGTTTTTCGTTGCCGTAATGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCg  <  1:473899/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTAAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAATTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGtt                           >  2:200482/1‑226 (MQ=255)                                                                                                                                                                                                                                                 gTTCCAGGGTTGCGTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCTCGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGtt                           <  1:200482/226‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCG  >  NC_000913/4509992‑4510482

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