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,9120GT93.9% 100.0 / 1.6 33W30C (TGG→TGTfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (13/18);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

AGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49829‑49996
                                                                                   |                                                                                    
aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa                                                                                >  1:156840/1‑90 (MQ=255)
 gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAAc                                                                               >  2:566524/1‑90 (MQ=255)
 gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAAc                                                                               >  2:379798/1‑90 (MQ=255)
     gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAACCGCCATGCCTTGGACCCTCCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  2:477529/90‑1 (MQ=255)
      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac                                                                          <  2:486706/90‑1 (MQ=255)
      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac                                                                          <  2:248643/90‑1 (MQ=255)
      aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac                                                                          <  1:732008/90‑1 (MQ=255)
          cggcgTGAGCGGTAGATCGCTTTATCGGCATGGAAACCCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      <  2:146296/90‑1 (MQ=255)
           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  2:620321/1‑90 (MQ=255)
           ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  1:124725/1‑90 (MQ=255)
               ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCTCCTGTTTTAAAGACAAAACCTtaaat                                                                 >  1:544117/1‑90 (MQ=255)
               ttAGCGGTAGATCGCGTTATCGGCATGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTtaaat                                                                 <  2:568918/90‑1 (MQ=255)
                       aGATCGGGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt                                                         <  1:613999/90‑1 (MQ=255)
                       aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt                                                         <  1:292839/90‑1 (MQ=255)
                       aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt                                                         <  2:307986/90‑1 (MQ=255)
                       aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt                                                         <  1:395535/90‑1 (MQ=255)
                            gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  2:466369/1‑90 (MQ=255)
                                tATCGGCATGGAAAACGCCATGCCGGGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                <  1:453288/90‑1 (MQ=255)
                                   cGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAAccc                                                           <  1:484552/76‑1 (MQ=255)
                                   cGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAAccc                                                           >  2:484552/1‑76 (MQ=255)
                                     gCATGAAAACCCCCATGCCGGGGACCCTCCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                                           <  2:424616/90‑1 (MQ=255)
                                       aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         >  2:75075/1‑90 (MQ=255)
                                        tGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTg                                                        >  2:3213/1‑74 (MQ=255)
                                        tGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTg                                                        <  1:3213/74‑1 (MQ=255)
                                              aCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACAcc                                                                       >  2:511237/1‑53 (MQ=255)
                                              aCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACAcc                                                                       <  1:511237/53‑1 (MQ=255)
                                                   aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                             >  2:267711/1‑90 (MQ=255)
                                                        gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg                        >  1:652895/1‑90 (MQ=255)
                                                          ggACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAACCCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  1:264181/90‑1 (MQ=255)
                                                                     gccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                 <  2:587065/52‑1 (MQ=255)
                                                                     gccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                 >  1:587065/1‑52 (MQ=255)
                                                                         aTCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc       <  2:143338/90‑1 (MQ=255)
                                                                              gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  1:75075/90‑1 (MQ=255)
                                                                                   |                                                                                    
AGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49829‑49996

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