Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 613,699 A→T 8.3% L321F (TTA→TTT)  fes → enterobactin/ferrienterobactin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913613,6990AT 8.3% 123.3 / 2.6 48L321F (TTA→TTT) fesenterobactin/ferrienterobactin esterase
Reads supporting (aligned to +/- strand):  ref base A (24/20);  new base T (2/2);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATGCGCAATTACACCCCATAAAAGAATCCATTTTCTGGC  >  NC_000913/613576‑613841
                                                                                                                           |                                                                                                                                              
gagTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAg                                                                                                                       <  8:203158/149‑1 (MQ=255)
        gtggGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGcc                                                                                                               >  3:3955/1‑149 (MQ=255)
         tggGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGaaaagct                                                                                                                                     <  4:93031/126‑1 (MQ=255)
         tggGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGaaaagct                                                                                                                                     >  3:93031/1‑126 (MQ=255)
         tggGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGcc                                                                                                               >  7:45063/1‑148 (MQ=255)
          gggCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCg                                                                                                              >  8:270410/1‑148 (MQ=255)
          gggCTTTCCGCGCTGTATGCCGGACTGCACTGGCCGGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCg                                                                                                              >  4:289718/1‑148 (MQ=255)
            gCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTTCTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAg                                                                                                           >  1:65368/1‑149 (MQ=255)
            gCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAg                                                                                                           <  4:3955/149‑1 (MQ=255)
              tttCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTTCTTGAAAAGCTAAAAGCTg                                                                                                                             <  6:359622/129‑1 (MQ=255)
              tttCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTTCTTGAAAAGCTAAAAGCTg                                                                                                                             >  5:359622/1‑129 (MQ=255)
              tttCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGt                                                                                                         >  7:29517/1‑149 (MQ=255)
                 ccGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCt                                                                                                       <  5:86228/148‑1 (MQ=255)
                  cgcgCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTgc                                                                                                     <  4:144118/149‑1 (MQ=255)
                                     cACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCggg                                                                                  >  8:17508/1‑149 (MQ=255)
                                            cTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTTCTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTcgc                                                                           <  2:65368/149‑1 (MQ=255)
                                             tGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAg                                                                                                                       >  8:120367/1‑104 (MQ=255)
                                             tGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAg                                                                                                                       <  7:120367/104‑1 (MQ=255)
                                                  gCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGcc                                                                      >  5:142422/1‑148 (MQ=255)
                                                  gCTTTGGCTGTGGATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAACAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCCGTGTGCTGGAAGCGGGTATTCGCGCGcc                                                                      >  4:306591/1‑148 (MQ=255)
                                                      tGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCgatg                                                                  <  1:207523/148‑1 (MQ=255)
                                                        gCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCgatgat                                                                >  4:19242/1‑148 (MQ=255)
                                                         cTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCa                                                              <  6:142422/149‑1 (MQ=255)
                                                          tgtgTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGCGGGCGTGTTACTAGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGCTGATCAt                                                             >  4:308836/1‑149 (MQ=255)
                                                          tgtgTATTAAGCCAGTCAGGATCGTACTGGTGACCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCAt                                                             >  6:68675/1‑149 (MQ=255)
                                                              ttttAAGCCAGTCTGGTTCGTTGTGGTGGCCGCTTCGGGGCGGGCAGCAAGTGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCg                                                          <  2:264163/146‑1 (MQ=255)
                                                              tATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCg                                                          >  1:264163/1‑148 (MQ=255)
                                                                 tAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGc                                                       >  2:143969/1‑148 (MQ=255)
                                                                    gacAGTCAGGATCGTACTCGTGGCCGCATCGGGGCGGGCAGCAACAGGGCTTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCAATGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAAt                                                   >  2:28220/3‑149 (MQ=255)
                                                                     ccAGTCAGGATGGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTAGTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGAATTGTTCTGGAAGCGTGTATTATCGAGCCGATGATCATg                                                            <  1:58889/139‑1 (MQ=255)
                                                                     ccAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATg                                                            >  2:58889/1‑139 (MQ=255)
                                                                           aGGATCGTACTGGTGACCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGgcgc                                            <  5:68675/149‑1 (MQ=255)
                                                                             gATCGTACTGGTGGTCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTg                                          >  5:34127/1‑149 (MQ=255)
                                                                               tCGTAGTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTATATGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCTATTATCATGCGGGCCAATGAGGCGCTGTa                                        <  1:28220/149‑1 (MQ=255)
                                                                                  tACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCgatga                                                                 <  2:111269/121‑1 (MQ=255)
                                                                                  tACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCgatga                                                                 >  1:111269/1‑121 (MQ=255)
                                                                                    cTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATgcg                                    >  6:247474/1‑148 (MQ=255)
                                                                                       gtggCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGGTAGCGCCGAAGGTCTGCGCATTGTGCTGGCAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATGCGCa                                  <  3:289718/147‑1 (MQ=255)
                                                                                            cGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATGCGCAATTacac                           >  4:254168/1‑149 (MQ=255)
                                                                                                tCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATGCGCAATTACACCCCa                       <  3:19242/149‑1 (MQ=255)
                                                                                                           agcaAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGTCAATCAGGCGCTGTATGCGCAATTACAcccc                        >  6:214778/1‑137 (MQ=255)
                                                                                                           agcaAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGTCAATCAGGCGCTGTATGCGCAATTACAcccc                        <  5:214778/137‑1 (MQ=255)
                                                                                                           agcaAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGAGCATTGTGCTGGAAGCGGGTATTCTCGAGCCTATGATCCTGCGGGCCAATCAGGCGCTGTATTCTCAATTACACACCAAAAAAGAATc             >  6:25397/1‑148 (MQ=255)
                                                                                                             caAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATGCGCAATTACACCCCATAAAAGAATCCAt          <  1:143969/149‑1 (MQ=255)
                                                                                                                  gggCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCgatga                                                                 <  1:112192/89‑1 (MQ=255)
                                                                                                                  gggCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCgatga                                                                 >  2:112192/1‑89 (MQ=255)
                                                                                                                   ggCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATGCGCAATTACACCCCATAAAAGAATCCATTTTCTg    <  7:17508/149‑1 (MQ=255)
                                                                                                                     cGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGCAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATGCGCAATTACACCCCATAAAAGAATCCATTTTCTGGc  <  8:29517/149‑1 (MQ=255)
                                                                                                                           |                                                                                                                                              
GAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTACTTGAAAAGCTAAAAGCTGGTGAAGTTAGCGCCGAAGGTCTGCGCATTGTGCTGGAAGCGGGTATTCGCGAGCCGATGATCATGCGGGCCAATCAGGCGCTGTATGCGCAATTACACCCCATAAAAGAATCCATTTTCTGGC  >  NC_000913/613576‑613841

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