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% 183.6 / ‑0.8 55W30C (TGG→TGTfolAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (27/26);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49827‑49996
                                                                                     |                                                                                    
tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTa                                                                                  <  2:977996/90‑1 (MQ=255)
tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTa                                                                                  >  1:911515/1‑90 (MQ=255)
 cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                                                                 >  2:701688/1‑90 (MQ=255)
 cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCCGCCGATCTCGCCTGTTTTaa                                                                                 >  2:1173843/1‑90 (MQ=255)
  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa                                                                                >  1:1311938/1‑90 (MQ=255)
  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa                                                                                >  1:877904/1‑90 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  2:911515/90‑1 (MQ=255)
       gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa                                                                           <  1:299835/90‑1 (MQ=255)
          tGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTGCCGGTCTCGCCTGTTTTAAACGCAacac                                                                        >  1:844023/1‑90 (MQ=255)
            cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      <  1:1508269/90‑1 (MQ=255)
            cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      <  1:446392/90‑1 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  1:1092251/1‑90 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  2:753906/1‑90 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  2:1132412/1‑90 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt                                                                     >  1:1320507/1‑90 (MQ=255)
             ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGGACCTTCCTGCCGATCTCGTCTGGTTTAAACGCAACACCtt                                                                     >  2:1423020/1‑90 (MQ=255)
                         aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt                                                         <  1:254681/90‑1 (MQ=255)
                          gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTg                                                        >  1:37640/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  1:1084272/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  1:123694/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa                                                    >  1:1483363/1‑90 (MQ=255)
                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCCGCCTGGCGATCTCGCCCGTTTTAAACGCAAACCCTTAAATAAACCCGTGATTa                                                    >  1:1555265/1‑90 (MQ=255)
                                 ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg                                                 <  2:576227/90‑1 (MQ=255)
                                 ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAAACCGTGATTATgg                                                 >  2:666312/1‑90 (MQ=255)
                                  tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                <  2:533413/90‑1 (MQ=255)
                                  tATCGGCATGAAAACCCCCATCCGGGGGAACCCCCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                <  1:986972/90‑1 (MQ=255)
                                    tgggCATGAAAACCCCCCTGCGGTGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc                                              <  1:1173843/88‑1 (MQ=255)
                                     cGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                                             >  1:1239813/1‑90 (MQ=255)
                                         aggaaaaaCGCCTTGCCGGGAAACCGGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  1:1344316/86‑1 (MQ=255)
                                         aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      >  2:588766/1‑61 (MQ=255)
                                         aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  2:675854/90‑1 (MQ=255)
                                         aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         >  1:968360/1‑90 (MQ=255)
                                         aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt                                         <  2:550933/90‑1 (MQ=255)
                                         aTGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt                                                                      <  1:588766/61‑1 (MQ=255)
                                          tGGAAACCGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa                                        <  1:587977/90‑1 (MQ=255)
                                             aaaaCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt                                     <  2:1320507/90‑1 (MQ=255)
                                                     aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca                             >  2:659921/1‑90 (MQ=255)
                                                      tGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa                            >  2:23080/1‑90 (MQ=255)
                                                          gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg                        >  2:434503/1‑90 (MQ=255)
                                                            ggACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACCCCTTAAATAACCCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  1:197816/90‑1 (MQ=255)
                                                            ggAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  2:713312/90‑1 (MQ=255)
                                                            ggAACCTCCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc                      <  1:378179/90‑1 (MQ=255)
                                                             gaccctccctgCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcg                     <  2:1081447/87‑1 (MQ=255)
                                                                 ctgcctgcCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg                 <  2:1239813/90‑1 (MQ=255)
                                                                 ctgcctgcCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg                 >  1:442184/1‑90 (MQ=255)
                                                                  tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                                           <  2:1173665/63‑1 (MQ=255)
                                                                  tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                                           >  1:1173665/1‑63 (MQ=255)
                                                                  tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt                >  2:1037403/1‑90 (MQ=255)
                                                                  tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt                >  1:1515171/1‑90 (MQ=255)
                                                                  tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt                >  2:907079/1‑90 (MQ=255)
                                                                       gccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCa           >  2:627518/1‑90 (MQ=255)
                                                                                gCCTGTTTTAAACGCAACCCCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  2:904921/90‑1 (MQ=255)
                                                                                gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  1:161925/90‑1 (MQ=255)
                                                                                gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  2:1518945/90‑1 (MQ=255)
                                                                                gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTAGGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa  <  1:884225/90‑1 (MQ=255)
                                                                                     |                                                                                    
TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA  >  NZ_CP009273/49827‑49996

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