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,7780TC97.9% 150.2 / ‑3.6 48E215G (GAA→GGA) folCbifunctional tetrahydrofolate synthase/dihydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (31/16);  total (32/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

GTGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCTTTCTCG  >  NZ_CP009273/2425691‑2425864
                                                                                       |                                                                                      
gTGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTccc                                                                                      <  1:1290357/90‑1 (MQ=255)
gTGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTccc                                                                                      >  2:1412082/1‑90 (MQ=255)
gTGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTccc                                                                                      <  1:469388/90‑1 (MQ=255)
 tGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCt                                                                                     <  2:409761/90‑1 (MQ=255)
   ccGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTgc                                                                                   <  1:1128921/90‑1 (MQ=255)
    cGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTgcg                                                                                  <  2:301096/90‑1 (MQ=255)
     gTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTgcgc                                                                                 <  1:1412082/90‑1 (MQ=255)
        aGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCcac                                                                              >  2:1370244/1‑90 (MQ=255)
         gCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCcaca                                                                             <  2:1452568/90‑1 (MQ=255)
            tCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCa                                                                          >  2:452731/1‑90 (MQ=255)
             cGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAg                                                                         >  1:74286/1‑90 (MQ=255)
             cGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAg                                                                         >  1:679912/1‑90 (MQ=255)
              gCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGc                                                                        <  2:114290/90‑1 (MQ=255)
                tAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCaa                                                                      <  1:232564/90‑1 (MQ=255)
                tAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCaa                                                                      <  2:1123332/90‑1 (MQ=255)
                tAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCaa                                                                      <  2:1030143/90‑1 (MQ=255)
                 aaaCGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAAt                                                                     <  1:571297/90‑1 (MQ=255)
                      cccAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAg                                                                >  2:664092/1‑90 (MQ=255)
                      cccAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAg                                                                >  2:603451/1‑90 (MQ=255)
                         aGTTATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTCTTCCCTGCGCCACATCAGCAATGGTAGAAg                                                             >  2:6647/1‑90 (MQ=255)
                         aGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAg                                                             >  1:112940/1‑90 (MQ=255)
                               gATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCAttt                                                       >  2:570054/1‑90 (MQ=255)
                                    gTGACGGAATAGTTCCACTCGACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGc                                                  >  1:620379/1‑90 (MQ=255)
                                    gTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGc                                                  >  2:1001668/1‑90 (MQ=255)
                                    gTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGc                                                  >  1:933135/1‑90 (MQ=255)
                                    gTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGc                                                  >  1:883564/1‑90 (MQ=255)
                                    gTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGc                                                  >  2:258990/1‑90 (MQ=255)
                                         ggAATAGTTCCACTCAACGCCCCGACGTTGTAGCAGTGCCCCTTTTCCCTGCGCCACATCAGCAATGGTAGCAGGCATTTCCGGCTCAcc                                             >  1:1033231/1‑90 (MQ=255)
                                                   cACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGcc                                   <  2:160050/90‑1 (MQ=255)
                                                    aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCg                                  >  2:660399/1‑90 (MQ=255)
                                                    aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCa                                               <  2:115041/77‑1 (MQ=255)
                                                    aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCa                                               >  1:115041/1‑77 (MQ=255)
                                                         aCGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGtttt                             >  1:89419/1‑90 (MQ=255)
                                                               cGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTg                       >  2:1399033/1‑90 (MQ=255)
                                                                           aGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCt           >  1:417456/1‑90 (MQ=255)
                                                                            gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAgg                                                            <  2:62305/40‑1 (MQ=38)
                                                                            gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAgg                                                            >  1:62305/1‑40 (MQ=38)
                                                                            gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGcc                                   <  1:19036/65‑1 (MQ=255)
                                                                            gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGcc                                   >  2:19036/1‑65 (MQ=255)
                                                                            gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg          >  2:614237/1‑90 (MQ=255)
                                                                            gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg          >  2:798923/1‑90 (MQ=255)
                                                                            gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg          >  2:885552/1‑90 (MQ=255)
                                                                             tGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGc         >  1:577507/1‑90 (MQ=255)
                                                                              gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCt                        <  1:824378/74‑1 (MQ=255)
                                                                              gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCt                        >  2:824378/1‑74 (MQ=255)
                                                                              gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCt        >  2:162418/1‑90 (MQ=255)
                                                                              gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCt        >  2:141725/1‑90 (MQ=255)
                                                                              gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCt        >  1:57912/1‑90 (MQ=255)
                                                                                aCCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCttt      >  2:1080271/1‑90 (MQ=255)
                                                                                 ccTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGc                                                  >  1:921811/1‑45 (MQ=255)
                                                                                 ccTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGc                                                  <  2:921811/45‑1 (MQ=255)
                                                                                    tttcCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCTTTCTcg  <  2:59509/90‑1 (MQ=255)
                                                                                       |                                                                                      
GTGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCTTTCTCG  >  NZ_CP009273/2425691‑2425864

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