Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,883 C→T P21S (CCG→TCG)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,8830CT97.7% 132.7 / NA 44P21S (CCG→TCG) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/22);  minor base G (1/0);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAT  >  NZ_CP009273/49796‑49957
                                                                                       |                                                                          
gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCg                                                                          >  1:279632/1‑90 (MQ=255)
gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCg                                                                          >  1:86642/1‑90 (MQ=255)
         ttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAAcctg                                                                 <  1:80845/90‑1 (MQ=255)
              tCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACctgcctgc                                                            >  2:293044/1‑90 (MQ=255)
              tCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACctgcctgc                                                            <  2:76046/90‑1 (MQ=255)
                        tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGcc                                                  <  2:377551/90‑1 (MQ=255)
                        tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGcc                                                  <  1:293044/90‑1 (MQ=255)
                        tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGcc                                                  <  2:439622/90‑1 (MQ=255)
                        tcAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGcc                                                  <  1:216474/90‑1 (MQ=255)
                               tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                           >  2:316042/1‑90 (MQ=255)
                               tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                           <  1:391069/90‑1 (MQ=255)
                                cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                          >  1:59268/1‑90 (MQ=255)
                                cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                          >  2:420730/1‑90 (MQ=255)
                                cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGGGGAACCTGCCTGCCGATCTCGCCTGGTTTaa                                          >  2:380134/1‑90 (MQ=255)
                                 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGcc                                                  <  2:391368/81‑1 (MQ=255)
                                 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGcc                                                  >  1:391368/1‑81 (MQ=255)
                                 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTg                                                <  1:435178/83‑1 (MQ=255)
                                 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTg                                                >  2:435178/1‑83 (MQ=255)
                                 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                           >  2:249072/1‑88 (MQ=255)
                                 aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                           <  1:249072/88‑1 (MQ=255)
                                  gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCAGGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAAc                                        <  2:130070/90‑1 (MQ=255)
                                      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                    <  1:580563/90‑1 (MQ=255)
                                      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                    <  1:564449/90‑1 (MQ=255)
                                      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                    <  1:421179/90‑1 (MQ=255)
                                      gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                    <  1:328971/90‑1 (MQ=255)
                                       aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAac                                   >  1:283323/1‑90 (MQ=255)
                                           cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCt                               <  1:250152/90‑1 (MQ=255)
                                            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                              >  1:446077/1‑90 (MQ=255)
                                            ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                              >  2:378469/1‑90 (MQ=255)
                                              cgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaa                            >  1:25363/1‑90 (MQ=255)
                                                 tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                         >  2:11882/1‑90 (MQ=255)
                                                           tCGCGTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGAt               >  2:534022/1‑90 (MQ=255)
                                                             gcgTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                    >  2:85657/1‑67 (MQ=255)
                                                             gcgTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                    <  1:85657/67‑1 (MQ=255)
                                                             gcgTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa             >  2:229766/1‑90 (MQ=255)
                                                             gcgTTATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa             >  1:505235/1‑90 (MQ=255)
                                                                ttATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg          <  1:69618/90‑1 (MQ=255)
                                                                ttATCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGGTTAAACGCAACACCTTAAATAAACCCGTGATTATgg          >  2:386995/1‑90 (MQ=255)
                                                                  aTCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAACCCCGTGATTATGGgc        <  2:166244/90‑1 (MQ=255)
                                                                  aTCGGCATGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc        <  2:579779/90‑1 (MQ=255)
                                                                        aTGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                           >  2:39253/1‑49 (MQ=255)
                                                                        aTGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTa                                           <  1:39253/49‑1 (MQ=255)
                                                                        aTGGAAAACGCCATGTCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt  <  1:332854/90‑1 (MQ=255)
                                                                        aTGGAAAACGCCATGGCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGGTTATGGgccgccct  >  1:52664/1‑88 (MQ=255)
                                                                                       |                                                                          
GACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAT  >  NZ_CP009273/49796‑49957

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