Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,910	T→G	W30G (TGG→GGG)	folA →	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,910	0	T	G	100.0%	70.2 / NA	21	W30G (TGG→GGG)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (9/12);  total (9/12)

CGGTAGATCGCGTTATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGG  >  NC_000913/49848‑49964                                                                |                                                       cgttAGATCGCGTTATCGGCATGGAAACCGCCATG‑CGGGGAACCCTGCCTGCCGATCTCGCCggg                                                      <  1:2403200/62‑1 (MQ=255)   ttAGATCGCGTTACCGG‑ATGGGAAAGCCCGTGCCCGTGCACCCTCCCTGCCGATCTCGCCGGGtt                                                    <  1:2191399/64‑1 (MQ=255)       aTCGCGTTATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGCTCTCGCCGGGTTTaaa                                                >  1:1645108/1‑65 (MQ=255)       aTCGCGTTATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTaaa                                                >  1:503712/1‑65 (MQ=255)       aTCGCGTTATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCGCCGGGGTTaaa                                                >  1:1288555/1‑65 (MQ=255)       aTCGCGTTATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCCCCGGGTTTaaa                                                >  1:1466804/1‑65 (MQ=255)       aTCGCGTTATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCCCCGGGTTTaaa                                                >  1:2287155/1‑65 (MQ=255)          gcgTTATCGGCATGGAAAACGCCATG‑CCGTGTTACCTGCGTGCCGATAGATACTGGTTTAAACGc                                             >  1:3272703/1‑65 (MQ=255)          gcgTCATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGc                                             >  1:2172194/1‑65 (MQ=255)            gTTATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa                                           >  1:2936014/1‑65 (MQ=255)               aTCGGCATAAAAACCCCCCTG‑CCGGGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAacac                                        <  1:489539/65‑1 (MQ=255)                tCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                       <  1:3228228/65‑1 (MQ=255)                tCGGCATGCAACCCGCTTCG‑CTGTGTAACCGCCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                       <  1:2649822/65‑1 (MQ=255)                tCGGCATGAAAACCGCCATG‑CCGTGGACCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                       <  1:1323479/65‑1 (MQ=255)                gcGGCATGGAAAGCGCCGTG‑CCGTGGACCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                       <  1:978545/64‑1 (MQ=255)                gcGGCATGAAACCCCCCAGG‑CGGGGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                       <  1:2689853/64‑1 (MQ=255)                acGGCTTGAAAAACCCCATG‑CCGTGGACCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc                                       <  1:4578/64‑1 (MQ=255)                    gaTGGAAACCGCCATG‑CCGTGAAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaa                                   <  1:2525766/64‑1 (MQ=255)                    cATGAAAAACCCCAGG‑CCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaa                                   <  1:2740362/65‑1 (MQ=255)                           accgcgATG‑CCGTGGACCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAcc                            <  1:2985404/63‑1 (MQ=255)                           accGCCATG‑CGGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAcc                            <  1:2532177/63‑1 (MQ=255)                             cGCCATG‑CCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCg                          >  1:2906024/1‑65 (MQ=255)                               ccATG‑CCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAccctgt                        >  1:246736/1‑62 (MQ=255)                                     ccGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt                   >  1:532025/1‑65 (MQ=255)                                           aaCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg             >  1:2546056/1‑65 (MQ=255)                                           aaCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg             >  1:1504409/1‑65 (MQ=255)                                            aCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc            >  1:2568960/1‑65 (MQ=255)                                            aCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc            >  1:348494/1‑65 (MQ=255)                                                gcctgccGATCTCGCCGGGCTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa        >  1:2618873/1‑65 (MQ=255)                                                      cGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg  >  1:2665158/1‑65 (MQ=255)                                                                |                                                       CGGTAGATCGCGTTATCGGCATGGAAAACGCCATG‑CCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGG  >  NC_000913/49848‑49964

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