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%	132.0 / NA	37	intergenic (+136/‑56)	kefC/folA	potassium:proton antiporter/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/17);  total (20/17)

AGAGTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGA  >  NC_000913/49635‑49872                                                                                                                                     |                                                                                                          agagTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACgg                                                                                                     >  1:589727/1‑139 (MQ=255)  gagTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGt                                                                                                    >  2:1060965/1‑139 (MQ=255)     tGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                 <  1:1060965/139‑1 (MQ=255)     tGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                 >  1:316556/1‑139 (MQ=255)      gACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTAc                                                                                                <  1:1006253/139‑1 (MQ=255)            aaTCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                             <  2:629287/136‑1 (MQ=255)            aaTCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                             >  1:629287/1‑136 (MQ=255)                 cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTAc                                                                                                <  1:27661/128‑1 (MQ=255)                 cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTAc                                                                                                >  2:27661/1‑128 (MQ=255)                 cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACg                                                                                               <  2:936156/129‑1 (MQ=255)                 cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACg                                                                                               >  1:936156/1‑129 (MQ=255)                              aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                 <  2:1078486/114‑1 (MQ=255)                              aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTa                                                                                                 >  1:1078486/1‑114 (MQ=255)                              aCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAAtt                                                                        >  2:1174705/1‑139 (MQ=255)                                tgtttgttTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                             <  1:577318/116‑1 (MQ=255)                                tgtttgttTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCt                                                                                             >  2:577318/1‑116 (MQ=255)                                tgtttgttTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAAtttt                                                                      >  1:937353/1‑139 (MQ=255)                                      tttttgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTAt                                                                <  1:1174705/139‑1 (MQ=255)                                                 ttGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc                                                     <  2:316556/139‑1 (MQ=255)                                                         cggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACa                                                                           <  2:211966/109‑1 (MQ=255)                                                         cggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACa                                                                           >  1:211966/1‑109 (MQ=255)                                                         cggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCa                                             >  2:1105798/1‑139 (MQ=255)                                                           gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAt                                               >  1:795079/1‑135 (MQ=255)                                                           gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGAt                                               <  2:795079/135‑1 (MQ=255)                                                                    gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGTACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg                                  <  1:110405/139‑1 (MQ=255)                                                                        gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGt                              <  2:937353/139‑1 (MQ=255)                                                                          gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTAc                                                                                                <  2:432445/71‑1 (MQ=255)                                                                          gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTAc                                                                                                >  1:432445/1‑71 (MQ=255)                                                                                  aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGa                    >  2:650786/1‑139 (MQ=255)                                                                                     cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATcg                 >  2:353096/1‑139 (MQ=255)                                                                                                tataAAATTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCa      >  2:228384/1‑139 (MQ=255)                                                                                                        tttCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCttt                                                                                           >  1:524480/1‑46 (MQ=255)                                                                                                        tttCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCttt                                                                                           <  2:524480/46‑1 (MQ=255)                                                                                                          tCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                           >  2:275610/1‑76 (MQ=255)                                                                                                          tCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                                           <  1:275610/76‑1 (MQ=255)                                                                                                                           gCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa  <  1:864285/116‑1 (MQ=255)                                                                                                                           gCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGa  >  2:864285/1‑116 (MQ=255)                                                                                                                                     |                                                                                                          AGAGTGACGTAAATCACACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGA  >  NC_000913/49635‑49872

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