Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,831	T→C	S3S (AGT→AGC)	folA →	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,831	0	T	C	100.0%	227.5 / NA	62	S3S (AGT→AGC)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/33);  total (29/33)

AGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCC  >  NC_000913/49696‑49955                                                                                                                                        |                                                                                                                             aGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTg                                                                                                                           >  1:654280/1‑139 (MQ=255)         gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcggc                                                                                                                   <  1:1060373/139‑1 (MQ=255)         gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcggc                                                                                                                   <  2:1336231/139‑1 (MQ=255)          agagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcggcg                                                                                                                  <  1:604424/139‑1 (MQ=255)           gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                 <  2:654280/139‑1 (MQ=255)             gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTa                                                                                                               >  2:33479/1‑139 (MQ=255)              agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAg                                                                                                              >  1:1192094/1‑139 (MQ=255)                 gTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCgg                                                                                                           >  1:407067/1‑139 (MQ=255)                   ggACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTa                                                                                                         >  1:1129746/1‑139 (MQ=255)                   ggACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTa                                                                                                         <  2:1307008/139‑1 (MQ=255)                      cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGAt                                                                                                      <  1:692854/139‑1 (MQ=255)                       tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATc                                                                                                     <  2:1192094/139‑1 (MQ=255)                         gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGc                                                                                                             <  1:272908/129‑1 (MQ=255)                         gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGc                                                                                                             >  2:272908/1‑129 (MQ=255)                            agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGtt                                                                                                <  1:386752/139‑1 (MQ=255)                            agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGtt                                                                                                <  2:466152/139‑1 (MQ=255)                              cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                              <  2:1444210/139‑1 (MQ=255)                              cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                              >  1:351072/1‑139 (MQ=255)                              cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                              >  1:458298/1‑139 (MQ=255)                                       aaaTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                                                                                     <  1:33479/139‑1 (MQ=255)                                                         tcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                   >  1:432363/1‑139 (MQ=255)                                                         tcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                   <  2:432363/139‑1 (MQ=255)                                                           cTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcc                                                                 <  2:1129746/139‑1 (MQ=255)                                                              gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTa                                                                                               <  1:975278/106‑1 (MQ=255)                                                              gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTa                                                                                               >  2:975278/1‑106 (MQ=255)                                                              gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgc                                                              >  2:1320303/1‑139 (MQ=255)                                                                 ccAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcctg                                                           <  1:1320303/139‑1 (MQ=255)                                                                  cAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcg                                                                                                                     <  2:1171558/80‑1 (MQ=255)                                                                  cAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcg                                                                                                                     >  1:1171558/1‑80 (MQ=255)                                                                    gttgacgacGGTTTACGCTTTACGTATAGTGGCGGCAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                 <  1:1071720/79‑1 (MQ=255)                                                                    gttgacgacGGTTTACGCTTTACGTATAGTGGCGGCAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                 >  2:1071720/4‑82 (MQ=255)                                                                    gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCg                                                        >  2:629825/4‑139 (MQ=255)                                                                    gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCg                                                        >  1:619533/4‑139 (MQ=255)                                                                           acgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGc                                                 <  2:351072/139‑1 (MQ=255)                                                                           acgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGc                                                 <  1:250105/139‑1 (MQ=255)                                                                            cgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGcc                                                >  2:935923/1‑139 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                 >  2:753955/1‑71 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                 <  1:753955/71‑1 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                                                                                     >  1:1450527/1‑99 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                                                                                     <  2:1450527/99‑1 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTgg                                             >  2:546040/1‑139 (MQ=255)                                                                                         aCGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                   <  1:629825/139‑1 (MQ=255)                                                                                          cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGAACGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGTCTGCCGATCTCGCCTGGTTTAAACGCaa                                  <  1:935923/139‑1 (MQ=255)                                                                                                gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGcc                                                                              >  2:774481/1‑89 (MQ=255)                                                                                                gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGcc                                                                              <  1:774481/89‑1 (MQ=255)                                                                                                   gCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaa                         >  1:734980/1‑139 (MQ=255)                                                                                                    cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTTTCGGCATGGAAAACGCCATGCCGTGGaa                                                                   >  1:790542/1‑96 (MQ=255)                                                                                                    cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTTTCGGCATGGAAAACGCCATGCCGTGGaa                                                                   <  2:790542/96‑1 (MQ=255)                                                                                                    cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                        >  1:1460328/1‑139 (MQ=255)                                                                                                     gACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                   <  2:1161565/95‑1 (MQ=255)                                                                                                     gACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                   >  1:1161565/1‑95 (MQ=255)                                                                                                           tttttttATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                   <  1:162259/89‑1 (MQ=255)                                                                                                           tttttttATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa                                                                   >  2:162259/1‑89 (MQ=255)                                                                                                               tttATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                   <  2:1288834/117‑1 (MQ=255)                                                                                                               tttATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                   >  1:1288834/1‑117 (MQ=255)                                                                                                                        aaaTCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg    >  1:248886/1‑139 (MQ=255)                                                                                                                        aaaTCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg    <  1:546040/139‑1 (MQ=255)                                                                                                                          aTCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc  <  2:407067/139‑1 (MQ=255)                                                                                                                                aTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                   >  2:1078041/1‑100 (MQ=255)                                                                                                                                aTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                   <  1:1078041/100‑1 (MQ=255)                                                                                                                                   aTCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                                                                             >  1:1417897/1‑55 (MQ=255)                                                                                                                                   aTCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                                                                             <  2:1417897/55‑1 (MQ=255)                                                                                                                                        |                                                                                                                             AGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCC  >  NC_000913/49696‑49955

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