Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,229,956 A→G K405E (AAA→GAA)  preA → NAD‑dependent dihydropyrimidine dehydrogenase subunit PreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,229,9560AG90.9% 55.4 / ‑0.1 22K405E (AAA→GAA) preANAD‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/9);  total (12/10)
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.91e-01

GAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGG  >  NZ_CP009273/2229869‑2230041
                                                                                       |                                                                                     
gagaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaa                                                                                     >  2:509351/1‑90 (MQ=255)
 agaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaag                                                                                    <  2:214278/90‑1 (MQ=255)
  gaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaaga                                                                                   >  1:129916/1‑90 (MQ=255)
      tgtgtgGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAacac                                                                               <  2:253119/90‑1 (MQ=255)
       gtgtgGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACAcc                                                                              >  1:60032/1‑90 (MQ=255)
        tgtgGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACAccc                                                                             <  2:27868/90‑1 (MQ=255)
                 tCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGtt                                                                    >  2:652034/1‑90 (MQ=255)
                   tGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGt                                                                  >  1:568989/1‑90 (MQ=255)
                      tttGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTaaa                                                               <  2:231965/90‑1 (MQ=255)
                      tttGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTaaa                                                               <  2:162998/90‑1 (MQ=255)
                                          gTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCAt                                           >  1:10428/1‑90 (MQ=255)
                                               ttGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGt                                      <  2:10428/90‑1 (MQ=255)
                                                      gAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTc                               >  2:200846/1‑90 (MQ=255)
                                                      gAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTc                               >  2:567124/1‑90 (MQ=255)
                                                        gCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGc                             >  1:551147/1‑90 (MQ=255)
                                                        gCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGc                             >  2:349178/1‑90 (MQ=255)
                                                           cGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCgg                          >  1:194482/1‑90 (MQ=255)
                                                           cGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCgg                          >  1:678845/1‑90 (MQ=255)
                                                                aagtgaagtTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCg                     >  2:605659/1‑90 (MQ=255)
                                                                      agtTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGa               <  2:223980/90‑1 (MQ=255)
                                                                      agtTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGa               <  1:712804/90‑1 (MQ=255)
                                                                      agtTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGa               <  1:262162/90‑1 (MQ=255)
                                                                                 ggCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATcc    >  2:558528/1‑90 (MQ=255)
                                                                                   cGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCgg  <  1:459029/90‑1 (MQ=255)
                                                                                   cGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCgg  <  2:678845/90‑1 (MQ=255)
                                                                                       |                                                                                     
GAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGG  >  NZ_CP009273/2229869‑2230041

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