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	96.0%	72.0 / ‑4.1	25	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);  major base G (12/12);  minor base . (0/1);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCC  >  NC_000913/4509992‑4510481                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     gCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAATAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGt                                                                                                                                                                                                                                                   <  1:167524/250‑1 (MQ=255) gCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGt                                                                                                                                                                                                                                                   <  1:15313/250‑1 (MQ=255)        gcTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTg                                                                                                                                                                                                                                          <  1:140473/251‑1 (MQ=255)                         aaaGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCCGCCGGTGATTAAGCAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTg                                                                                                                                                                                                                         >  1:324/1‑251 (MQ=255)                                         gCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCt                                                                                                                                                                                                         <  1:187255/251‑1 (MQ=255)                                           ccGTCCGCCTTCACATTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCAACAACTGGGCTGTATAACATGGGAGCCGGTGAGTAAGAAAGCTGGATCTACCTAAGTAGTAGCAAACGGACACTTTGTAGAAAAATCGATGGGCAGAAGCTGGAGGAGCAATTCCGCTCGA‑GGTTGGCTTTTGATGCACGTTGGCGTCTGATCATTGATTTATCCTCa                                                                                                                                                                                                       <  1:156680/251‑1 (MQ=255)                                                       aCTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTGGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCt                                                                                                                                                                                           <  1:77904/251‑1 (MQ=255)                                                        cTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTc                                                                                                                                                                                          >  1:157612/1‑251 (MQ=255)                                                                       ggTAAGCTGGGTCAGACACCGCCTCATGAAGAAAGCATCGATAACAATCGATCCAACAACTGGCCTTTAAAACGTGGCAGCCGGTGATTAATAAAGATGCACCGACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGAGTAGCAAATCCGTTCCAGGGTTG‑CTTTTGATCTACGTTGTCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGATCCa                                                                                                                                                                           <  2:48434/251‑1 (MQ=255)                                                                                                                   cAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCAGTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGAtc                                                                                                                               <  2:216880/251‑1 (MQ=255)                                                                                                                           ccATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCt                                                                                                                       >  1:82765/1‑251 (MQ=255)                                                                                                                                aaCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAAt                                                                                                                  >  2:261988/1‑251 (MQ=255)                                                                                                                                          tataACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACtt                                                                                                          >  1:108460/1‑249 (MQ=255)                                                                                                                                                       ccGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTAt                                                                                                                                                                                                            >  1:162593/1‑138 (MQ=255)                                                                                                                                                       ccGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTAt                                                                                                                                                                                                            <  2:162593/138‑1 (MQ=255)                                                                                                                                                              ttAAGAAAGCTGCACCTACATAAATAGTAGCAAACCCACACTTTTTAGAAACATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTACGCTCTGGCACAGAGCGAACTTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTAATGTCACGTTTTTGCGtt                                                                                    >  2:237731/1‑251 (MQ=255)                                                                                                                                                                         gCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCCTTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGCATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAtttt                                                                         >  1:151190/1‑251 (MQ=255)                                                                                                                                                                                cctaAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcg                                                                  >  1:164445/1‑251 (MQ=255)                                                                                                                                                                                                    cacTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAgg                                              <  2:164445/251‑1 (MQ=255)                                                                                                                                                                                                            aGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGGTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCAGTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATAc                                       <  1:232839/250‑1 (MQ=255)                                                                                                                                                                                                                 aaTCGATGGTCAGCAACTGGATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATAATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGCATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATATAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAAt                                 >  2:250131/1‑251 (MQ=255)                                                                                                                                                                                                                                    gATTAGCAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTGCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCa              <  1:112931/251‑1 (MQ=255)                                                                                                                                                                                                                                          cAATTCCGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAattg        >  1:10558/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATGGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAattgattg    <  2:324/250‑1 (MQ=255)                                                                                                                                                                                                                                                cGTTCCAGGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACAAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGcc  >  2:250417/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTG‑CTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCC  >  NC_000913/4509992‑4510481

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