Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,884	C→T	P21L (CCG→CTG)	folA →	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,884	0	C	T	100.0%	56.6 / NA	17	P21L (CCG→CTG)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/4);  total (13/4)

AATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTG  >  NC_000913/49822‑49942                                                               |                                                           aaTGATCAGTCTGAGTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGt                                                          <  1:2420084/65‑1 (MQ=255)       cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAcc                                                    <  1:802342/65‑1 (MQ=255)       cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAcc                                                    <  1:2743098/65‑1 (MQ=255)         gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAAcctg                                                  <  1:3092444/65‑1 (MQ=255)                    gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAAAGCCAAGCTGTGGAAACTTCCTGCCGAtcg                                       >  1:3074042/1‑64 (MQ=255)                     cgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGCtcgc                                      >  1:2660363/1‑63 (MQ=255)                      tttAGCGGTAGACCGCGTTATCGGCATGGAAACCCCCATGCTGTGGAACCTGCCTGCCGATcgcg                                     <  1:1004693/64‑1 (MQ=255)                      gTTAGCGGTAGATCGCGTTATCGGCATGGAAACCGCCATGCTGTGGACCCTGCCTGCCGATcgcg                                     <  1:3866678/65‑1 (MQ=255)                      gTTAGCGGTAGATCGCGTTAGCGGCATGGAAAACGCCTTGCTGTGGAACCTGCCTGCCGATcgcg                                     <  1:796317/65‑1 (MQ=255)                            gctAGACCGGGTTAGGGGAATGGAAAACGCGTTGCGGTGGAGCCTGCCTGCCGATCGCGCCTGGt                               <  1:2693585/63‑1 (MQ=255)                             gTAGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGtt                              >  1:817160/1‑65 (MQ=255)                               aGATCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTa                            >  1:1223018/1‑65 (MQ=255)                                 aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                          >  1:285277/1‑65 (MQ=255)                                 aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                          >  1:2865174/1‑65 (MQ=255)                                 aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                          >  1:2104965/1‑65 (MQ=255)                                 aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                          >  1:3109264/1‑65 (MQ=255)                                 aTCGCGTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTaaa                          >  1:487093/1‑65 (MQ=255)                                   cgcgTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACg                        >  1:3237357/1‑65 (MQ=255)                                      gTTATCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACgtct                     >  1:2375853/1‑62 (MQ=255)                                         aTCGGCATGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAacac                  <  1:461149/65‑1 (MQ=255)                                         aTCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAacac                  <  1:2237580/65‑1 (MQ=255)                                          tCGGCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACAcc                 <  1:3231433/65‑1 (MQ=255)                                            ggCATGGAAAACGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCtt               >  1:3378340/1‑65 (MQ=255)                                              cATGGAAACCGCCATGCTTTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaa             <  1:13646/65‑1 (MQ=255)                                              cATGGAAAACGCCATGCTGTGGACCCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTtaa             <  1:2793909/65‑1 (MQ=255)                                                   aaaaCGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAAACCCTtcaataac        >  1:3683377/1‑64 (MQ=255)                                                    aaccGCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAAc       <  1:1322103/62‑1 (MQ=255)                                                        gCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt   >  1:1613565/1‑65 (MQ=255)                                                        gCCATGCTGTGGAACCTGCCTGCCGATCGCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGt   >  1:3185976/1‑65 (MQ=255)                                                         ccATGCTGTGGAACCTGCCTGCCGATCGCCCCTGGTTTAAACGCAACACCTTACATAAACTCGTg  >  1:1795755/1‑65 (MQ=255)                                                               |                                                           AATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTG  >  NC_000913/49822‑49942

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