Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,568,153 A→T 11.4% intergenic (‑46/+171) eutH ← / ← eutG putative inner membrane protein/putative alcohol dehydrogenase in ethanolamine utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,568,1530AT11.4% 117.5 / 6.9 46intergenic (‑46/+171)eutH/eutGputative inner membrane protein/putative alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base A (18/22);  new base T (2/3);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GAGCAGAACCGCCGAACTGCGACAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCT  >  NC_000913/2568014‑2568291
                                                                                                                                           |                                                                                                                                          
gAGCAGAACCGCCGCACTGCCACAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCAGGATGTACATCATTATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAAc                                                                                                                                         <  1:578819/143‑1 (MQ=255)
caGCAGATCCGCTGAACTGCGACAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCTTCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAAc                                                                                                                                         <  2:601262/142‑1 (MQ=255)
 accAGAACCGCCGAACTGCGACAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTGACTCCCTGTGAGTTGTTATTTATTGGCGGATGCGGCGTAAACg                                                                                                                                        <  1:526990/141‑1 (MQ=255)
 aGCAGAACCGCCGAACTGCGACAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACg                                                                                                                                        <  2:287227/143‑1 (MQ=255)
      aaCCGCCGAACTGCGACAGGATCCTGTCTACGGCAGCTCTCAGCATAAAGAACATCATGACGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCCGTTAGTTTTTATTTATTCGCGGATGCGCCGTAACCGCCTTa                                                                                                                                   >  2:236261/1‑143 (MQ=255)
           ccgAACTGCGACAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATAAGGCGTTTACGCCGCATCCGc                                                                                                                              <  2:364887/143‑1 (MQ=255)
            cgAACTGCGACAGCATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATAAGGCGTTTACGCCGCATCCGcc                                                                                                                             <  1:408949/143‑1 (MQ=255)
              aaCTGCGACAGGATCCTGTCTACGGCAGCTCTCAGCGTAAAGAACATCATGATGTACATGATGATTTCGTTAATGCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATAAGGCGTTTACGCCGCATCCGccaa                                                                                                                           <  2:489502/143‑3 (MQ=255)
                    caCAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACAtgt                                                                                                                     <  1:218722/142‑1 (MQ=255)
                     aCAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACAtctg                                                                                                                    >  2:534362/1‑143 (MQ=255)
                               ttCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTa                                                                                                          <  1:279113/142‑1 (MQ=255)
                               gTCGACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTa                                                                                                          <  2:241555/143‑1 (MQ=255)
                                                gCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCACTACATGTGCAATCCCGTAGGCTGGATAAGACGCgg                                                                                         >  1:86630/1‑143 (MQ=255)
                                                 cATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGc                                                                                        >  1:311129/1‑143 (MQ=255)
                                                  aTAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCa                                                                                       >  2:211285/1‑143 (MQ=255)
                                                    aaaGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAg                                                                                     >  2:197128/1‑143 (MQ=255)
                                                       gAACAGCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGt                                                                                  <  1:104113/143‑1 (MQ=255)
                                                         aCATCATGATGTACATGATGATTTCGTTAATTCCCAGCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCg                                                                                <  2:473207/143‑1 (MQ=255)
                                                           atcatGATGTACATGATGATTTCGTTAATTCTCATCGCCTTTAATAAAAAATAAATGTAAATTATTGGCGGATGAGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTATGCTGGATAAGCCGCGGCAAGCGTcaca                                                                              >  1:138087/1‑143 (MQ=18)
                                                           atcatGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                                              >  2:123282/1‑143 (MQ=255)
                                                            tcatGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCAt                                                                             >  1:99954/1‑143 (MQ=255)
                                                              atgatgTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCATTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATcc                                                                           <  1:135480/143‑1 (MQ=255)
                                                                gatgTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCgg                                                                         >  1:627298/1‑143 (MQ=255)
                                                                   gTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCaa                                                                      <  1:603911/143‑1 (MQ=255)
                                                                   gTACATGATGATTTCGTTAATTCCCATCGCCGTTGCTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCaa                                                                      <  2:222524/143‑1 (MQ=255)
                                                                     aCATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAGACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAtt                                                                    <  1:578855/143‑1 (MQ=255)
                                                                             atTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCTTAGGCTGGATAAGACGCGGACAGAGTTGCATCCGGCAGTCGCATcgcg                                                            >  1:158083/1‑143 (MQ=25)
                                                                              tttCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCgcgc                                                           <  1:142816/143‑1 (MQ=255)
                                                                               ttCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTACACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGGCGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGcc                                                          >  1:551268/1‑143 (MQ=255)
                                                                                   ttAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCCCTg                                                      >  1:516203/1‑143 (MQ=255)
                                                                                            agcGCCTTTCCTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcg                                             <  1:444604/141‑1 (MQ=255)
                                                                                                 cTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGAGGGTGCGGCGTg                                        >  2:35224/1‑143 (MQ=255)
                                                                                                   ttACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGaa                                      >  2:230706/1‑143 (MQ=255)
                                                                                                       tCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGcc                                  <  2:615034/143‑1 (MQ=255)
                                                                                                            gTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATc                             <  1:197128/143‑1 (MQ=255)
                                                                                                               aGTGGTTATTTATTGGCGGATGCGGCGTAAACGCCTCACGCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGc                          <  1:610117/143‑1 (MQ=255)
                                                                                                                    ttatttatTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGGGAACGCCTTATCCGCCCTAc                     >  1:594534/1‑143 (MQ=255)
                                                                                                                        tcagtGGCTGATGCTGTGTAATCGCCTTGTCCGCCCGACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAtgt                 <  2:549875/139‑1 (MQ=17)
                                                                                                                          tagggCGGATAAGGCGTTTACGCCGCATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGCTCCGGCGCGAACGCCTCCTCCGCCCCCCCTGTGc               >  2:246532/4‑143 (MQ=17)
                                                                                                                          attGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATc                                                                            >  2:640042/1‑82 (MQ=255)
                                                                                                                          attGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGc               <  1:304109/143‑1 (MQ=255)
                                                                                                                            gggCGGATAAGGCGTTTACGCCGCATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGGGAACGCCGTATCCGCGCCACAGGTGCaa             >  1:635569/2‑143 (MQ=17)
                                                                                                                              gCGGATGCGGCGTATACGCCTTACCCGCCCTACATGTGCAATCCCGTAGGTTGGGTGCGGCTCCCCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATc           <  1:147330/143‑1 (MQ=17)
                                                                                                                                   ggcggcgTGAACGCCTTATCAGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTa      <  1:496387/142‑1 (MQ=255)
                                                                                                                                     cggcgTAAACGCCTTATCCGCCCTACATGTGCACTCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg    <  1:519982/143‑1 (MQ=255)
                                                                                                                                       gcgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGCGAAGGCCTTATGCGCCCTCCAGGTGCACTCCCGCAGGca  >  2:474824/1‑142 (MQ=18)
                                                                                                                                           |                                                                                                                                          
GAGCAGAACCGCCGAACTGCGACAGGATCCTGTCTACGGCAGCTATCAGCATAAAGAACATCATGATGTACATGATGATTTCGTTAATTCCCATCGCCTTTACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCT  >  NC_000913/2568014‑2568291

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