Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 628,905 T→C 8.1% A431A (GCA→GCG)  eutB ← ethanolamine ammonia‑lyase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947628,9050TC 8.1% 138.8 / 2.7 49A431A (GCA→GCG) eutBethanolamine ammonia‑lyase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (25/20);  new base C (2/2);  total (27/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 = 9.86e-01

CCAAGGGTTGTCGGGGGTAGGCGTGCGGTGGTCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATGAAGTTGATGCCGGCCACCCCCAGCAGCGT  >  NC_002947/628759‑629051
                                                                                                                                                   |                                                                                                                                                  
ccaAGGGTTGTCGGGGGTAGGCGTGCGGTGGTCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGc                                                                                                                                                  <  3:9414/149‑1 (MQ=255)
ccaAGGGTTGTCGGGGGTAGGCGTGCGGTGGTCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGc                                                                                                                                                  >  6:50289/1‑149 (MQ=255)
  aaGGGTTGTCGGGGGTAGGCGTGCGGTGGTCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCt                                                                                                                                                >  1:106382/1‑149 (MQ=255)
          tCGGGGGTAGGCGTGCGGTGGTCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATCCCGGTGCGGGCCAGCCATGCCTCGAATTcc                                                                                                                                        >  5:287526/1‑149 (MQ=255)
                    gCGTGCGGTGGTCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCt                                                                                                                              >  6:87090/1‑149 (MQ=255)
                      gTGCGGTGGTCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTc                                                                                                                            <  5:66326/149‑1 (MQ=255)
                              gtCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCa                                                                                                                     <  2:106382/148‑1 (MQ=255)
                                    ggTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCgg                                                                                                                                 <  6:211039/130‑1 (MQ=255)
                                    ggTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCgg                                                                                                                                 >  5:211039/1‑130 (MQ=255)
                                      tCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGgc                                                                                                            >  3:109222/1‑149 (MQ=255)
                                         tCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGgcgcg                                                                                                         <  5:85389/149‑1 (MQ=255)
                                            cTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTTGGTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCtt                                                                                                                             >  1:270006/1‑127 (MQ=255)
                                            cTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCtt                                                                                                                             <  2:270006/126‑1 (MQ=255)
                                            cTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGt                                                                                                      <  1:338423/149‑1 (MQ=255)
                                                gCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCACGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACaa                                                                                                  >  4:206530/1‑149 (MQ=255)
                                                gCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCACGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACaa                                                                                                  >  4:201667/1‑149 (MQ=255)
                                                     cTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGt                                                                                             <  1:231900/149‑1 (MQ=255)
                                                     cTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCACGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGt                                                                                             <  3:201667/149‑1 (MQ=255)
                                                     cTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCACGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGt                                                                                             <  3:206530/149‑1 (MQ=255)
                                                        tgCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCAt                                                                                          <  7:30510/149‑1 (MQ=255)
                                                          ccAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATg                                                                                         >  2:116816/1‑148 (MQ=255)
                                                                           gCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTa                                                                       <  4:109222/149‑1 (MQ=255)
                                                                            ccGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAg                                                                      >  2:257352/1‑149 (MQ=255)
                                                                             cggcggCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGt                                                                     <  5:50289/149‑1 (MQ=255)
                                                                              ggcggcAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGtt                                                                    <  5:87090/149‑1 (MQ=255)
                                                                                 ggcAGGTT‑GTCACCGAAGCGCACCCGGCCGTCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAg                                                                                  >  5:249244/1‑132 (MQ=255)
                                                                                 ggcAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAg                                                                                  <  6:249244/132‑1 (MQ=255)
                                                                                  gcAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAg                                                                 >  3:36664/1‑148 (MQ=255)
                                                                                      gTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTc                                                                                                                                               >  2:66784/1‑66 (MQ=255)
                                                                                      gTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTc                                                                                                                                               <  1:66784/66‑1 (MQ=255)
                                                                                       tt‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGa                                                                  >  6:80534/1‑142 (MQ=255)
                                                                                       tt‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGa                                                                  <  5:80534/142‑1 (MQ=255)
                                                                                       tt‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCatgat                                                           <  1:116816/149‑1 (MQ=255)
                                                                                        t‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCCCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCatgatg                                                          >  4:154372/1‑149 (MQ=255)
                                                                                              ccGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATgtcgtc                                                     >  1:353535/1‑149 (MQ=255)
                                                                                                    cgcACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGcc                                               >  3:298012/1‑149 (MQ=255)
                                                                                                           ggCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATg                                        >  8:199393/1‑149 (MQ=255)
                                                                                                              cATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGcc                                      >  2:90978/1‑148 (MQ=255)
                                                                                                                  ggCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCatga                                 <  1:257352/149‑1 (MQ=255)
                                                                                                                    ccTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGGCGGAGCCGGGGATGCCCatgatg                               <  4:36664/149‑1 (MQ=255)
                                                                                                                         gTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATGAAGtt                          >  7:290078/1‑149 (MQ=255)
                                                                                                                           gAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATGAAGTTGa                        >  4:303609/1‑149 (MQ=255)
                                                                                                                              gATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCatga                                                            <  7:197935/110‑1 (MQ=255)
                                                                                                                              gATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCatga                                                            >  8:197935/1‑110 (MQ=255)
                                                                                                                                 aCCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATGAAGTTGATGCCgg                  >  3:354840/1‑149 (MQ=255)
                                                                                                                                 aCCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATGAAGTTGATGCCgg                  >  2:84699/1‑149 (MQ=255)
                                                                                                                                         gggCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATGAAGTTGATGCCGGCCAcccc           >  5:184474/1‑148 (MQ=255)
                                                                                                                                               gccaTGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATTAAGTTGATGCCGGCCACCCCCagcagc    <  4:298012/149‑1 (MQ=255)
                                                                                                                                                  aTGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATGAAGTTGATGCCGGCCACCCCCAGCAGCGt  >  2:117605/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAAGGGTTGTCGGGGGTAGGCGTGCGGTGGTCCATGGTCATCCCTATGCAAGCTGTGCCAAGGCCTGGCGGAAGGCCGGCGGCAGGTT‑GTCACCGAAGCGCACCCGGCCATCGGCCTGGGTGAAGATACCGGTGCGGGCCAGCCATGCCTCGAATTCCGGCCCGGGCTTCAGGCCCAGGGTCTGGCGCGCGTACAACGCGTCATGGAACGAGGTGGTCTGGTAGTTGAGCATGATGTCGTCGGAGCCGGGGATGCCCATGATGAAGTTGATGCCGGCCACCCCCAGCAGCGT  >  NC_002947/628759‑629051

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: ‑