Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,767	A→G	intergenic (+136/‑56)	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,767	0	A	G	100.0%	246.8 / NA	67	intergenic (+136/‑56)	kefC/folA	potassium:proton antiporter/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/33);  total (34/33)

AGTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTG  >  NC_000913/49637‑49894                                                                                                                                   |                                                                                                                                agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                         >  1:785458/1‑139 (MQ=255) agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                         >  2:1615998/1‑139 (MQ=255) agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                         >  2:1704902/1‑139 (MQ=255)  gTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGttt                                                                                                                        <  1:1704902/139‑1 (MQ=255)   tGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGGACCAGTTGGCGACGGTTTa                                                                                                                       <  1:1154783/139‑1 (MQ=255)   tGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                                       >  2:1154783/1‑139 (MQ=255)    gACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTAc                                                                                                                      <  2:785458/139‑1 (MQ=255)      cGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGc                                                                                                                    <  1:1615998/139‑1 (MQ=255)        tAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCtt                                                                                                                  <  2:199743/139‑1 (MQ=255)        tAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCtt                                                                                                                  <  2:834415/139‑1 (MQ=255)         aaaTCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCttt                                                                                                                 <  2:1079089/139‑1 (MQ=255)              acacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGta                                                                                                            <  2:1721182/139‑1 (MQ=255)               cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGtat                                                                                                           >  2:723709/1‑139 (MQ=255)                 cTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAg                                                                                                         >  2:94151/1‑139 (MQ=255)                         cTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACa                                                                                                 >  1:537835/1‑139 (MQ=255)                           aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGc                                                                                                                    >  1:86016/1‑118 (MQ=255)                           aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGc                                                                                                                    <  2:86016/118‑1 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCttt                                                                                                                 >  2:52676/1‑120 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCttt                                                                                                                 <  1:52676/120‑1 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCtt                                                                                                                  >  2:967489/1‑119 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCtt                                                                                                                  <  1:967489/119‑1 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAAtt                                                                                              >  2:820271/1‑139 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAAtt                                                                                              >  1:113079/1‑139 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATCTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAAtt                                                                                              >  2:474579/1‑139 (MQ=255)                              tgtttgttTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAAtttt                                                                                            >  1:1167707/1‑139 (MQ=255)                                                       cggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCa                                                                   >  2:464509/1‑139 (MQ=255)                                                         gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACg                                                                                                                     <  1:1321024/87‑1 (MQ=255)                                                         gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACg                                                                                                                     >  2:1321024/1‑87 (MQ=255)                                                         gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGt                                                                 >  2:1584021/1‑139 (MQ=255)                                                         gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGt                                                                 >  1:1448313/1‑139 (MQ=255)                                                          cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                         >  1:75954/1‑82 (MQ=255)                                                          cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGtt                                                                                                                         <  2:75954/82‑1 (MQ=255)                                                          cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTc                                                                >  2:1195226/1‑139 (MQ=255)                                                                  gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg                                                        >  2:261829/1‑139 (MQ=255)                                                                  gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATCTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg                                                        <  1:474579/139‑1 (MQ=255)                                                                    gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcggc                                                      <  2:537835/139‑1 (MQ=255)                                                                    gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcggc                                                      <  1:94151/139‑1 (MQ=255)                                                                    gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcggc                                                      <  2:1448313/139‑1 (MQ=255)                                                                      gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                    <  1:1195226/139‑1 (MQ=255)                                                                      gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCAACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                    <  2:1167707/139‑1 (MQ=255)                                                                       agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt                                                   >  2:1580959/1‑139 (MQ=255)                                                                        gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                                                   >  1:1551710/1‑74 (MQ=255)                                                                        gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                                                   <  2:1551710/74‑1 (MQ=255)                                                                           cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                    >  2:989110/1‑134 (MQ=255)                                                                           cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                                                    <  1:989110/134‑1 (MQ=255)                                                                            gTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGc                                                                                                                    >  1:435140/1‑69 (MQ=255)                                                                            gTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGc                                                                                                                    <  2:435140/69‑1 (MQ=255)                                                                               gACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                                       <  2:12716/63‑1 (MQ=255)                                                                               gACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                                       >  1:12716/1‑63 (MQ=255)                                                                                    gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                                       >  2:1213599/1‑58 (MQ=255)                                                                                    gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                                       <  1:1213599/58‑1 (MQ=255)                                                                                    gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATcgc                                      >  1:1077186/1‑139 (MQ=255)                                                                                        gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTa                                  <  1:1580959/139‑1 (MQ=255)                                                                                         cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                 <  1:820271/139‑1 (MQ=255)                                                                                                      tttCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaa                      >  2:1687429/1‑137 (MQ=255)                                                                                                      tttCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaaa                      <  1:1687429/137‑1 (MQ=255)                                                                                                      tttCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAAc                    >  2:1138231/1‑139 (MQ=255)                                                                                                        tCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                 >  1:1005942/1‑76 (MQ=255)                                                                                                        tCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg                                                        >  2:1409350/1‑101 (MQ=255)                                                                                                        tCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg                                                        <  1:1409350/101‑1 (MQ=255)                                                                                                        tCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTGTTTTAtagcga                                                                                 <  2:1005942/76‑6 (MQ=255)                                                                                                               catcatCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCg           <  2:113079/139‑1 (MQ=255)                                                                                                                  catcCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTgg        >  2:864155/1‑139 (MQ=255)                                                                                                                        cGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg  <  1:723709/139‑1 (MQ=255)                                                                                                                        cGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctg  <  1:464509/139‑1 (MQ=255)                                                                                                                         gCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGc                                                >  1:917161/1‑92 (MQ=255)                                                                                                                         gCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGc                                                <  2:917161/92‑1 (MQ=255)                                                                                                                                   |                                                                                                                                AGTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTG  >  NC_000913/49637‑49894

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