Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,910 T→A W30R (TGG→AGG)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9100TA97.2% 121.5 / NA 37W30R (TGG→AGG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base A (15/20);  minor base C (1/0);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49828‑49996
                                                                                  |                                                                                      
cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                                                 >  1:218705/1‑90 (MQ=255)
cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTaa                                                                                 >  2:795128/1‑90 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                           <  1:827998/90‑1 (MQ=255)
      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                           <  2:47646/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                          <  1:238508/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                          <  1:318354/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                          >  1:92375/1‑90 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                          <  1:667745/90‑1 (MQ=255)
       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGCCATGGAAAACACCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac                                                                          <  1:572397/90‑1 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCtt                                                                     >  2:704592/1‑90 (MQ=255)
            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCtt                                                                     >  2:704058/1‑90 (MQ=255)
                 tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTtaaata                                                                <  2:807634/90‑1 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCCTGCCGGGGAACCTGCCTGCCCAACTCGCCAGGGTTAAAAGCAAAACCTTAAATAAACCCGGGATTa                                                    >  2:702727/1‑90 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGGACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  2:872875/1‑90 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                           >  1:123371/1‑67 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                           <  2:123371/67‑1 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGGTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  1:614736/1‑90 (MQ=255)
                             gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCGGCCTGCCGAACTCGCCCGGGTTAAACGCAAAACCTTAAATAAACCCCTGATTa                                                    >  1:862355/1‑90 (MQ=255)
                              cgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                                                   >  1:179450/1‑90 (MQ=255)
                                ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                           <  1:5491/64‑1 (MQ=255)
                                ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa                                                                           >  2:5491/1‑64 (MQ=255)
                                  aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc                                               <  2:92375/90‑1 (MQ=255)
                                       cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                                          <  2:460768/90‑1 (MQ=255)
                                        aTGGAAACCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  1:587290/90‑1 (MQ=255)
                                        aTGGAAACCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  1:702727/90‑1 (MQ=255)
                                        aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         >  1:479641/1‑90 (MQ=255)
                                         tGGAAAACGCCATGCCGTGGACCCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                                        <  1:639210/90‑1 (MQ=255)
                                                   cATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc                              >  2:553807/1‑90 (MQ=255)
                                                         gTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg                        <  1:872875/90‑1 (MQ=255)
                                                           ggAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  1:245875/90‑1 (MQ=255)
                                                                 tgcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  2:346723/84‑1 (MQ=255)
                                                                 tgcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      >  1:346723/1‑84 (MQ=255)
                                                                  gcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGtt               >  1:654516/1‑90 (MQ=255)
                                                                      gccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCa           >  2:183819/1‑90 (MQ=255)
                                                                              cGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa   <  1:493963/90‑1 (MQ=255)
                                                                               gCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  2:196527/90‑1 (MQ=255)
                                                                               gCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  2:393592/90‑1 (MQ=255)
                                                                                  |                                                                                      
CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49828‑49996

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