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%	203.2 / NA	56	S3S (AGT→AGC)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/28);  total (28/28)

AGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAA  >  NC_000913/49696‑49936                                                                                                                                        |                                                                                                          aGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTg                                                                                                        >  2:1052367/1‑139 (MQ=255)         gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcggc                                                                                                <  2:69766/139‑1 (MQ=255)           gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  1:113874/139‑1 (MQ=255)           gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  1:369347/139‑1 (MQ=255)             gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTa                                                                                            >  1:529577/1‑139 (MQ=255)               gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTg                                                                                                    >  2:115232/1‑129 (MQ=255)               gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTg                                                                                                    <  1:115232/129‑1 (MQ=255)                      cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGAt                                                                                   <  2:974964/139‑1 (MQ=255)                                   tataAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  1:333532/115‑1 (MQ=255)                                   tataAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              >  2:333532/1‑115 (MQ=255)                                           tttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGTCAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGGAGATCGCGTTATCGGCa                                                                      <  2:772365/131‑1 (MQ=255)                                           tttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCa                                                                      >  1:772365/1‑131 (MQ=255)                                            ttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACg                                                             <  2:529577/139‑1 (MQ=255)                                              ccTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGcc                                                           <  2:1335598/139‑1 (MQ=255)                                                         tcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGt                                                                                       >  1:358783/1‑100 (MQ=255)                                                         tcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGt                                                                                       <  2:358783/100‑1 (MQ=255)                                                              gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgc                                           >  2:128969/1‑139 (MQ=255)                                                                    gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCg                                     >  2:166958/4‑139 (MQ=255)                                                                       gacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAtc                                  >  2:760496/1‑139 (MQ=255)                                                                          gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGAt                                                                                   <  1:1116588/87‑1 (MQ=255)                                                                          gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGAt                                                                                   >  2:1116588/1‑87 (MQ=255)                                                                           acgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGc                              <  1:128969/139‑1 (MQ=255)                                                                            cgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCAt                                                         <  2:40127/111‑1 (MQ=255)                                                                            cgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCAt                                                         >  1:40127/1‑111 (MQ=255)                                                                              gTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              >  2:1141132/1‑72 (MQ=255)                                                                              gTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  1:1141132/72‑1 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGAtt                                                                                                     <  2:738231/64‑1 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGAtt                                                                                                     >  1:738231/1‑64 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATcg                                                                                 <  2:688901/84‑1 (MQ=255)                                                                               tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATcg                                                                                 >  1:688901/1‑84 (MQ=255)                                                                                    gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcg                                                                                                  <  2:729484/62‑1 (MQ=255)                                                                                    gCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcg                                                                                                  >  1:729484/1‑62 (MQ=255)                                                                                         aCGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGttt                       <  1:325612/132‑1 (MQ=255)                                                                                         aCGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                >  2:299883/1‑139 (MQ=255)                                                                                         aCGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATAAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGttt                       >  2:325612/1‑132 (MQ=255)                                                                                              taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              >  2:1032555/1‑56 (MQ=255)                                                                                              taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  1:1032555/56‑1 (MQ=255)                                                                                                gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  2:1021894/54‑1 (MQ=255)                                                                                                gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              >  1:1021894/1‑54 (MQ=255)                                                                                                gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              >  2:799671/1‑54 (MQ=255)                                                                                                gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  1:799671/54‑1 (MQ=255)                                                                                                 tGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTta        >  1:452557/1‑139 (MQ=255)                                                                                                  ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  1:1138014/52‑1 (MQ=255)                                                                                                  ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              >  2:1138014/1‑52 (MQ=255)                                                                                                  ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaa       >  2:1099518/1‑139 (MQ=255)                                                                                                    cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGttc                                                                            >  2:784981/1‑67 (MQ=255)                                                                                                    cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTa                                                                            <  1:784981/68‑1 (MQ=255)                                                                                                    cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAt                                   >  2:738860/1‑109 (MQ=255)                                                                                                    cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAt                                   <  1:738860/109‑1 (MQ=255)                                                                                                     gACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat     <  1:592289/138‑1 (MQ=255)                                                                                                     gACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat     >  2:592289/1‑138 (MQ=255)                                                                                                      aCAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              <  2:1142846/48‑1 (MQ=255)                                                                                                      aCAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                              >  1:1142846/1‑48 (MQ=255)                                                                                                       cAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaataaa  <  1:1099518/139‑1 (MQ=255)                                                                                                                               aaTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACAcc           <  1:639088/106‑1 (MQ=255)                                                                                                                               aaTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACAcc           >  2:639088/1‑106 (MQ=255)                                                                                                                                        |                                                                                                          AGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAA  >  NC_000913/49696‑49936

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