Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,450:1 +A 8.0% intergenic (+22/+57) ftnB → / ← yecJ putative ferritin‑like protein/DUF2766 domain‑containing protein YecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,4501.A 8.0% 181.2 / 11.3 50intergenic (+22/+57)ftnB/yecJputative ferritin‑like protein/DUF2766 domain‑containing protein YecJ
Reads supporting (aligned to +/- strand):  ref base . (23/23);  new base A (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

GATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGG  >  NC_000913/1987312‑1987592
                                                                                                                                           |                                                                                                                                              
gatgGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGACACGTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCAttagcgc                                                                                                                                            <  1:232388/143‑6 (MQ=255)
     tCTGTTGGTGCAAACCAGTCTTGATGAAGTGCGGAGTGCGATACTTGCGGGTATGTGCACTGTGAAGCCCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCATTAGCGCCGATg                                                                                                                                       <  2:53822/143‑1 (MQ=255)
      cggtggctgcAAACCATTCTTGATGTAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCAATGCATTAGCGCCGATGa                                                                                                                                      <  2:519271/138‑1 (MQ=255)
       tGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGaa                                                                                                                                     <  2:490076/143‑1 (MQ=255)
        gTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGACGAAg                                                                                                                                    >  2:516794/1‑143 (MQ=255)
             tgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGCAGGtttt                                                                                                                               >  1:733963/1‑143 (MQ=255)
             tgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGtttt                                                                                                                               <  2:98066/143‑1 (MQ=255)
               cAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTga                                                                                                                             <  2:20411/143‑1 (MQ=255)
               cAAACCATTCTTGATGAAGCGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTga                                                                                                                             <  1:235728/143‑1 (MQ=255)
                     aTTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAAcc                                                                                                                       >  2:333902/1‑143 (MQ=255)
                      ttCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCg                                                                                                                      >  2:309045/1‑143 (MQ=255)
                        cTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCt                                                                                                                    <  2:410154/143‑1 (MQ=255)
                               aaGTGTGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCAt                                                                                                             <  1:198732/143‑1 (MQ=255)
                                    cgcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACGGCGCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAAGCGCTGCCTCATCGGtt                                                                                                        >  2:442274/1‑143 (MQ=255)
                                      cAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGCGACCCCGCTGCCTCATCTGTTTg                                                                                                      >  2:126354/1‑143 (MQ=255)
                                       aGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGa                                                                                                     <  1:352035/143‑1 (MQ=255)
                                           cGAAACTTGCGGGTATGTGCCCTGTGCAGACAGACCAACATGTACTGAATGTCGTGTCACACCAGCTGGATGGATCATCATCGGCGCTGATGCTTTGCGTCGCTGAAAGGTCTGAGCAAACGCAACCTCATCTGTTTGAagca                                                                                                 >  2:46175/1‑143 (MQ=255)
                                               aCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCgg                                                                                             <  1:593325/143‑1 (MQ=255)
                                                   gCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCCTTGATCAGCAGCGGCGCTAATGCATTGCGCCGATGAAGGGTTTGAGAAACCGCTGCCTCACCTGTTTGAAGCAGCGGtttt                                                                                         >  2:744510/1‑143 (MQ=255)
                                                   gCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGtttt                                                                                         <  1:651266/143‑1 (MQ=255)
                                                    cGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCGTCTGTTTGAAGCAGCGGttttt                                                                                        >  2:616855/1‑143 (MQ=255)
                                                         aTGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAAt                                                                                   <  1:333902/143‑1 (MQ=255)
                                                                      aGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtg                                                                      >  1:633797/1‑143 (MQ=255)
                                                                        actgaccAACACGTTCTTGATGTCGTGTCACACCAGTTGCTGGGCTCATCATCGGCGCTAATGCATTGCTCCGATGAAGGTTTTGAGAAACCGCTGCCTCAGCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCtgtg                                                                    <  2:320217/140‑1 (MQ=255)
                                                                           gaccAACATGTTCTGAATGGCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTgggg                                                                 <  1:83649/143‑1 (MQ=255)
                                                                               aaCATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGAATTCCCCCTGTGGGGTaca                                                             >  1:556982/1‑141 (MQ=255)
                                                                                      tCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGGAGCAGCGGTTTTTTTAATGGGATTCCCCCTGTGGGGTAAGCGTGAGt                                                      >  1:197561/1‑143 (MQ=255)
                                                                                       cTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGtt                                                     <  2:238198/143‑1 (MQ=255)
                                                                                        tGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTc                                                    <  1:233909/143‑1 (MQ=255)
                                                                                         gAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCa                                                   <  1:445924/143‑1 (MQ=255)
                                                                                              tCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTGTTTTGAGGGGAATCACCCAGGGGGGGGAGCTGGAGTTGaaggag                                              >  1:256189/1‑139 (MQ=255)
                                                                                               cGTGTCACACCAGCTGCATTGCTTATTATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAgc                                             <  2:128221/143‑1 (MQ=255)
                                                                                               cGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTGTTAATGGGATTCACCCTGTGGGGTAAACGCGAGGTCAATAAgc                                             >  2:30248/1‑143 (MQ=255)
                                                                                                    cacaGCAGCTGCATTGATCATCATCGGCGCTAATGCATGGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGAt                                        <  1:207196/143‑1 (MQ=255)
                                                                                                      cacCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGGTTCACCCTGTGGGGTCAACTTGAGTTCAATGAGCGCGATgg                                      >  1:132628/1‑143 (MQ=255)
                                                                                                           gCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCtttt                                 <  1:126354/143‑1 (MQ=255)
                                                                                                           gCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCtttt                                 >  1:794487/1‑143 (MQ=255)
                                                                                                               cATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTCTATCATCGGCTTCTCCCGGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGa                             <  1:766439/143‑1 (MQ=255)
                                                                                                                aTTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCGGTGTGAAGCAGCGGGTTTTTTAATGGGATTCACCCTGTGGGGTAAACGTGAGTTCAATAAGCGCGATGGCTTTTGGGAt                            >  2:643043/1‑143 (MQ=255)
                                                                                                                    atcatcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTGGAGTTCAATAAGCGCGCGGGCTTTTTGGAtgcca                        >  2:84437/1‑139 (MQ=255)
                                                                                                                     tcatcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGccc                       >  2:454197/1‑143 (MQ=255)
                                                                                                                      cttcatcGGCGCAATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCg                      >  1:39625/3‑143 (MQ=255)
                                                                                                                       atcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTAACTAAAGGAGCGGCTTTTTTATTGGGATTCACCCTGTGGGGCAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGc                     >  1:373596/1‑143 (MQ=255)
                                                                                                                         catcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCAt                   <  1:406344/143‑1 (MQ=255)
                                                                                                                          atcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGCAGCAGCGGTGTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAGGCGCGATGGCTTTTTGGCTTGCCCGCGTg                  >  2:552731/1‑143 (MQ=255)
                                                                                                                          atcGGCGCAATGCATTAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATg                  >  1:733310/1‑143 (MQ=255)
                                                                                                                            cGGCGCTAATGCATTGCGCGGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGt                <  2:384307/143‑1 (MQ=255)
                                                                                                                             ggCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg               >  1:523831/1‑143 (MQ=255)
                                                                                                                             ggCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTg               >  1:202125/1‑143 (MQ=255)
                                                                                                                                gccaaTGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTATGCTGTGATTGCCGAAGGAAAAATGAAGTTTTTGCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGAcc            <  2:169777/140‑1 (MQ=255)
                                                                                                                                    aTGCTTTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGt        <  2:792755/143‑1 (MQ=255)
                                                                                                                                         gtGCGCCGATGAAGGTGTTGGGAAGCCGCGGCCTCGTCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAGAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcg   <  2:625442/142‑1 (MQ=255)
                                                                                                                                          tGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcgg  <  1:653651/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGG  >  NC_000913/1987312‑1987592

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