Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,566,464 A→G 7.4% S141P (TCT→CCT)  eutA ← reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,566,4640AG 7.4% 160.4 / 3.5 54S141P (TCT→CCT) eutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base A (32/18);  new base G (2/2);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

CCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGAT  >  NC_000913/2566324‑2566611
                                                                                                                                            |                                                                                                                                                   
cccGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGttt                                                                                                                                             >  1:102815/1‑149 (MQ=255)
                    gcagGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGaaa                                                                                                                                                   <  6:266138/123‑1 (MQ=255)
                    gcagGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGaaa                                                                                                                                                   >  5:266138/1‑123 (MQ=255)
                              accGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACg                                                                                                                >  4:19023/1‑148 (MQ=255)
                                   cGTTGCGGCCGGCAGTGACGCTGATTTTTACGGCATCGAACAGGGCGTAGTTCGCGGTGCAACCGCCGATTTAGCTATTCAGTAACCGACACAGCCGTTGTTAAGAAAGGGTTTGTGCCCCGGCTCCGGGACCGGCGATCACGGAGTc                                                                                                           >  4:19881/1‑148 (MQ=255)
                                    gTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGa                                                                                                         <  2:292960/149‑1 (MQ=255)
                                          gcaggcagTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCg                                                                                                   >  6:323270/1‑149 (MQ=255)
                                           caggcagTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCgg                                                                                                  >  8:93609/1‑149 (MQ=255)
                                            aggcagTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGAtt                                                                                                            <  3:203654/138‑1 (MQ=255)
                                            aggcagTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGAtt                                                                                                            >  4:203654/1‑138 (MQ=255)
                                                    ccGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCCATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCt                                                                                         >  4:47391/1‑149 (MQ=255)
                                                                 cGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                            >  5:90380/1‑149 (MQ=255)
                                                                 cGGCATCGAACAGGGCGTAGTTCGCGCTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTGGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaa                                                                             >  6:308538/1‑148 (MQ=255)
                                                                      tCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGcc                                                                       <  2:102815/149‑1 (MQ=255)
                                                                          aCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGc                                                                   >  7:380352/1‑149 (MQ=255)
                                                                          aCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGc                                                                   <  6:367651/149‑1 (MQ=255)
                                                                          aCAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGc                                                                   >  4:230834/1‑149 (MQ=255)
                                                                            aGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCTGCGCTGGCAACGACAAAATCGCCCAGCGa                                                                 >  2:54882/1‑149 (MQ=255)
                                                                              ggCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTTGAGGTGCGGCCCGGCGCTgg                                                                                       >  6:221068/1‑125 (MQ=255)
                                                                              ggCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTTGAGGTGCGGCCCGGCGCTgg                                                                                       <  5:221068/125‑1 (MQ=255)
                                                                                 gTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                            >  1:311998/1‑133 (MQ=255)
                                                                                 gTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                            <  2:311998/133‑1 (MQ=255)
                                                                                   aGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTga                                                             >  6:141324/1‑146 (MQ=255)
                                                                                        gcgGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgcgc                                                     >  1:43055/1‑149 (MQ=255)
                                                                                        gcgGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgcgc                                                     >  5:346123/1‑149 (MQ=255)
                                                                                        gcgGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAgcgc                                                     <  5:141324/149‑1 (MQ=255)
                                                                                               cACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAcc                                              <  3:47391/149‑1 (MQ=255)
                                                                                                aCCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCGGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGAGTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccg                                             >  6:148043/1‑149 (MQ=255)
                                                                                                 ccgccgATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCAccgc                                            >  1:23381/1‑149 (MQ=255)
                                                                                                  cgccgATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                            >  8:129984/1‑116 (MQ=255)
                                                                                                  cgccgATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACaaaa                                                                            <  7:129984/116‑1 (MQ=255)
                                                                                                    ccgATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGGTTGGGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCgg                                         <  3:19881/149‑1 (MQ=255)
                                                                                                      gaggTCGATATTCAGTACCCGACACATCCGTTGTTCAGAACGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGc                                       <  7:174403/146‑1 (MQ=255)
                                                                                                            gATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTgagaga                                                         <  4:204235/125‑1 (MQ=255)
                                                                                                            gATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTgagaga                                                         >  3:204235/1‑125 (MQ=255)
                                                                                                            gATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCAt                                 >  4:170795/1‑149 (MQ=255)
                                                                                                                   aGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGt                          >  8:356177/1‑149 (MQ=255)
                                                                                                                       cccGACACAGCCGTTGTTCAGGAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGccc                                                                                               >  5:296837/1‑76 (MQ=255)
                                                                                                                       cccGACACAGCCGTTGTTCAGGAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGccc                                                                                               <  6:302101/76‑1 (MQ=255)
                                                                                                                       cccGACACAGCCGTTGTTCAGGAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGccc                                                                                               <  6:296837/76‑1 (MQ=255)
                                                                                                                       cccGACACAGCCGTTGTTCAGGAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGccc                                                                                               >  5:302101/1‑76 (MQ=255)
                                                                                                                        ccGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTcg                     >  7:374089/1‑149 (MQ=255)
                                                                                                                         cGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTcgc                    <  3:19023/149‑1 (MQ=255)
                                                                                                                          gACACAGCCGTTGTTCAGAAAGGGTTTGTGACCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTcgcg                   >  8:192964/1‑149 (MQ=255)
                                                                                                                           acacAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCgcgc                  <  6:90380/149‑1 (MQ=255)
                                                                                                                               aGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCttt               >  2:257511/1‑148 (MQ=255)
                                                                                                                               aGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAACCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCtat               >  8:309534/1‑146 (MQ=255)
                                                                                                                                  cGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGCCCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCAcc           >  8:141760/1‑149 (MQ=255)
                                                                                                                                    ttgttCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCgg         >  5:39498/1‑149 (MQ=255)
                                                                                                                                          aGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGATGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatg    <  5:308538/148‑1 (MQ=255)
                                                                                                                                           gAAAGGGTTTGTGCCCCGGCTCCGTGCCCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGCGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgat  >  4:152372/1‑149 (MQ=255)
                                                                                                                                           gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCgg                                         <  4:53821/110‑1 (MQ=255)
                                                                                                                                           gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCAc            <  2:47242/139‑1 (MQ=255)
                                                                                                                                           gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCAc            >  1:47242/1‑139 (MQ=255)
                                                                                                                                           gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgat  >  3:101477/1‑149 (MQ=255)
                                                                                                                                           gAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgat  <  1:54882/149‑1 (MQ=255)
                                                                                                                                            aaaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgat  >  5:336521/1‑148 (MQ=255)
                                                                                                                                            aaaGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGtgatgat  >  5:360168/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CCCGTGGCTGTCGGTTTCCAGCAGGCGACCACCGACGTTGAGGCAGGCAGTGCCGCTGATTTTTCCGGCATCGAACAGGGCGTAGTTCGCGGTGCCACCGCCGATGTCGATATTCAGTACCCGACACAGCCGTTGTTCAGAAAGGGTTTGTGCCCCGGCTCCGTGACCGGCGATCACGGATTCGAGGTGCGGCCCGGCGCTGGCAACGACAAAATCGCCCAGCGACTGAGAGAGCGCCATCACCGCCGGGCGAGCATTGCGGGTTTTCGCGCTTTCACCGGTGATGAT  >  NC_000913/2566324‑2566611

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