Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,605,702	A→C	intergenic (‑39/‑102)	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,702	0	A	C	94.1%	97.8 / ‑3.5	34	intergenic (‑39/‑102)	fhuF/yjjZ	ferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (16/16);  total (17/17)
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.43e-01

TAGCACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATC  >  NC_000913/4605466‑4605948                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       tAGCACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGttt                                                                                                                                                                                                                                           <  1:80228/250‑1 (MQ=255)  aGCACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTa                                                                                                                                                                                                                                          <  2:86254/250‑1 (MQ=255)            gTGATTACCAGTGCGGTAGCGTCATGGCATTAAGGGGGGCAGTTTCATCCAGGCGGATAAACGCCAGCAAAGGTGTACGATGGTTTGCGATCGTCGCCCGCACGGCTATCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCATCATAGAGCGGTGCGTAACGATAGGCCATAATCAGGATAGTAAGCTAAATGATAAGGATTGCTAATCCTAGCGATAGGTTTACCCGagcgcaa                                                                                                                                                                                                                               <  2:425160/251‑7 (MQ=255)               atgaCCATTGCGCCAGCGTCATGGCGTTAAGGGGGGCAGGTTCATCCAGGCGGATAAAATCAAGCAAATGTTGACGATGTTTTGCGCTCGTCGCCAGCCCGGCTTGCGCAACCGGGAGATCCTGCGGCTGGAGGTGCGCTCGCCAGATGACGTCTTCATCGAGCGGTGAGGCACGATAGGCCATAATCGGGGTAGTAATATAAATGATAAGGATTGCTAATCCTAGGGCTAGGTTTACCCGATAGCAAggg                                                                                                                                                                                                                            >  1:18369/1‑251 (MQ=255)                       tGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCt                                                                                                                                                                                                                    >  1:173910/1‑251 (MQ=255)                             aGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTg                                                                                                                                                                                                              <  1:176475/251‑1 (MQ=255)                               cGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGGTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCGAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCa                                                                                                                                                                                                            <  2:8291/251‑1 (MQ=255)                                     ggCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGAt                                                                                                                                                                                                      <  1:147396/251‑1 (MQ=255)                                         gtAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCGAAATGATAATGATGGCTAATCCTAGCGATAGGTTTACCAGATAGCAAGGGATTTATCTGGCTTGCAAATGATaaaa                                                                                                                                                                                                  <  1:181871/250‑1 (MQ=255)                                                   cAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGCTAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCatat                                                                                                                                                                                        >  1:501099/1‑251 (MQ=255)                                                   cAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCatat                                                                                                                                                                                        >  1:302575/1‑251 (MQ=255)                                                                                    gTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATg                                                                                                                                                        >  1:404438/1‑250 (MQ=255)                                                                                         cGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTACTCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATc                                                                                                                                                   >  2:245961/1‑250 (MQ=255)                                                                                                  gCGATCGTCGCCCTCACGGCTTTCGCAAGCGTGGGAGCCGGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAAGTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAAc                                                                                                                                         <  2:210789/251‑1 (MQ=255)                                                                                                     aTCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTGCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGTGATTTATCTGGCTTGCAAATGATAAAAATTATCTTATGTTATTTGTTATCCTTATCAATGTGAGAGATGTAATCATGTTGCACCGt                                                                                                                                       >  2:112615/1‑250 (MQ=255)                                                                                                       cgtcgCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTAc                                                                                                                                     <  1:139222/250‑1 (MQ=255)                                                                                                        gtcgCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGCTAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGc                                                                                                                                   >  2:80330/1‑251 (MQ=255)                                                                                                                    gCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCt                                                                                                                       >  2:30049/1‑251 (MQ=255)                                                                                                                        gcgcAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTggg                                                                                                                    <  1:245961/250‑1 (MQ=255)                                                                                                                            aaGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAgtg                                                                                                               <  1:398948/251‑1 (MQ=255)                                                                                                                                         tGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCATTGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCggggt                                                                                                  >  2:372292/1‑249 (MQ=255)                                                                                                                                           cGGCTGGGGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGCTTACCAGATAGCAAGTGATTTATCTGGCTTGCAAATGATAAAAATTTTCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTATGATGGGCAGTGGCTGGGGGGTGTTTCTGCCTGGATTg                                                                                                >  2:188358/1‑251 (MQ=255)                                                                                                                                                          tCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTg                                                                                 >  2:169419/1‑251 (MQ=255)                                                                                                                                                                           cATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATc                                                                >  1:433244/1‑251 (MQ=255)                                                                                                                                                                              agagCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGt                                                              <  1:372292/250‑1 (MQ=255)                                                                                                                                                                              agagCGGGGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATCTTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTc                                                             <  2:404438/251‑1 (MQ=255)                                                                                                                                                                                   ggttcgtaactattggCCATAATCGGGATAGTAAACTAAATGATAATGAATGCTAATCCTAGCGATAGGTTTTCCCGATAGCAAGGGATTTATCGGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCACTGGCTTGGGAGTGTTGCTGCCGGGATTGCTGATTGCATGGCTGATGTATGCGGATTTATCGTCAGATc                                                        <  2:272404/237‑1 (MQ=255)                                                                                                                                                                                              aTAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGAtt                                             <  1:224396/251‑1 (MQ=255)                                                                                                                                                                                                 ggCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCCAGGGATTTATCTGGCTTGCAAATGATCAACATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAGTCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGGTGCTGAATGCGGATTTATCGTCAGATCAGTGCGGGATTGTc                                          >  1:45611/1‑251 (MQ=255)                                                                                                                                                                                                            ggATAGTAATCTAGATGATAATGATTGCTAAGCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAGATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGGCGGGATGGCGGATGGCAGGGCTGATGGATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATggg                               <  2:84921/251‑1 (MQ=255)                                                                                                                                                                                                                     tCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTAtt                      >  1:217276/1‑251 (MQ=255)                                                                                                                                                                                                                       tAAATGAAAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGGTGGGATTAGTATTg                     >  1:343941/1‑250 (MQ=255)                                                                                                                                                                                                                                      gCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGt     >  1:427483/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaTCCTAGCGATAGGTTTAACCGATAGCAAGGGAGTTATCTGGCTTGCAAATGATAAAAAGTATCATATGATATTGGTTATCAGTATCAATGAGAGAGATGAAATCATGTTGAAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGTGATTAGTATTGACGCCGATGATGCTGTATc  <  2:263246/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       TAGCACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATC  >  NC_000913/4605466‑4605948

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