Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,765	C→T	intergenic (+134/‑58)	kefC → / → folA	K(+) : H(+) antiporter KefC/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%	153.0 / NA	41	intergenic (+134/‑58)	kefC/folA	K(+) : H(+) antiporter KefC/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/20);  total (21/20)

GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA  >  NC_000913/49703‑49826                                                               |                                                              ggAGAGAGCGTGGTCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             <  1:2507514/65‑1 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:162009/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:1741451/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCAACAGTtga                                                             >  1:482051/1‑65 (MQ=255)  gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac                                                            >  1:2449403/1‑65 (MQ=255)  gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac                                                            >  1:3256026/1‑65 (MQ=255)  gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac                                                            >  1:220426/1‑65 (MQ=255)  gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac                                                            >  1:1538292/1‑65 (MQ=255)     agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTgacgac                                                         >  1:2975860/1‑65 (MQ=255)     agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTgacgac                                                         >  1:3424305/1‑65 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:3035919/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2370480/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2081681/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:1807195/1‑65 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:625399/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:665828/65‑1 (MQ=255)        gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:1635120/65‑1 (MQ=255)        gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:1783504/65‑1 (MQ=255)        gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:2433497/65‑1 (MQ=255)         cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGtt                                                     >  1:1525435/1‑65 (MQ=255)         cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACAGCATCCTCGCACCAGTTGACGACGGtt                                                     >  1:102983/1‑65 (MQ=255)             gACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACg                                                 >  1:154292/1‑65 (MQ=255)              aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                <  1:1994364/65‑1 (MQ=255)              aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                <  1:693811/65‑1 (MQ=255)              aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:12624/1‑65 (MQ=255)              aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:710782/1‑65 (MQ=255)              aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTATACGc                                                >  1:1266134/1‑65 (MQ=255)               cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                               >  1:994795/1‑65 (MQ=255)                tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                              <  1:1814046/65‑1 (MQ=255)                  gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTa                                            <  1:744249/65‑1 (MQ=255)                   ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:3156022/65‑1 (MQ=255)                   ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:1000419/65‑1 (MQ=255)                   ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:542345/65‑1 (MQ=255)                   ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:2375779/65‑1 (MQ=255)                     agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGt                                         <  1:833997/65‑1 (MQ=255)                     agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGt                                         <  1:665953/65‑1 (MQ=255)                     agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGt                                         <  1:704386/65‑1 (MQ=255)                      gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGta                                        >  1:754812/1‑65 (MQ=255)                        agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:3134235/1‑65 (MQ=255)                        agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:848045/1‑65 (MQ=255)                        agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:1570878/1‑65 (MQ=255)                        agAATATAAAATTTTCCTCAACATCATCCACGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:763191/1‑65 (MQ=255)                                 aaTTTTCCTCAACATCATCCTCGCACAAGTTGACGACGGGTTACGCTTTACGTATAGTGGCGACa                             >  1:511417/1‑65 (MQ=255)                                  aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:61025/1‑65 (MQ=255)                                  aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:3448265/1‑65 (MQ=255)                                                   ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa           >  1:1012803/1‑65 (MQ=255)                                                   ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa           >  1:1406770/1‑65 (MQ=255)                                                    cTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt          >  1:587898/1‑65 (MQ=255)                                                         aCCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCaa     >  1:1313607/1‑65 (MQ=255)                                                          ccAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAAt    >  1:610598/1‑65 (MQ=255)                                                            aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa  <  1:1778105/65‑1 (MQ=255)                                                               |                                                              GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA  >  NC_000913/49703‑49826

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