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	88.6%	94.9 / 3.0	35	intergenic (‑39/‑102)	fhuF/yjjZ	ferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (18/13);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGT  >  NC_000913/4605469‑4605945                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    cacaTTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTAccc                                                                                                                                                                                                                                    >  1:252208/1‑251 (MQ=255)   caTTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATAGCTAATCATAGCGATAGGTTTACCCg                                                                                                                                                                                                                                   >  2:329180/1‑250 (MQ=255)                          aGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTg                                                                                                                                                                                                           <  1:68839/251‑1 (MQ=255)                                            ggggCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCTCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTAt                                                                                                                                                                                          <  1:153583/250‑1 (MQ=255)                                            ggggCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATc                                                                                                                                                                                         <  1:246914/251‑1 (MQ=255)                                              ggCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATACAAATTATCat                                                                                                                                                                                       >  1:138422/1‑251 (MQ=255)                                                       aTCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTGTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTg                                                                                                                                                                              <  2:241621/251‑1 (MQ=255)                                                              cGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGttatc                                                                                                                                                                        >  1:25953/1‑250 (MQ=255)                                                                 aTAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGttatcatta                                                                                                                                                                    <  1:180494/251‑1 (MQ=255)                                                                       tCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGCTGACATCTTCATAGAGCGGTGCGGCACGATAGGCCATAATCGGGATAGTAAACTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGATTGCAAATGATAAAAATTATCATATGGTATTGGTTATCATTATCAATg                                                                                                                                                              >  1:287870/1‑251 (MQ=255)                                                                       tCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATg                                                                                                                                                              >  2:10418/1‑251 (MQ=255)                                                                        ccAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGa                                                                                                                                                             <  1:244451/251‑1 (MQ=255)                                                                                       gATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATCGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGttatc                                                                                                                                                                        >  2:465056/1‑225 (MQ=255)                                                                                       gATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGttatc                                                                                                                                                                        <  1:465056/225‑1 (MQ=255)                                                                                                    cgtcgCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGa                                                                                                                                                             <  1:248889/223‑1 (MQ=255)                                                                                                    cgtcgCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGa                                                                                                                                                             >  2:248889/1‑223 (MQ=255)                                                                                                      tcgCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCt                                                                                                                               >  1:76124/1‑251 (MQ=255)                                                                                                          ccGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGc                                                                                                                           >  1:299925/1‑251 (MQ=255)                                                                                                                                             ggAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTa                                                                                                                                                                                           <  2:232206/152‑1 (MQ=255)                                                                                                                                             ggAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTa                                                                                                                                                                                           >  1:232206/1‑152 (MQ=255)                                                                                                                                                        cGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTg                                                                              <  1:478187/250‑1 (MQ=255)                                                                                                                                                         gCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATAATTATCAATGAGAGAGATGAAATCATGTTGCAACGTAAGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGAt                                                                            >  1:194560/1‑251 (MQ=255)                                                                                                                                                                gACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCg                                                                     <  1:417604/251‑1 (MQ=255)                                                                                                                                                                   aTCTTCATAGAGCGGTGCGGAACGCTAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGGGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGGTGGATGCGGAt                                                                  >  2:283049/1‑251 (MQ=255)                                                                                                                                                                     cTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGAttt                                                                >  1:60767/1‑251 (MQ=255)                                                                                                                                                                               cggggCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTc                                                          >  2:70257/5‑247 (MQ=255)                                                                                                                                                                                                    aaTCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGGTGTATGCGGATTTATCGTCAGATCAGTGGCGGGTTGTCATTCtt                                 >  2:450046/1‑250 (MQ=255)                                                                                                                                                                                                         ggATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGGTTgggatcggaagagg                            >  2:163162/1‑238 (MQ=255)                                                                                                                                                                                                         ggATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTg                                         <  1:163162/238‑1 (MQ=255)                                                                                                                                                                                                            tAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAACATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGCTTTATCGTCAGATCAGTGGCGGATTGTCCTTCTGGTggggtt                         >  2:130377/1‑248 (MQ=255)                                                                                                                                                                                                               tAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTACTTCTTATGGGGTTAGt                      >  2:284049/1‑251 (MQ=255)                                                                                                                                                                                                                        tgataatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGGTGGGGTTAGTATTGGCGcc             >  1:332461/1‑251 (MQ=255)                                                                                                                                                                                                                        tgataatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGcc             >  1:55007/1‑251 (MQ=255)                                                                                                                                                                                                                         gataatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCg            <  1:450046/251‑1 (MQ=255)                                                                                                                                                                                                                                   gCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGt  <  1:418404/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    CACATTCGGTGATGACCATTGCGCCAGCGTCATGGCGTTAAGTGGGGCAGGTTCATCCAGGCGGATAAACTCCAGCAAATGTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGT  >  NC_000913/4605469‑4605945

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