Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,550,898 T→C I381V (ATC→GTC)  eutB ← ethanolamine ammonia‑lyase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,550,8980TC97.1% 116.5 / ‑4.7 35I381V (ATC→GTC) eutBethanolamine ammonia‑lyase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (17/17);  total (17/18)
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.93e-01

ATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTT  >  NZ_CP009273/2550814‑2550982
                                                                                    |                                                                                    
aTGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGAccatca                                                                                 >  2:786046/1‑90 (MQ=255)
 tGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAg                                                                                <  2:202442/90‑1 (MQ=255)
 tGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAg                                                                                <  1:584180/90‑1 (MQ=255)
       ggtCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGATGGCGAGCAGGACCATCAGGTTTTc                                                                          <  1:459062/90‑1 (MQ=255)
        gtCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCg                                                                         >  1:788651/1‑90 (MQ=255)
            ggTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGa                                                                     <  2:802052/90‑1 (MQ=255)
                gTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGtt                                                                 <  2:532833/90‑1 (MQ=255)
                 ttGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTc                                                                >  1:290459/1‑90 (MQ=255)
                  tGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCt                                                               <  2:318194/90‑1 (MQ=255)
                      cATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGgtc                                                           <  2:583423/90‑1 (MQ=255)
                            gTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGgtc                                                           >  1:251358/1‑84 (MQ=255)
                            gTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGgtc                                                           <  2:251358/84‑1 (MQ=255)
                                tcaCCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGtcagcgtcagcg                                                 <  2:788651/90‑1 (MQ=255)
                                    ccAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                             >  2:581235/1‑90 (MQ=255)
                                     cAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggtt                                            >  2:849479/1‑90 (MQ=255)
                                     cAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggtt                                            >  1:916189/1‑90 (MQ=255)
                                     cAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggtt                                            <  2:916189/90‑1 (MQ=255)
                                      aGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGa                                                                     <  2:937458/64‑1 (MQ=255)
                                      aGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGa                                                                     >  1:937458/1‑64 (MQ=255)
                                      aGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGgtc                                                           >  1:671774/1‑74 (MQ=255)
                                      aGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGgtc                                                           <  2:671774/74‑1 (MQ=255)
                                      aGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggttg                                           <  2:290459/90‑1 (MQ=255)
                                      aGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTAAGCGtggttg                                           <  2:432486/90‑1 (MQ=255)
                                       gCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggttgg                                          >  1:892502/1‑90 (MQ=255)
                                       gCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggtt                                            >  1:1048919/1‑88 (MQ=255)
                                       gCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggtt                                            <  2:1048919/88‑1 (MQ=255)
                                                ccATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAAcagc                                 <  1:1059846/90‑1 (MQ=255)
                                                   tgatgTAGTTGCAGCCTGCGGTGGCTAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGt                              >  1:448733/1‑90 (MQ=255)
                                                      tgTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                           <  2:763436/90‑1 (MQ=255)
                                                      tgTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                           >  2:256733/1‑90 (MQ=255)
                                                      tgTAGTTGCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                           >  2:799911/1‑90 (MQ=255)
                                                        tAGTTGCATCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAg                         >  1:251503/1‑90 (MQ=255)
                                                             gCAGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCAt                    >  1:895231/1‑90 (MQ=255)
                                                               aGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATag                  >  1:954030/1‑90 (MQ=255)
                                                               aGCCTGCGGTGGCGAGCAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATag                  >  1:819490/1‑90 (MQ=255)
                                                                               cAGGACCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCtt  <  2:251503/90‑1 (MQ=255)
                                                                                    |                                                                                    
ATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTT  >  NZ_CP009273/2550814‑2550982

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