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%	147.6 / NA	41	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 (20/21);  total (20/21)

TGACAGGGAAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTTCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGC  >  NC_000913/850259‑850520                                                                                                                                   |                                                                                                                                    tgaCAGGGAAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTaa                                                                                                                             >  2:180770/1‑139 (MQ=255)         aaaaCATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                     <  2:1549698/139‑1 (MQ=255)                 cGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTaa                                                                                                             >  2:395554/1‑139 (MQ=255)                          tACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGa                                                                                                    <  1:1095720/139‑1 (MQ=255)                           aCGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                     >  2:1390976/1‑121 (MQ=255)                           aCGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                     <  1:1390976/121‑1 (MQ=255)                           aCGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                              <  2:327560/128‑1 (MQ=255)                           aCGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                              >  1:327560/1‑128 (MQ=255)                                     ctGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                         <  1:180770/139‑1 (MQ=255)                                         aaaaGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATcaca                                                                                     <  1:854068/139‑1 (MQ=255)                                           aaGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACAtt                                                                                   <  1:395554/139‑1 (MQ=255)                                                               gATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                              >  2:755759/1‑92 (MQ=255)                                                               gATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                              <  1:755759/92‑1 (MQ=255)                                                               gATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                         <  1:1402691/113‑1 (MQ=255)                                                               gATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                         >  2:1402691/1‑113 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTa                                                                                                                              <  2:177647/66‑1 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTa                                                                                                                              >  1:177647/1‑66 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                         >  1:1436684/1‑103 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                         <  2:1436684/103‑1 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                         <  2:1514109/103‑1 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                         >  1:1514109/1‑103 (MQ=255)                                                                                 cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACt                                                                                                         <  1:708647/79‑1 (MQ=255)                                                                                 cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACt                                                                                                         >  2:708647/1‑79 (MQ=255)                                                                                 cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCAc                                             >  2:827702/1‑139 (MQ=255)                                                                                   ttCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTg                                           >  2:1474473/1‑139 (MQ=255)                                                                                       gTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCg                                       <  1:784700/139‑1 (MQ=255)                                                                                                aaTGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCa                                                          <  1:612463/111‑1 (MQ=255)                                                                                                aaTGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCa                                                          >  2:612463/1‑111 (MQ=255)                                                                                                     cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                              <  2:726411/54‑1 (MQ=255)                                                                                                     cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                              >  1:726411/1‑54 (MQ=255)                                                                                                     cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGa                                                                                                    >  2:601506/1‑64 (MQ=255)                                                                                                     cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGa                                                                                                    <  1:601506/64‑1 (MQ=255)                                                                                                     cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTc                         >  1:1720669/1‑139 (MQ=255)                                                                                                              cATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGtt                                  <  2:898925/121‑1 (MQ=255)                                                                                                              cATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGtt                                  >  1:898925/1‑121 (MQ=255)                                                                                                                 cacaAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTaa                                                                                                             >  2:1641503/1‑43 (MQ=255)                                                                                                                 cacaAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTaa                                                                                                             <  1:1641503/43‑1 (MQ=255)                                                                                                                      aTCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATg                                                        >  2:1377709/1‑91 (MQ=255)                                                                                                                      aTCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATg                                                        <  1:1377709/91‑1 (MQ=255)                                                                                                                          cgAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCt    >  1:424795/1‑139 (MQ=255)                                                                                                                            aaGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGc  <  1:827702/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                    TGACAGGGAAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTTCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGC  >  NC_000913/850259‑850520

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