Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,425,778 T→C E215G (GAA→GGA)  folC ← bifunctional tetrahydrofolate synthase/dihydrofolate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,425,7780TC93.5% 85.8 / ‑0.5 31E215G (GAA→GGA) folCbifunctional tetrahydrofolate synthase/dihydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (20/9);  total (21/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.96e-01

GCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCT  >  NZ_CP009273/2425693‑2425858
                                                                                     |                                                                                
gCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTg                                                                              >  2:617604/1‑90 (MQ=255)
    tGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGcc                                                                          <  2:611310/90‑1 (MQ=255)
        cATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACAt                                                                      <  1:36818/90‑1 (MQ=255)
          tCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCa                                                                    <  1:112339/90‑1 (MQ=255)
          tCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCa                                                                    >  1:329940/1‑90 (MQ=255)
                     ccAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGa                                                         <  1:556507/90‑1 (MQ=255)
                     ccAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGa                                                         <  2:620753/90‑1 (MQ=255)
                       aGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAg                                                       <  1:689264/90‑1 (MQ=255)
                         tCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGc                                                     >  1:491151/1‑90 (MQ=255)
                            tGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCAtt                                                  >  2:492272/1‑90 (MQ=255)
                                  gTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTgc                                                                             >  2:836250/1‑57 (MQ=255)
                                  gTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTgc                                                                             <  1:836250/57‑1 (MQ=255)
                                               tCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTg                               >  2:797052/1‑90 (MQ=255)
                                               tCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTg                               >  2:122911/1‑90 (MQ=255)
                                                 cACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGcc                             <  1:162392/90‑1 (MQ=255)
                                                  aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCg                            >  2:239903/1‑90 (MQ=255)
                                                  aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCg                            >  2:718192/1‑90 (MQ=255)
                                                  aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCg                            >  1:249503/1‑90 (MQ=255)
                                                       aCGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGtttt                       >  2:443270/1‑90 (MQ=255)
                                                            ccGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCt                  >  1:644611/1‑90 (MQ=255)
                                                               aCGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCg               >  2:197332/1‑90 (MQ=255)
                                                               aCGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCg               >  2:132191/1‑90 (MQ=255)
                                                                          gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg    >  2:175710/1‑90 (MQ=255)
                                                                          gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg    >  2:65662/1‑90 (MQ=255)
                                                                          gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg    >  1:35104/1‑90 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGc                                                     >  2:595408/1‑39 (MQ=38)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGc                                                     <  1:595408/39‑1 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGcc      >  2:559010/1‑86 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGcc      <  1:559010/86‑1 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCt  >  2:55770/1‑90 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCt  >  2:562131/1‑90 (MQ=255)
                                                                                     |                                                                                
GCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCT  >  NZ_CP009273/2425693‑2425858

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