Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,605,791 A→G intergenic (‑128/‑13) 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_0009134,605,7910AG94.3% 92.7 / ‑2.6 35intergenic (‑128/‑13)fhuF/yjjZferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (20/13);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGG  >  NC_000913/4605549‑4606036
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
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)
                                                                                                                                                        cctAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATcac                                                                                         <  1:80330/248‑1 (MQ=255)
                                                                                                                                                                 aGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGc                                                                              >  1:234572/1‑251 (MQ=255)
                                                                                                                                                                                 aaGGGATTTATCTGGCTTGCAATTGATACAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGAAGTGGCTGGGGAGTGTTGCTGCCGGGATTGATGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGGTTGTACTTCTGAGGGGGTTCGTATTGAGACCGTTGGTGCTGTATCACAAAACATTGCGGCAATACATTTTGCt                                                              >  2:481009/1‑251 (MQ=255)
                                                                                                                                                                                 aaGGGATTTATCTGGCTTGCAAATGATAATAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTCTCACAAACAGTTGCGGCATTACATTTTGCt                                                              >  1:391711/1‑251 (MQ=255)
                                                                                                                                                                                                   gCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGATGCCGGGGTTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCCTTATGATGGGATTAGTATTTACCCCAATGATGCTGTATCACAAACAGTTGCGGCTATACATTTTGCTACCATCGTGCCTTGCACt                                            >  2:443098/1‑251 (MQ=255)
                                                                                                                                                                                                         gATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAAGCATGTTGCAACGTACGCTGGGCAGTGGCGGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTg                                       <  2:212321/250‑1 (MQ=255)
                                                                                                                                                                                                              aaaTTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGc                                  >  1:355036/1‑250 (MQ=255)
                                                                                                                                                                                                                   aTCATATGGTATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCCGATTGCAGGGATGATGTATGCGGATTTATCGTCAGATAAGTGGCGGATTGTCATTCTGATGGGGTTAGTATTGCCGCCGGAGAGGCTGTATCACACACAGATGCGGCATTGAATTTTGCTACCCTCGTGGCTGGCACTTTTTGCTGGcgcaatg                            >  2:84641/1‑245 (MQ=255)
                                                                                                                                                                                                                              ttGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCCTTCTGATGGGATTAGTATTGACGCCGATGATGCTGGATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGGGCCTGGGCCTTTTTGTTTGCATCATGCTTATGATaa                  >  2:311084/1‑250 (MQ=255)
                                                                                                                                                                                                                               tgggtagcataatcaAGGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGTCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGTATTAGTATTGACGCCTATGATGCTGTATCACAAACAGTTGCGGCATGACATTTTGCTACCATCGTGACTGGCACTGATTGCTGGCTTCATGCTGATGATAAt                 <  2:71047/239‑1 (MQ=255)
                                                                                                                                                                                                                                            caATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAg    <  1:169419/250‑1 (MQ=255)
                                                                                                                                                                                                                                             aaTGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCCTTCTGATGGGGTTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATAGTGCCAGGCCCTTTTTGCTGGCATCATGCTGTTGATAATGAATTTGAATCCggg  >  2:384687/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GTTCACGATGTTTTGCGATCGTCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGG  >  NC_000913/4605549‑4606036

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: