Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,442 Δ1 bp 15.1% intergenic (+14/+65) ftnB → / ← yecJ ferritin B, putative ferrous iron reservoir/DUF2766 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,4420T.15.1% 204.9 / 33.5 73intergenic (+14/+65)ftnB/yecJferritin B, putative ferrous iron reservoir/DUF2766 family protein
Reads supporting (aligned to +/- strand):  ref base T (34/28);  new base . (5/6);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

AAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTT  >  NC_000913/1987343‑1987541
                                                                                                   |                                                                                                    
aaGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCt                                                                                                      <  1:2131377/100‑1 (MQ=255)
 aGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTa                                                                                                     >  1:926520/1‑100 (MQ=255)
  gTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                    >  1:1094646/1‑100 (MQ=255)
  gTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTaa                                                                                                    <  1:286903/100‑1 (MQ=255)
     ctcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGc                                                                                                 >  1:1613023/3‑100 (MQ=255)
         gTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAtt                                                                                              >  1:1783470/1‑99 (MQ=255)
         gTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCAttag                                                                                            <  1:1818316/100‑3 (MQ=255)
              aaaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCg                                                                                       <  2:360797/100‑1 (MQ=255)
              aaaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTAACACCAGCTGCATTGATCATCATATGCGCTAATGCATT‑tcgc                                                                                         >  2:229744/1‑98 (MQ=255)
               aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGa                                                                                      <  2:1016284/100‑1 (MQ=39)
               aaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCGGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGa                                                                                      <  2:1700142/100‑1 (MQ=39)
                aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGAt                                                                                     <  1:733418/100‑1 (MQ=255)
                aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGC‑AATGCATTAGCGCCGAt                                                                                     <  1:1555669/100‑1 (MQ=39)
                aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCAATGATCATCATCGGCGC‑AATGCATTAGCGCCGAt                                                                                     <  2:2262770/100‑1 (MQ=38)
                 cTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATg                                                                                    >  2:1961514/1‑100 (MQ=255)
                 cTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATg                                                                                    >  2:292065/1‑100 (MQ=255)
                     cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAgg                                                                                >  1:386730/1‑100 (MQ=255)
                     cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAgg                                                                                >  2:1318362/1‑100 (MQ=255)
                     cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAgg                                                                                <  1:813652/100‑1 (MQ=255)
                           tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTga                                                                          <  1:1421966/100‑1 (MQ=255)
                           tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGATTTga                                                                          >  1:900650/1‑100 (MQ=255)
                           tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGATTTga                                                                          >  1:61116/1‑100 (MQ=255)
                           tgtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑ACGCCGATGAAAGTTTTga                                                                          >  1:809454/1‑100 (MQ=255)
                            gtgCTCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTgag                                                                         <  1:1016784/100‑1 (MQ=255)
                            gtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTTGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTgag                                                                         <  1:1688036/100‑1 (MQ=255)
                            gtgCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTgag                                                                         <  1:2169547/100‑1 (MQ=255)
                              gCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaa                                                                       >  1:312971/1‑100 (MQ=255)
                              gCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGaa                                                                       >  1:1098460/1‑100 (MQ=255)
                                ccTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAAc                                                                     <  1:1696754/100‑1 (MQ=255)
                                   gtgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATAATCGGCGCTAATGCATT‑GCACCGATGAAGGTATTGAGAAACCg                                                                   >  1:904566/1‑99 (MQ=255)
                                        gaccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCt                                                             <  2:1947851/100‑1 (MQ=255)
                                          ccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCa                                                           >  1:344572/1‑100 (MQ=255)
                                            gaccAACATGGTCTGTATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGCGCAGCCGCTGCctcctc                                                         >  2:1639910/1‑100 (MQ=255)
                                              ccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATc                                                         >  2:1346712/1‑98 (MQ=255)
                                                aaCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCAGCATCGGCGCTAATGCATT‑GAGCCGATGAAGGTTTTGAGAAACCGCTGCCGCATCTGt                                                      >  2:1303662/1‑99 (MQ=255)
                                                  cATGTTCTGAATGTCGTGTCACACCCGCTGCATTGCTCATCATCGGCGCTCATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTg                                                   <  2:705563/100‑1 (MQ=255)
                                                  cATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTg                                                   <  2:31905/100‑1 (MQ=255)
                                                  cATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTg                                                   >  2:379651/1‑100 (MQ=255)
                                                    tGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCTTCaaacagatgag                                                 <  2:386727/100‑11 (MQ=255)
                                                    tGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGaa                                                 >  1:84135/1‑100 (MQ=255)
                                                         tGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAagcagc                                            <  1:110615/100‑1 (MQ=255)
                                                            aTGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGt                                         >  1:1384379/1‑100 (MQ=255)
                                                              gTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttt                                       >  1:2136055/1‑100 (MQ=255)
                                                               tCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttt                                      <  1:306756/100‑1 (MQ=255)
                                                                 gtgtCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttttt                                    <  2:867889/100‑1 (MQ=255)
                                                                  tgtCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGttttttt                                   <  1:993519/100‑1 (MQ=255)
                                                                     cacaCCAGGTGCATTGATAATCATCGGAGCTAATGCATT‑GCGTCGATGAAGGTATTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGGTGTTTTAAt                                <  2:34622/100‑1 (MQ=255)
                                                                     cacaCCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                >  1:2398213/1‑100 (MQ=255)
                                                                     cacaCCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                >  1:2398061/1‑100 (MQ=255)
                                                                      acacCAGCTGCATTGATCATCATCGGCGCTAATGCGTT‑GCGTCGATGAAGGTTTTGTGAATTCGTTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                <  2:499339/99‑1 (MQ=255)
                                                                        acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATggg                             <  2:413678/100‑1 (MQ=255)
                                                                        acCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATgg                              <  1:610770/99‑1 (MQ=255)
                                                                         ccAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGTGa                            >  1:1945966/1‑100 (MQ=255)
                                                                         ccAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGa                            >  1:1685702/1‑100 (MQ=255)
                                                                          cAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAt                           >  1:1391980/1‑100 (MQ=255)
                                                                           aGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGAtt                          <  1:220691/100‑1 (MQ=255)
                                                                            gCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTc                         >  2:2316062/1‑100 (MQ=255)
                                                                               gCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc                      >  2:553649/1‑100 (MQ=255)
                                                                               ctcaaaaccttcatcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAcc                      >  1:1724361/11‑100 (MQ=25)
                                                                                cATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCAccc                     <  2:210714/100‑1 (MQ=255)
                                                                                    gATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                 >  1:1251392/1‑100 (MQ=255)
                                                                                        atcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTTTTTGAAGCAGCGGTTTTTTTAATGGGATACACCCTGTgggct             <  1:1622513/100‑3 (MQ=255)
                                                                                        atcatcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGt             <  2:1603446/100‑1 (MQ=255)
                                                                                          catcGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTaa           <  2:639244/100‑1 (MQ=255)
                                                                                          catcGGCGC‑AATGCATTAGCGCCGCTGAAGGTTCTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTa            >  1:2322685/1‑99 (MQ=38)
                                                                                           atcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTaaa          >  1:2260168/1‑100 (MQ=39)
                                                                                           atcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTaaa          >  1:790350/1‑100 (MQ=39)
                                                                                           atcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTaaa          >  2:493136/1‑100 (MQ=39)
                                                                                            tcGGCGC‑AATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAAc         <  2:1724340/100‑1 (MQ=39)
                                                                                              ggCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGAGTAAAc         >  1:18527/1‑98 (MQ=255)
                                                                                                cgcTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGa     >  1:2038869/1‑100 (MQ=255)
                                                                                                 gcTAATTCATT‑GTGGCAATGAAGATTTTGAGAAACCGCTGCCTCATCTGTTAGAAGAAGCGTGTTATTTAATGGGTGTCAGCCGCTGGGGTAAAc         <  2:1597504/95‑1 (MQ=255)
                                                                                                 gcTAATGCATT‑GCGCCGATGAAGGTTTTGTGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg    <  2:1444455/100‑1 (MQ=255)
                                                                                                  cTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAg    >  1:412531/1‑99 (MQ=255)
                                                                                                   tAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTTATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGt   <  2:409071/99‑1 (MQ=255)
                                                                                                   tAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt  >  1:1484912/1‑100 (MQ=255)
                                                                                                   tAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt  <  1:118671/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
AAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATT‑GCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTT  >  NC_000913/1987343‑1987541

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