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%	188.9 / NA	50	intergenic (+134/‑58)	kefC/folA	K(+) : H(+) antiporter KefC/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/32);  total (18/32)

GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA  >  NC_000913/49703‑49826                                                               |                                                              ggAGAGAGCGTGGACTCGCCAGCAGAATATAAACTCTTCCTCAAAATCATGCTAGCACCAGTcga                                                             >  1:3782970/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             <  1:2487943/65‑1 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:626622/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:1259157/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:1720670/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:3648358/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:3637928/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACAAGTtga                                                             >  1:3890993/1‑65 (MQ=255)  gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac                                                            >  1:2316046/1‑65 (MQ=255)  gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac                                                            >  1:1821474/1‑65 (MQ=255)  gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac                                                            >  1:2922314/1‑65 (MQ=255)   agagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacg                                                           <  1:1130243/65‑1 (MQ=255)       cgCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:1436083/2‑65 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:3063598/1‑65 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:3210171/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2763846/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:665303/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2384441/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2636637/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1572247/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1638016/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1658105/65‑1 (MQ=255)       agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGATgacgactg                                                       >  1:1193494/1‑63 (MQ=255)        gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:970717/65‑1 (MQ=255)        gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:3211907/65‑1 (MQ=255)        gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:2941220/65‑1 (MQ=255)        gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:2552295/65‑1 (MQ=255)        gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:1787682/65‑1 (MQ=255)         cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGtt                                                     >  1:2059132/1‑65 (MQ=255)            ggACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                  <  1:3141869/65‑1 (MQ=255)            ggACTCCCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                  <  1:1472667/65‑1 (MQ=255)              aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                <  1:355648/65‑1 (MQ=255)              aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                <  1:3550449/65‑1 (MQ=255)              aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACAAGCTGACGACGGTTTACGc                                                >  1:3407197/1‑65 (MQ=255)               cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                               >  1:2331223/1‑65 (MQ=255)               cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                               <  1:3290952/65‑1 (MQ=255)                tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                              <  1:1513835/65‑1 (MQ=255)                  gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTa                                            <  1:177553/65‑1 (MQ=255)                  gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTa                                            <  1:1250762/65‑1 (MQ=255)                   ccAGCAGAATGTAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:2444452/65‑1 (MQ=255)                   ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:3845091/65‑1 (MQ=255)                   ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           >  1:3302647/1‑65 (MQ=255)                   ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:2740236/65‑1 (MQ=255)                   ccAGCAGAATATAAAATTTTCCTCAACACCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:1881615/65‑1 (MQ=255)                    cagcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACg                                          <  1:967895/65‑1 (MQ=255)                     agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGt                                         <  1:437879/65‑1 (MQ=255)                      gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtc                                        >  1:2019889/1‑64 (MQ=255)                        agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:878226/1‑65 (MQ=255)                        agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:1558190/1‑65 (MQ=255)                        agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:3641578/1‑65 (MQ=255)                        agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:2286502/1‑65 (MQ=255)                         gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAg                                     <  1:2301675/65‑1 (MQ=255)                           atatAAAATTTTCCTCAACATCATCCTCGCCCCAGTTGACGACGGTTAACGCTTTACGTATAGTg                                   >  1:2264658/1‑65 (MQ=255)                                  aTTTTCCTCAACATCATCGTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:752151/1‑65 (MQ=255)                                  aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:2123198/1‑65 (MQ=255)                                  aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:3995232/1‑65 (MQ=255)                                              atcatcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCg                >  1:1537697/1‑65 (MQ=255)                                                catcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCggg              <  1:1993499/65‑1 (MQ=255)                                                 atcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATTGGGa             >  1:3807907/1‑65 (MQ=255)                                                   ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACCATTTTTTTTATCGGGaaa           >  1:1971196/1‑65 (MQ=255)                                                   ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa           >  1:2708160/1‑65 (MQ=255)                                                   ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa           >  1:556799/1‑65 (MQ=255)                                                   ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa           >  1:2339981/1‑65 (MQ=255)                                                    cTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt          >  1:433571/1‑65 (MQ=255)                                                          ccAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAAt    <  1:1009658/65‑1 (MQ=255)                                                           cAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATg   <  1:2240207/65‑1 (MQ=255)                                                           cAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATg   <  1:2700554/65‑1 (MQ=255)                                                            aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa  <  1:97072/65‑1 (MQ=255)                                                            aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa  <  1:2074006/65‑1 (MQ=255)                                                            aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa  <  1:3963253/65‑1 (MQ=255)                                                               |                                                              GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA  >  NC_000913/49703‑49826

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