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%	161.1 / NA	45	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 (24/21);  total (24/21)

TGACAGGGAAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTTCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGCGAC  >  NC_000913/850259‑850523                                                                                                                                   |                                                                                                                                       tgaCAGGGAAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTaa                                                                                                                                >  1:1519105/1‑139 (MQ=255)     aGGGAAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAAtttt                                                                                                                            >  2:246229/1‑139 (MQ=255)         aaaaCATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        <  2:544429/139‑1 (MQ=255)          aaaCATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGata                                                                                                                       >  1:809635/1‑139 (MQ=255)                gCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTa                                                                                                                 <  2:1519105/139‑1 (MQ=255)                       aCTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTa                                                                                                          <  2:809635/139‑1 (MQ=255)                          tACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGa                                                                                                       <  1:246229/139‑1 (MQ=255)                           aCGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGAc                                                                                                      >  2:458934/1‑139 (MQ=255)                           aCGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGACCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGAc                                                                                                      >  2:1102695/1‑139 (MQ=255)                                gTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACtt                                                                                                           >  2:45193/1‑129 (MQ=255)                                gTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACtt                                                                                                           <  1:45193/129‑1 (MQ=255)                                 ttCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTAttt                                                                                                >  1:1102742/1‑139 (MQ=255)                                     ctGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  2:1656965/139‑1 (MQ=255)                                                aTTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAgg                                                                                 >  2:1024704/1‑139 (MQ=255)                                                                aTAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAgg                                                                                 <  1:564517/123‑1 (MQ=255)                                                                aTAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAgg                                                                                 >  2:564517/1‑123 (MQ=255)                                                                    cGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCa                                                             <  1:1102695/139‑1 (MQ=255)                                                                      gtgtATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACt                                                                                                            >  1:427622/1‑90 (MQ=255)                                                                      gtgtATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACt                                                                                                            <  2:427622/90‑1 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCt                                                        >  1:986435/1‑139 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCt                                                        >  1:1544002/1‑139 (MQ=255)                                                                         tATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCt                                                        >  2:1258715/1‑139 (MQ=255)                                                                              aaaCGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        <  1:634468/70‑1 (MQ=255)                                                                              aaaCGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        >  2:634468/1‑70 (MQ=255)                                                                                 cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAAttt                                                                                                                             >  1:106510/1‑62 (MQ=255)                                                                                 cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAAttt                                                                                                                             <  2:106510/62‑1 (MQ=255)                                                                                 cGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCAc                                                >  1:1300505/1‑139 (MQ=255)                                                                                                  tGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        >  2:6809/1‑50 (MQ=255)                                                                                                  tGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGat                                                                                                                        <  1:6809/50‑1 (MQ=255)                                                                                                  tGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            >  2:1471344/1‑78 (MQ=255)                                                                                                  tGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  1:1471344/78‑1 (MQ=255)                                                                                                  tGGCCCGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCt                               >  1:355934/1‑139 (MQ=255)                                                                                                     cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            >  2:78099/1‑75 (MQ=255)                                                                                                     cccGTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  1:78099/75‑1 (MQ=255)                                                                                                     cccGTTAGACATCACAAATCGCGAAGAGTTCACCATTAATTTTTGATATATTTAAAACTTAGGACTTAATTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTAGCCGCAGTTCTGGCTTTc                            >  1:957166/1‑139 (MQ=255)                                                                                                        gTTAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCAcc                         <  1:1024704/139‑1 (MQ=255)                                                                                                          tAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAgg                                                                                 >  2:478437/1‑81 (MQ=255)                                                                                                          tAGACATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAgg                                                                                 <  1:478437/81‑1 (MQ=255)                                                                                                             aCATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAgg                    <  2:1544002/139‑1 (MQ=255)                                                                                                              cATCACAAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGt                   <  2:1300505/139‑1 (MQ=255)                                                                                                                   caAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            >  2:1547447/1‑61 (MQ=255)                                                                                                                   caAATCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAAt                                                                                            <  1:1547447/61‑1 (MQ=255)                                                                                                                      aTCGCGAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGt           <  1:1258715/139‑1 (MQ=255)                                                                                                                         gcgAAGAGTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGc        >  2:1426508/1‑139 (MQ=255)                                                                                                                               agTTCCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTCTGAAAAAAATTGCATGTCTTTCAGCACTGGCTGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGCGat  <  2:1102742/139‑2 (MQ=255)                                                                                                                                   |                                                                                                                                       TGACAGGGAAAACATGCGTAATACTTACGCAGTTCTCTGAAAAAGTGATTTAAATTTAGATGGATAGCGGTGTATGGAAACGTTCTGTTACATGAAATGGCCCGTTAGACATCACAAATCGCGAAGAGTTTCCCATTAATTTTTGATATATTTAAAACTTAGGACTTATTTGAATCACATTTGAGGTGGTTATGAAAAAAATTGCATGTCTTTCAGCACTGGCCGCAGTTCTGGCTTTCACCGCAGGTACTTCCGTAGCTGCGAC  >  NC_000913/850259‑850523

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