Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	850,389	T→C	intergenic (‑292/‑61)	rhtA ← / → ompX	threonine and homoserine efflux system/outer membrane protein X

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	850,389	0	T	C	100.0%	201.8 / NA	58	intergenic (‑292/‑61)	rhtA/ompX	threonine and homoserine efflux system/outer membrane protein X
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/28);  total (30/28)

AAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTTCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGCGAC  >  NC_000913/850267‑850523                                                                                                                           |                                                                                                                                       aaaaCATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        <  1:687184/139‑1 (MQ=255) aaaaCATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        <  2:1325255/139‑1 (MQ=255) aaaaCATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        >  2:687184/1‑139 (MQ=255) aaaaCATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        >  1:411136/1‑139 (MQ=255)         cGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        >  1:1374963/1‑131 (MQ=255)         cGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        <  2:1374963/131‑1 (MQ=255)                  tACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGa                                                                                                       <  2:411136/139‑1 (MQ=255)                     gCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACtt                                                                                                    >  2:962928/1‑139 (MQ=255)                         ttCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        <  1:174356/115‑1 (MQ=255)                         ttCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        >  2:174356/1‑115 (MQ=255)                         ttCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTAttt                                                                                                >  1:594321/1‑139 (MQ=255)                         ttCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTAttt                                                                                                >  2:758912/1‑139 (MQ=255)                           ctctGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                                 >  2:1294480/1‑120 (MQ=255)                           ctctGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                                 <  1:1294480/120‑1 (MQ=255)                             ctGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATTGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTTATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  1:724365/139‑1 (MQ=255)                             ctGAAAAAGTGATTTAAATCTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTCTTTGAAt                                                                                            <  1:962928/139‑1 (MQ=255)                              tGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATc                                                                                           >  2:305510/1‑139 (MQ=255)                                        aTTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAgg                                                                                 >  1:272304/1‑139 (MQ=255)                                                  aGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGaa                                                                       <  2:1119176/139‑1 (MQ=255)                                                        aTAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        >  2:1465535/1‑84 (MQ=255)                                                        aTAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        <  1:1465535/84‑1 (MQ=255)                                                             ggTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  1:315329/107‑1 (MQ=255)                                                             ggTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATAGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            >  2:315329/1‑107 (MQ=255)                                                              gtgtATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATg                                                           <  2:272304/139‑1 (MQ=255)                                                                 tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCAtt                                                                                                                                  <  2:1101191/65‑1 (MQ=255)                                                                 tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCAtt                                                                                                                                  >  1:1101191/1‑65 (MQ=255)                                                                 tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            >  2:703423/1‑103 (MQ=255)                                                                 tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  1:703423/103‑1 (MQ=255)                                                                      aaaCGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACtt                                                                                                           >  1:124673/1‑83 (MQ=255)                                                                      aaaCGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACtt                                                                                                           <  2:124673/83‑1 (MQ=255)                                                                         cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                                 >  2:56153/1‑74 (MQ=255)                                                                         cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                                 <  1:56153/74‑1 (MQ=255)                                                                         cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCAc                                                >  2:689058/1‑139 (MQ=255)                                                                           ttCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  1:152436/93‑1 (MQ=255)                                                                           ttCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            >  2:152436/1‑93 (MQ=255)                                                                           ttCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACAtt                                                                                      <  2:340646/99‑1 (MQ=255)                                                                           ttCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACAtt                                                                                      >  1:340646/1‑99 (MQ=255)                                                                           ttCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTg                                              >  1:431585/1‑139 (MQ=255)                                                                             cTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGc                                            >  1:1384193/1‑139 (MQ=255)                                                                             cTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGc                                            >  2:1017214/1‑139 (MQ=255)                                                                                        aaTGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            >  2:850803/1‑80 (MQ=255)                                                                                        aaTGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  1:850803/80‑1 (MQ=255)                                                                                             cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                                 >  2:660398/1‑54 (MQ=255)                                                                                             cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                                 <  1:660398/54‑1 (MQ=255)                                                                                             cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAgg                                                                                 <  1:737997/86‑1 (MQ=255)                                                                                             cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAgg                                                                                 >  2:737997/1‑86 (MQ=255)                                                                                                gTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCAcc                         <  1:689058/139‑1 (MQ=255)                                                                                                       aTCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATg                                                           >  2:427691/1‑98 (MQ=255)                                                                                                       aTCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATg                                                           <  1:427691/98‑1 (MQ=255)                                                                                                          acaAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAgg                    >  2:494757/1‑134 (MQ=255)                                                                                                          acaAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAgg                    <  1:494757/134‑1 (MQ=255)                                                                                                           caAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  2:360621/61‑1 (MQ=255)                                                                                                           caAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            >  1:360621/1‑61 (MQ=255)                                                                                                                   gAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGc                                                              <  2:55798/83‑1 (MQ=255)                                                                                                                   gAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGc                                                              >  1:55798/1‑83 (MQ=255)                                                                                                                    aaGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATg                                                           <  1:1390921/85‑1 (MQ=255)                                                                                                                    aaGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATg                                                           >  2:1390921/1‑85 (MQ=255)                                                                                                                       agTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGCGAc  <  2:431585/139‑1 (MQ=255)                                                                                                                           |                                                                                                                                       AAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTTCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGCGAC  >  NC_000913/850267‑850523

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