Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,905	T→G	L28R (CTC→CGC)	folA →	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,905	0	T	G	96.3%	86.8 / ‑3.2	27	L28R (CTC→CGC)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (20/6);  total (21/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTG  >  NC_000913/49843‑49962                                                               |                                                          gTTAGCGGTCTATCGTTTTATCGGCATGGAAACCGCCATGTTGTGGACCCTCCCTGCCGATcgcg                                                         <  1:3209717/65‑1 (MQ=255)       ggTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATCGCGCCTGGt                                                   <  1:1425485/65‑1 (MQ=255)          aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTa                                                >  1:3500661/1‑65 (MQ=255)            aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                              >  1:3154710/1‑65 (MQ=255)            aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                              >  1:1671916/1‑65 (MQ=255)            aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                              >  1:250723/1‑65 (MQ=255)            aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGGTTaaa                                              >  1:1322577/1‑65 (MQ=255)            aTCGCGTTATCGGCATGGAAAACGCCATGCTGTAGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                                              >  1:2600114/1‑65 (MQ=255)               gcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGc                                           >  1:2656662/1‑65 (MQ=255)                   tATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAaca                                       <  1:306497/65‑1 (MQ=255)                    aTCGTCATGGAAAACGCGGTGCTGTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGGAacac                                      <  1:696106/65‑1 (MQ=255)                    aTCGGCATGGAAACCGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAacac                                      <  1:1313131/65‑1 (MQ=255)                         cATGGAAAACGCCATGCTGGGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaa                                 <  1:1767867/65‑1 (MQ=255)                         cATGAAAACCGCCGTGCTGTGGATCCAGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaa                                 <  1:794985/65‑1 (MQ=255)                                        gctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                  >  1:2243239/4‑65 (MQ=255)                                         ctgTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACCCAACACCTTAAATAAACCCGTGATTAt                 >  1:1876265/3‑65 (MQ=255)                                               aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACTCGTGATTATGGgccg           >  1:1646752/1‑65 (MQ=255)                                               aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg           >  1:269807/1‑65 (MQ=255)                                               aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg           >  1:3132515/1‑65 (MQ=255)                                               aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg           >  1:1471002/1‑65 (MQ=255)                                               aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg           >  1:370193/1‑65 (MQ=255)                                               aaCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg           >  1:3757374/1‑65 (MQ=255)                                                aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc          >  1:1465306/1‑65 (MQ=255)                                                aCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc          >  1:1403024/1‑65 (MQ=255)                                                 cctgcctgCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc         >  1:2422897/1‑65 (MQ=255)                                                 cctgcctgCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc         >  1:1327088/1‑65 (MQ=255)                                                  ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa        >  1:2511036/1‑65 (MQ=255)                                                  ctgcctgcCGATCGCGCCTGGTTTAAACGCAACACATTAAATAAACCCGTGATTATGGGCCGCCa        >  1:2881117/1‑65 (MQ=255)                                                        gccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGCTTATGGGCCCCCATACCTg  >  1:2579346/1‑65 (MQ=255)                                                        gccGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg  >  1:2310583/1‑65 (MQ=255)                                                               |                                                          GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTG  >  NC_000913/49843‑49962

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑