Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 50,280 T→C F153S (TTT→TCT)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927350,2800TC89.3% 68.6 / 2.4 28F153S (TTT→TCT) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (5/20);  total (6/22)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCTA  >  NZ_CP009273/50197‑50366
                                                                                   |                                                                                      
cccGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGAt                                                                                  <  1:140533/90‑1 (MQ=255)
           aGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcggc                                                                       <  2:59928/90‑1 (MQ=255)
              cGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGt                                                                    <  2:544030/90‑1 (MQ=255)
              cGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGt                                                                    <  1:63846/90‑1 (MQ=255)
                gATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTaa                                                                  <  2:304225/90‑1 (MQ=255)
                   gACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt                                                               <  2:29220/90‑1 (MQ=255)
                             cGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt                                                               <  1:53465/80‑1 (MQ=255)
                             cGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt                                                               <  1:560754/80‑1 (MQ=255)
                             cGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt                                                               >  2:53465/1‑80 (MQ=255)
                             cGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGAAAttt                                                               >  2:560754/1‑80 (MQ=255)
                                tATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTa                                                  <  1:355360/90‑1 (MQ=255)
                                            tCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCg                                      <  1:836273/90‑1 (MQ=255)
                                                     cTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTTCGGCTAACCCCCGATGCGACg                             >  1:521429/1‑90 (MQ=255)
                                                               gAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGt                   <  2:394907/90‑1 (MQ=255)
                                                                 aCTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTTCGGCTAGCGCCGGGTGCGGCGCCGGTCGCGTCt                 >  2:55649/1‑90 (MQ=255)
                                                                          gCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:418934/90‑1 (MQ=255)
                                                                          gCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:188844/90‑1 (MQ=255)
                                                                          gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:100285/90‑1 (MQ=255)
                                                                          gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:500025/90‑1 (MQ=255)
                                                                          gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:243547/90‑1 (MQ=255)
                                                                          gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:740628/90‑1 (MQ=255)
                                                                           cTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCt       <  2:367833/90‑1 (MQ=255)
                                                                            tATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGGTGGGACGCCGGGCGCGTCTTTTCCGGCCtt      >  1:653008/1‑90 (MQ=255)
                                                                              ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATATAATTTACGGCTAGCGCCGGAGGCGACGCCGGGCGCGTCTTATCCGGCCCTcc    >  1:549904/1‑90 (MQ=38)
                                                                              ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc    <  1:740430/90‑1 (MQ=255)
                                                                              ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc    <  1:412143/90‑1 (MQ=255)
                                                                              ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc    <  2:776429/90‑1 (MQ=255)
                                                                                gCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCta  <  2:785102/90‑1 (MQ=255)
                                                                                   |                                                                                      
CCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCTA  >  NZ_CP009273/50197‑50366

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