Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,912 G→T W30C (TGG→TGT folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9120GT96.4% 89.3 / NA 28W30C (TGG→TGTfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (12/15);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACG  >  NZ_CP009273/49828‑49992
                                                                                    |                                                                                
cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaa                                                                             >  2:339174/1‑90 (MQ=255)
cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTaa                                                                             >  2:273264/1‑90 (MQ=255)
 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTaaa                                                                            >  2:597666/1‑90 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                       <  1:263292/90‑1 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                       <  2:252812/90‑1 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                       >  2:609358/1‑90 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                       <  1:615660/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac                                                                      <  1:278934/90‑1 (MQ=255)
           cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                  <  2:580375/90‑1 (MQ=255)
            ggcgTTAGCGTTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                 >  2:48952/1‑90 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                 >  2:70926/1‑90 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                 >  2:449422/1‑90 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACAGCtt                                                                 >  2:87782/1‑90 (MQ=255)
                 tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTtaaata                                                            <  1:449422/90‑1 (MQ=255)
                      gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAAccc                                                       <  2:211302/90‑1 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac                                                                      <  1:565461/68‑1 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac                                                                      >  2:565461/1‑68 (MQ=255)
                                  aTCGGCATGAAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                           <  1:48952/90‑1 (MQ=255)
                                   tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc                                          <  1:273264/90‑1 (MQ=255)
                                        aTGAAAACCGCCATGCCGGGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                     <  1:123074/90‑1 (MQ=255)
                                            aaaaCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                 <  2:114258/90‑1 (MQ=255)
                                             aaccGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                                <  2:649658/87‑1 (MQ=255)
                                                    aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                     >  2:238496/1‑78 (MQ=255)
                                                    aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                     <  1:238496/78‑1 (MQ=255)
                                                         gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg                    >  1:620514/1‑90 (MQ=255)
                                                           ggAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                  <  1:339174/90‑1 (MQ=255)
                                                                 tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt            >  2:36013/1‑90 (MQ=255)
                                                                           tctcGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACg  >  1:322768/1‑90 (MQ=255)
                                                                                    |                                                                                
CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACG  >  NZ_CP009273/49828‑49992

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