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

GAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCA  >  NC_000913/49706‑49956                                                                                                                              |                                                                                                                              gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  <  1:1198995/139‑1 (MQ=255) gagagCGTGCACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  <  2:978694/139‑1 (MQ=255)  agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  >  2:576077/1‑138 (MQ=255)  agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  <  1:576077/138‑1 (MQ=255)   gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTa                                                                                                                >  2:1220564/1‑139 (MQ=255)                    cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                               >  1:1024176/1‑139 (MQ=255)                                   tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  >  2:963726/1‑105 (MQ=255)                                   tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  <  1:963726/105‑1 (MQ=255)                                   tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  <  1:1296815/105‑1 (MQ=255)                                   tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  >  2:1296815/1‑105 (MQ=255)                                   tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  <  2:472606/105‑1 (MQ=255)                                   tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  >  1:472606/1‑105 (MQ=255)                                   tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGc                                                                                >  2:454243/1‑139 (MQ=255)                                                   cGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg                                                                <  2:1402621/139‑1 (MQ=255)                                                          gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcgg                                                                                                                     <  1:230905/76‑1 (MQ=255)                                                          gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcgg                                                                                                                     >  2:230905/4‑79 (MQ=255)                                                                             ttACGTATAGTGGCGACACTTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  <  2:1012579/63‑1 (MQ=255)                                                                             ttACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  >  1:1012579/1‑63 (MQ=255)                                                                                      gTGGCGAGAATTTTCTTCATTGGGAAATCTCATTGTCCAGCCGCATGGTGGCTTTACCGGTCGCTCGCGTTATCGTCATGGACAACGCCATACCGTCGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACCCCtt                             >  2:684469/1‑139 (MQ=255)                                                                                      gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                               >  2:446162/1‑73 (MQ=255)                                                                                      gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                                                                               <  1:446162/73‑1 (MQ=255)                                                                                      gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcct                                                                 <  2:44044/103‑1 (MQ=255)                                                                                      gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcct                                                                 >  1:44044/1‑103 (MQ=255)                                                                                      gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCtt                             >  1:52512/1‑139 (MQ=255)                                                                                        ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  <  1:733900/52‑1 (MQ=255)                                                                                        ggCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                                                                                  >  2:733900/1‑52 (MQ=255)                                                                                         gCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGc                                                                                                              >  1:120944/1‑55 (MQ=255)                                                                                         gCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGc                                                                                                              <  2:120944/55‑1 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTa                                                                                                >  1:1429440/1‑68 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTa                                                                                                <  2:1429440/68‑1 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                                                                                      >  2:82365/1‑78 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa                                                                                      <  1:82365/78‑1 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGcc                                                                               <  1:915122/85‑1 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                         >  2:45249/1‑139 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                         >  1:451189/1‑139 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                         >  1:1155484/1‑139 (MQ=255)                                                                                          cGACAATTTTTTTTATCGGGAAATCTCAATGACCCGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGcc                                                                               >  2:915122/1‑85 (MQ=255)                                                                                           gACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                   <  2:356648/128‑1 (MQ=255)                                                                                           gACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCaa                                   >  1:356648/1‑128 (MQ=255)                                                                                           gACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                        <  1:1220564/139‑1 (MQ=255)                                                                                           gACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaata                        >  1:463223/1‑139 (MQ=255)                                                                                                 tttttttATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg                  <  2:953509/139‑1 (MQ=255)                                                                                                          cGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg                  <  2:580656/130‑1 (MQ=255)                                                                                                          cGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCg                  >  1:580656/1‑130 (MQ=255)                                                                                                            ggAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc       <  1:45249/139‑1 (MQ=255)                                                                                                               aaTCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc    <  2:1155484/139‑1 (MQ=255)                                                                                                                 tctcAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa  <  2:451189/139‑1 (MQ=255)                                                                                                                      aTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGgaa                                                                                    <  2:828840/52‑1 (MQ=255)                                                                                                                      aTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaa                                                                                    >  1:828840/1‑52 (MQ=255)                                                                                                                      aTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                         >  2:526571/1‑111 (MQ=255)                                                                                                                      aTGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTtaaat                         <  1:526571/111‑1 (MQ=255)                                                                                                                          tCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                                                                              >  1:738941/1‑54 (MQ=255)                                                                                                                          tCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                                                                              <  2:738941/54‑1 (MQ=255)                                                                                                                              |                                                                                                                              GAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCA  >  NC_000913/49706‑49956

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