Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,059,964	C→T	*317* (TAG→TAA)	cbl ←	ssuEADCB/tauABCD operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,059,964	0	C	T	100.0%	152.4 / NA	47	*317* (TAG→TAA)	cbl	ssuEADCB/tauABCD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/23);  total (24/23)

TGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTC  >  NC_000913/2059836‑2060093                                                                                                                                 |                                                                                                                                  tgCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGAt                                                                                                                         <  2:1031159/139‑1 (MQ=255) tgCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGAt                                                                                                                         <  2:880922/139‑1 (MQ=255)   ccTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                >  1:175023/1‑130 (MQ=255)   ccTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                <  2:175023/130‑1 (MQ=255)   ccTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATaa                                                                                                                       >  1:604779/1‑139 (MQ=255)    cTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                      <  1:655443/139‑1 (MQ=255)              aTTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCctc                                                                                                            >  2:208152/1‑139 (MQ=255)                  ctctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGATCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATACTCAATTTCCTCTTCa                                                                                                        <  2:705106/139‑1 (MQ=255)                        aGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCtt                                                                                                  <  2:644827/139‑1 (MQ=255)                          ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTg                                                                                                     <  1:214742/134‑1 (MQ=255)                          ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTg                                                                                                     >  2:214742/1‑134 (MQ=255)                            cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCa                                                                                              >  2:396339/1‑139 (MQ=255)                               tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                           >  1:564833/1‑139 (MQ=255)                                cGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAc                                                                                          <  1:396339/139‑1 (MQ=255)                                 ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAcc                                                                                         >  2:873875/1‑139 (MQ=255)                                 ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAcc                                                                                         >  2:24826/1‑139 (MQ=255)                                 ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAcc                                                                                         >  1:999650/1‑139 (MQ=255)                                   tAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTg                                                                                       >  1:204216/1‑139 (MQ=255)                                   tAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTg                                                                                       >  2:137815/1‑139 (MQ=255)                                      aaCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCttt                                                                                                 <  1:813531/126‑1 (MQ=255)                                      aaCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCttt                                                                                                 >  2:813531/1‑126 (MQ=255)                                           cggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGa                                                                               >  1:741833/1‑139 (MQ=255)                                                 gTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCt                                                                         <  1:91540/139‑1 (MQ=255)                                                        cGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                  <  1:97598/139‑1 (MQ=255)                                                        cGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                  <  2:564833/139‑1 (MQ=255)                                                             gTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                             <  2:122126/139‑1 (MQ=255)                                                               cACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCt                                                           <  2:204216/139‑1 (MQ=255)                                                                 cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGc                                                         >  2:188822/1‑139 (MQ=255)                                                                       cGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGAt                                                                                                                         <  2:196351/69‑1 (MQ=255)                                                                       cGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGAt                                                                                                                         >  1:196351/1‑69 (MQ=255)                                                                          gTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTTACTGCTTTCCATCACCt                                                                                        <  1:375093/99‑1 (MQ=255)                                                                          gTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTTACTGCTTTCCATCACCt                                                                                        >  2:375093/1‑99 (MQ=255)                                                                                   gaggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTc                                                                                                               >  1:859195/1‑67 (MQ=255)                                                                                   gaggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTc                                                                                                               <  2:859195/67‑1 (MQ=255)                                                                                      gagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcg                                    <  2:741833/139‑1 (MQ=255)                                                                                      gagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcg                                    <  1:208152/139‑1 (MQ=255)                                                                                                 cTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt                         <  1:188822/139‑1 (MQ=255)                                                                                                   gAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTa                                                     <  2:232669/109‑1 (MQ=255)                                                                                                   gAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTa                                                     >  1:232669/1‑109 (MQ=255)                                                                                                                  ataGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCt        <  2:999650/139‑1 (MQ=255)                                                                                                                       gggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGcc                                  <  2:521889/108‑1 (MQ=255)                                                                                                                       gggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGcc                                  >  1:521889/1‑108 (MQ=255)                                                                                                                        ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTc  >  2:435989/1‑139 (MQ=255)                                                                                                                        ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTc  >  2:158382/1‑139 (MQ=255)                                                                                                                        ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTc  >  1:871464/1‑139 (MQ=255)                                                                                                                        ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTc  >  2:817647/1‑139 (MQ=255)                                                                                                                        ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTc  >  1:595819/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                  TGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTC  >  NC_000913/2059836‑2060093

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: ‑