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%	199.1 / NA	62	intergenic (+134/‑58)	kefC/folA	potassium:proton antiporter/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/30);  total (32/30)

AGTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGC  >  NC_000913/49637‑49895                                                                                                                                 |                                                                                                                                   agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGtt                                                                                                                          >  2:1289353/1‑139 (MQ=255) agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGtt                                                                                                                          >  1:953302/1‑139 (MQ=255) agTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGATGACGACGGtt                                                                                                                          >  1:853382/1‑139 (MQ=255)     aCGTAAATCACACTTTACTGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACg                                                                                                                      <  1:883718/139‑1 (MQ=255)     aCGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGTGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACg                                                                                                                      >  2:794871/1‑139 (MQ=255)      cGTAAATCACACTTTAGAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTTCGc                                                                                                                     <  1:894106/139‑1 (MQ=255)      cGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                                                                                     >  2:255563/1‑139 (MQ=255)      cGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                                                                                     <  2:1128978/139‑1 (MQ=255)               cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                                                                                     <  2:972346/130‑1 (MQ=255)               cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                                                                                     >  1:972346/1‑130 (MQ=255)                acTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                                                                                           >  2:607234/1‑139 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                                                                                                    >  1:630527/1‑118 (MQ=255)                            aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                                                                                                    <  2:630527/118‑1 (MQ=255)                                tttgtttttgttTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                                                                                       >  2:718214/1‑111 (MQ=255)                                tttgtttttgttTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                                                                                       <  1:718214/111‑1 (MQ=255)                                        tgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGc                                                                                                      <  1:912515/120‑1 (MQ=255)                                        tgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGc                                                                                                      >  2:912515/1‑120 (MQ=255)                                        tgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg                                                                                   >  2:1076735/1‑139 (MQ=255)                                        tgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg                                                                                   >  1:249119/1‑139 (MQ=255)                                           ttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa                                                                                <  1:1076735/139‑1 (MQ=255)                                              aTTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtct                                                                             >  1:56260/1‑139 (MQ=255)                                                   aaTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATg                                                                        <  2:853382/139‑1 (MQ=255)                                                      gcggtgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                                                                     <  1:794871/139‑1 (MQ=255)                                                         gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg                                                                                   <  1:214996/122‑1 (MQ=255)                                                         gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg                                                                                   >  2:214996/1‑122 (MQ=255)                                                          cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                                                                                     >  2:367867/1‑87 (MQ=255)                                                          cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                                                                                     <  1:367867/87‑1 (MQ=255)                                                          cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAgcc                                                                 >  1:397984/1‑137 (MQ=255)                                                                    gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTgcggc                                                       <  2:916420/139‑1 (MQ=255)                                                                      gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                     <  1:255563/139‑1 (MQ=255)                                                                      gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                     <  2:953302/139‑1 (MQ=255)                                                                       agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGtt                                                    >  1:725836/1‑139 (MQ=255)                                                                       agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGtt                                                    >  2:1221306/1‑139 (MQ=255)                                                                        gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATc                                                                                      >  1:11938/1‑104 (MQ=255)                                                                        gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATc                                                                                      <  2:11938/104‑1 (MQ=255)                                                                        gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg                                                                                   <  1:1166738/107‑1 (MQ=255)                                                                        gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg                                                                                   >  2:1166738/1‑107 (MQ=255)                                                                        gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa                                                                                <  2:948880/110‑1 (MQ=255)                                                                        gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa                                                                                >  1:948880/1‑110 (MQ=255)                                                                         agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAg                                                  >  1:425800/1‑139 (MQ=255)                                                                           cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc                                                                            <  2:823432/111‑1 (MQ=255)                                                                           cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc                                                                            >  1:823432/1‑111 (MQ=255)                                                                               gACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                  >  1:570239/1‑101 (MQ=255)                                                                               gACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                                                  <  2:570239/101‑1 (MQ=255)                                                                                    gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                                                                                       <  1:344711/59‑1 (MQ=255)                                                                                    gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                                                                                       >  2:344711/1‑59 (MQ=255)                                                                                    gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                     >  2:1090905/1‑125 (MQ=255)                                                                                    gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGt                                                     <  1:1090905/125‑1 (MQ=255)                                                                                         cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                  <  1:1289353/139‑1 (MQ=255)                                                                                         cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTAt                                  >  2:250530/1‑139 (MQ=255)                                                                                             atatAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGc                              <  1:1221306/139‑1 (MQ=255)                                                                                              tataAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCa                             <  2:56260/139‑1 (MQ=255)                                                                                                   aaTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGaa                        <  2:293826/139‑1 (MQ=255)                                                                                                        tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                                                                                       <  1:789088/39‑1 (MQ=255)                                                                                                        tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                                                                                       >  2:789088/1‑39 (MQ=255)                                                                                                        tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg                                                                                   >  2:1078294/1‑75 (MQ=255)                                                                                                        tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg                                                                                   <  1:1078294/75‑1 (MQ=255)                                                                                                          cTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCa                 <  2:397984/139‑1 (MQ=255)                                                                                                           tCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCAt                <  2:676133/139‑1 (MQ=255)                                                                                                             aaCATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGc              >  1:750986/1‑139 (MQ=255)                                                                                                                    tcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGaa       >  2:314138/1‑139 (MQ=255)                                                                                                                         gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGCCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAAcctgc  >  2:1123568/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                   AGTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGC  >  NC_000913/49637‑49895

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