Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,765	C→T	intergenic (+134/‑58)	kefC → / → folA	potassium:proton antiporter/dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,765	0	C	T	100.0%	181.9 / NA	56	intergenic (+134/‑58)	kefC/folA	potassium:proton antiporter/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/27);  total (29/27)

AGAGTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCT  >  NC_000913/49635‑49897                                                                                                                                   |                                                                                                                                     agagTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACgg                                                                                                                              >  1:772791/1‑139 (MQ=255)   agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                            >  2:916167/1‑139 (MQ=255)   agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                            >  2:757478/1‑139 (MQ=255)   agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                            >  1:263733/1‑139 (MQ=255)   agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                            >  1:978672/1‑139 (MQ=255)   agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                            >  1:55645/1‑139 (MQ=255)   agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                            >  1:433057/1‑139 (MQ=255)       aCGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACg                                                                                                                        >  2:575719/1‑139 (MQ=255)       aCGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACg                                                                                                                        >  2:180858/1‑139 (MQ=255)       aCGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACg                                                                                                                        >  2:148138/1‑139 (MQ=255)       aCGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCAGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACg                                                                                                                        >  1:844772/1‑139 (MQ=255)        cGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGTAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGc                                                                                                                       <  2:754476/139‑1 (MQ=255)         gTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                                                      <  1:803789/139‑1 (MQ=255)          tAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCtt                                                                                                                     <  1:304279/139‑1 (MQ=255)                       aCAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGc                                                                                                        <  2:433057/139‑1 (MQ=255)                            tAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACaa                                                                                                   >  1:179255/1‑139 (MQ=255)                              aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                                                      <  2:72407/118‑1 (MQ=255)                              aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                                                      >  1:72407/1‑118 (MQ=255)                              aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAAtt                                                                                                 >  2:216960/1‑139 (MQ=255)                                tgtttgttTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAAtttt                                                                                               >  1:355778/1‑139 (MQ=255)                                     gtttttgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTa                                                                                          >  1:883601/1‑139 (MQ=255)                                           gtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                    <  2:772791/139‑1 (MQ=255)                                           gtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                    <  1:383300/139‑1 (MQ=255)                                                 ttGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc                                                                              <  1:180858/139‑1 (MQ=255)                                                 ttGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc                                                                              <  1:529192/139‑1 (MQ=255)                                                     aaTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATg                                                                          <  2:390183/139‑1 (MQ=255)                                                      aTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa                                                                         <  1:575719/139‑1 (MQ=255)                                                       tGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa                                                                                  <  2:322403/129‑1 (MQ=255)                                                       tGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa                                                                                  >  1:322403/1‑129 (MQ=255)                                                        gcggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                                                       <  1:804305/139‑1 (MQ=255)                                                         cggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAt                                                                        >  2:725875/1‑137 (MQ=255)                                                         cggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAt                                                                        <  1:725875/137‑1 (MQ=255)                                                            cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATg                                                                          >  1:272460/1‑132 (MQ=255)                                                            cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATg                                                                          <  2:272460/132‑1 (MQ=255)                                                                   aGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTAc                                                                                                                         >  1:541522/1‑78 (MQ=255)                                                                   aGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTAc                                                                                                                         <  2:541522/78‑1 (MQ=255)                                                                   aGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgc                                                            <  1:1019047/139‑1 (MQ=255)                                                                      gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcggc                                                         <  2:883601/139‑1 (MQ=255)                                                                        gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                       <  2:844772/139‑1 (MQ=255)                                                                        gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                       <  1:355715/139‑1 (MQ=255)                                                                        gaaagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                       <  2:100633/136‑1 (MQ=255)                                                                         agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                                                      >  1:250262/1‑139 (MQ=255)                                                                             cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                       <  2:33485/134‑1 (MQ=255)                                                                             cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                       >  1:33485/1‑134 (MQ=255)                                                                               tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                    >  2:358489/1‑103 (MQ=255)                                                                               tGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                    <  1:358489/103‑1 (MQ=255)                                                                                     cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAt                                                                        >  2:442364/1‑109 (MQ=255)                                                                                     cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAt                                                                        <  1:442364/109‑1 (MQ=255)                                                                                      gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa                                                                                   <  2:470756/97‑1 (MQ=255)                                                                                      gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa                                                                                   >  1:470756/1‑97 (MQ=255)                                                                                           cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                    <  2:263733/139‑1 (MQ=255)                                                                                                                aCATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGcc               <  1:216960/139‑1 (MQ=255)                                                                                                                    catcCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTgg           >  2:851185/1‑139 (MQ=255)                                                                                                                      tcCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa         >  2:815374/1‑139 (MQ=255)                                                                                                                             aCCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgcct  >  1:1028613/1‑139 (MQ=255)                                                                                                                             aCCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAACAGGCTGATTGCGGCGTTAGCCGTAGATTGCGTTATCGGCATGGAAAACCCCATGCCGCGGAAccttcct  >  1:685836/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                     AGAGTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCT  >  NC_000913/49635‑49897

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: ‑