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

GGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGG  >  NC_000913/49702‑49964                                                                                                                                  |                                                                                                                                      gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcg                                                                                                                              <  1:1249896/139‑1 (MQ=255)     gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          <  2:855348/139‑1 (MQ=255)        agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAg                                                                                                                       >  2:1581554/1‑139 (MQ=255)                  cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATcg                                                                                                             >  2:125006/1‑139 (MQ=255)                        cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                                       <  1:1581554/139‑1 (MQ=255)                                     tttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGc                                                                                                                      <  2:1605604/111‑1 (MQ=255)                                     tttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGc                                                                                                                      >  1:1605604/1‑111 (MQ=255)                                       tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCAt                                                                                                 >  2:203694/1‑130 (MQ=255)                                       tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAGCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCAt                                                                                                 <  1:203694/130‑1 (MQ=255)                                          tCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCAt                                                                                     <  1:467446/139‑1 (MQ=255)                                              catcatCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCg                                                                                 >  1:1554858/1‑139 (MQ=255)                                                        gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          >  1:1238453/1‑88 (MQ=255)                                                        gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          <  2:1238453/88‑1 (MQ=255)                                                         cACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                                       >  1:1164023/1‑106 (MQ=255)                                                         cACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                                       <  2:1164023/106‑1 (MQ=255)                                                         cACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcc                                                                      >  2:922615/1‑139 (MQ=255)                                                              gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGAt                                                                                                                                  <  2:814911/71‑1 (MQ=255)                                                              gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGAt                                                                                                                                  >  1:814911/4‑74 (MQ=255)                                                              gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCg                                                                 >  2:535910/4‑139 (MQ=255)                                                                     acgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGc                                                          <  1:140377/139‑1 (MQ=255)                                                                         tttACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTgg                                                      >  2:1226178/1‑139 (MQ=255)                                                                                 ttACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          >  1:625144/1‑63 (MQ=255)                                                                                 ttACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          <  2:625144/63‑1 (MQ=255)                                                                                   aCGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                            <  1:125006/139‑1 (MQ=255)                                                                                       ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcg                                                                                                                              >  1:53862/1‑53 (MQ=255)                                                                                       ataGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcg                                                                                                                              <  2:53862/53‑1 (MQ=255)                                                                                        taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaa                                                                                             <  1:283371/85‑1 (MQ=255)                                                                                        taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaa                                                                                             >  2:283371/1‑85 (MQ=255)                                                                                        taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACAcc                                       >  2:625746/1‑139 (MQ=255)                                                                                         aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          <  1:857594/55‑1 (MQ=255)                                                                                         aGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          >  2:857594/1‑55 (MQ=255)                                                                                          gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                                     >  1:530754/1‑139 (MQ=255)                                                                                            ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGTCTGCCGATCTCGCCTGGTTTAAACGCa                                            >  2:1014454/1‑130 (MQ=255)                                                                                            ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAGCGCAACACCTtaa                                   >  1:379659/1‑139 (MQ=255)                                                                                            ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCa                                            <  1:1014454/130‑1 (MQ=255)                                                                                            ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaa                                   >  2:689834/1‑139 (MQ=255)                                                                                             gCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg                                                                        <  2:253806/101‑1 (MQ=255)                                                                                             gCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg                                                                        >  1:253806/1‑101 (MQ=255)                                                                                              cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          <  1:20130/50‑1 (MQ=255)                                                                                              cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          >  2:20130/1‑50 (MQ=255)                                                                                              cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          >  1:1188810/1‑50 (MQ=255)                                                                                              cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGc                                                                                                                      >  1:1570187/1‑54 (MQ=255)                                                                                              cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGc                                                                                                                      <  2:1570187/54‑1 (MQ=255)                                                                                              cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                                 >  2:1150476/1‑139 (MQ=255)                                                                                              cGACAATTTATTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                          <  2:1188810/50‑1 (MQ=255)                                                                                                     tttttttATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAGCACCTTAAATAAACCCg                          <  1:922615/139‑1 (MQ=255)                                                                                                     tttttttATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg                          <  1:1226178/139‑1 (MQ=255)                                                                                                               gggAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAGCGCAACACCTTAAATAAACCCGTGATTATggg                <  2:379659/139‑1 (MQ=255)                                                                                                                   aaTCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc            <  1:689834/139‑1 (MQ=255)                                                                                                                         aaTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAcc      <  1:625746/139‑1 (MQ=255)                                                                                                                             aTCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg  <  1:368826/139‑1 (MQ=255)                                                                                                                               cAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCa                                                                                                  <  1:863709/41‑1 (MQ=255)                                                                                                                               cAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCa                                                                                                  >  2:863709/1‑41 (MQ=255)                                                                                                                                  |                                                                                                                                      GGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGG  >  NC_000913/49702‑49964

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