Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,903 T→A D27E (GAT→GAA folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9030TA95.5% 63.9 / ‑4.2 22D27E (GAT→GAAfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (8/13);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

AATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTC  >  NZ_CP009273/49822‑49976
                                                                                 |                                                                         
aaTGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTg                                                                   <  2:13537/90‑1 (MQ=255)
      cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTaa                                                             >  1:34488/1‑90 (MQ=255)
            gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                       <  2:34488/90‑1 (MQ=255)
            gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                       <  1:516926/90‑1 (MQ=255)
            gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCaa                                                       <  1:563254/90‑1 (MQ=255)
             aTTGCGGCGTTAGCGGTAGCTCGCGTTATCGGCATGGAAAACGCCACGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAac                                                      <  2:532342/90‑1 (MQ=255)
             aTTGCGGCGTTAGCGGGAGATCGCGTTATCGGCATGGGAAACCCCAATCCGCGGAAGCCCCCCGCCGAACTCCCCTTGGTTAAACggcaa                                                      >  1:382916/1‑86 (MQ=255)
                      ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaat                                             >  2:146638/1‑90 (MQ=255)
                      ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGGTTAAACGCAACACCTtaaat                                             >  1:377135/1‑90 (MQ=255)
                       tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaata                                            <  1:416148/90‑1 (MQ=255)
                        aGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTtaaataa                                           <  1:176479/90‑1 (MQ=255)
                               gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                    >  2:566787/1‑90 (MQ=255)
                               gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTg                                    >  1:510300/1‑90 (MQ=255)
                                   gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCCGCCTGCCGAACTCGCCTGGGTTAAACGCAACACCTTAAATAAACCCGTGATTa                                >  1:533818/1‑90 (MQ=255)
                                       taccGGCATGGAAACCGCCATGCCGTGGACCCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                            <  2:274902/87‑1 (MQ=255)
                                         tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc                          <  2:235893/90‑1 (MQ=255)
                                             cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa                      <  2:469770/90‑1 (MQ=255)
                                                  aaaaCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                 <  1:184666/90‑1 (MQ=255)
                                                          aTGCCGTGGAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca         >  1:261144/1‑90 (MQ=255)
                                                                 ggAACCTGCCTGCCGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc  <  1:172553/90‑1 (MQ=255)
                                                                       tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc  >  2:128484/1‑84 (MQ=255)
                                                                       tgcctgccGAACTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc  <  1:128484/84‑1 (MQ=255)
                                                                                 |                                                                         
AATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTC  >  NZ_CP009273/49822‑49976

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: