Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	883,647	G→T	intergenic (‑259/‑26)	ybjG ← / → mdfA	undecaprenyl pyrophosphate phosphatase/multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	883,647	0	G	T	100.0%	211.1 / NA	56	intergenic (‑259/‑26)	ybjG/mdfA	undecaprenyl pyrophosphate phosphatase/multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (35/21);  total (35/21)

TTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCGCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTGGACG  >  NC_000913/883585‑883710                                                               |                                                                ttAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGctcg                                                               <  1:924498/65‑1 (MQ=255) ttAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGctcg                                                               <  1:3303185/65‑1 (MQ=255) ttAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGctcg                                                               <  1:3370324/65‑1 (MQ=255)   aaCCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGat                                                             <  1:43093/65‑1 (MQ=255)    aCCGTGGTTTCAGCTGCAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGatt                                                            >  1:1534277/1‑65 (MQ=255)    aCCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGatt                                                            >  1:371234/1‑65 (MQ=255)    aCCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGatt                                                            >  1:1291571/1‑65 (MQ=255)    aCCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGatt                                                            >  1:3352966/1‑65 (MQ=255)    aCCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGatt                                                            >  1:1254050/1‑65 (MQ=255)    aCCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGatt                                                            >  1:119249/1‑65 (MQ=255)    aCCGTGGTTTCAGCTACAAAATTCGCATTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGatt                                                            >  1:201840/1‑65 (MQ=255)    aCCGTGGTTTCACCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGatt                                                            >  1:2257945/1‑65 (MQ=255)      cGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGattat                                                          >  1:39794/1‑65 (MQ=255)      cGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGattat                                                          >  1:1575730/1‑65 (MQ=255)           tttCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGattattattc                                                     >  1:2949248/1‑64 (MQ=255)           tttCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGATTATTATTg                                                     >  1:3343395/1‑65 (MQ=255)           tttCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGATTATTATTg                                                     >  1:2227758/1‑65 (MQ=255)           tttCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGATTATTATTg                                                     >  1:1647128/1‑65 (MQ=255)                gCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTGGATTATTATTGGCgaa                                                <  1:1743405/65‑1 (MQ=255)                gCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCgaa                                                <  1:3302898/65‑1 (MQ=255)                gCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCgaa                                                <  1:1292456/65‑1 (MQ=255)                  tACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAATCTCTGCTCGATTATTATTGGCgaaga                                              <  1:2473720/65‑1 (MQ=255)                               ttCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCaaa                                 <  1:1208532/65‑1 (MQ=255)                               ttCTCGTTAGCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCaaa                                 <  1:2833030/65‑1 (MQ=255)                                  tcGTTAGCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAaaata                              <  1:2611674/65‑1 (MQ=255)                                        gCTGCGCTTTTATTAAGCTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        <  1:2943130/65‑1 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:49737/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:3256456/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:3192033/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        <  1:117195/65‑1 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:3036828/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:27376/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:619084/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:1426011/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:2430836/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:2376020/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:235370/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:1503470/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:516143/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:2199464/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:2068321/1‑65 (MQ=255)                                        gCTGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAg                        >  1:1854366/1‑65 (MQ=255)                                          tGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCt                      >  1:2556465/1‑65 (MQ=255)                                          tGCGCTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCt                      <  1:1367247/65‑1 (MQ=255)                                           gcgcTTTTATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCtt                     <  1:2243901/65‑1 (MQ=255)                                                tttATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGt                <  1:828809/65‑1 (MQ=255)                                                tttATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGt                <  1:1128357/65‑1 (MQ=255)                                                tttATTAAACTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGt                <  1:2465064/65‑1 (MQ=255)                                                    ttaAACTCTGCTCGATTATTATTGGCGAAGAAACTGCATGCAAAATAAATTAGCTTCCGGTGCCa            >  1:9394/1‑65 (MQ=255)                                                      aaaCTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAgg          <  1:3036206/65‑1 (MQ=255)                                                      aaaCTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAgg          <  1:629275/65‑1 (MQ=255)                                                       aaCTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGc         <  1:3088154/65‑1 (MQ=255)                                                       aaCTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGc         >  1:730365/1‑65 (MQ=255)                                                       aaCTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGc         <  1:1121637/65‑1 (MQ=255)                                                        aCTCTGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCt        <  1:1615866/65‑1 (MQ=255)                                                          tctGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTg      >  1:3166057/1‑65 (MQ=255)                                                          tctGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTg      >  1:2728879/1‑65 (MQ=255)                                                          tctGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTg      >  1:169419/1‑65 (MQ=255)                                                          tctGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTg      >  1:1424748/1‑65 (MQ=255)                                                            tGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTGGa    <  1:2906658/65‑1 (MQ=255)                                                            tGCTCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTGGa    <  1:1544961/65‑1 (MQ=255)                                                              agcgATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTGGACg  <  1:1267432/64‑1 (MQ=255)                                                               |                                                                TTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAACTCTGCGCGATTATTATTGGCGAAGAAATTGCATGCAAAATAAATTAGCTTCCGGTGCCAGGCTTGGACG  >  NC_000913/883585‑883710

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