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%	133.6 / NA	42	*317* (TAG→TAA)	cbl	ssuEADCB/tauABCD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/20);  total (22/20)

TGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGTCCC  >  NC_000913/2059856‑2060099                                                                                                             |                                                                                                                                        tGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTg                                                                                                           <  1:810536/139‑1 (MQ=255)  gTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGc                                                                                                          <  1:907812/139‑1 (MQ=255)     gTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCttt                                                                                                       >  2:1107863/1‑139 (MQ=255)      ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGc                                                                                                          <  1:506529/135‑1 (MQ=255)      ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGc                                                                                                          >  2:506529/1‑135 (MQ=255)      ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTc                                                                                                      >  1:40935/1‑139 (MQ=255)           tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                      <  1:726489/102‑1 (MQ=255)           tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTtat                                                                                                                                      >  2:726489/1‑102 (MQ=255)           tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTc                                                                                                      >  1:474372/1‑134 (MQ=255)           tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTc                                                                                                      <  2:474372/134‑1 (MQ=255)           tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                                 >  1:764218/1‑139 (MQ=255)           tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                                 >  1:781541/1‑139 (MQ=255)           tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                                 >  2:100091/1‑139 (MQ=255)             ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGAt                                                                                                                               >  2:500082/1‑107 (MQ=255)             ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGAt                                                                                                                               <  1:500082/107‑1 (MQ=255)                         gACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGata                                                                                   >  1:1356571/1‑139 (MQ=255)                                            aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTc                                                                                                                     >  2:1359948/1‑86 (MQ=255)                                            aCTGGTTCGAGTCCAGACAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTc                                                                                                                     <  1:1359948/86‑1 (MQ=255)                                                gTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCAtt                                                            >  1:1437010/1‑139 (MQ=255)                                                gTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGAGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCAtt                                                            >  2:783401/1‑139 (MQ=255)                                                   cGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                         <  1:1107863/139‑1 (MQ=255)                                                           gTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                         >  2:466926/1‑131 (MQ=255)                                                           gTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                         <  1:466926/131‑1 (MQ=255)                                                           gTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCa                                                 >  1:385208/1‑139 (MQ=255)                                                                aggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTg                                                                                                           >  2:864451/1‑76 (MQ=255)                                                                aggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTg                                                                                                           <  1:864451/76‑1 (MQ=255)                                                                   agCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                                   >  2:1140996/1‑113 (MQ=255)                                                                   agCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                                   <  1:1140996/113‑1 (MQ=255)                                                                    gCCAAATTCTTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcg                                          >  2:471911/1‑137 (MQ=255)                                                                    gCCAAATTCTTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcg                                          <  1:471911/137‑1 (MQ=255)                                                                    gCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGcc                                        <  1:100091/139‑1 (MQ=255)                                                                     ccAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcgc                                         <  2:1445034/137‑1 (MQ=255)                                                                     ccAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcgc                                         >  1:1445034/1‑137 (MQ=255)                                                                                 agtgCCGGCTTTTATAGCGGTATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCGGCTTGATTTCCTCAACCGACAGTCCTGCATTACAAAGTTCCTGAAAGTGCCAGACATAGTTACg                           <  1:166875/136‑1 (MQ=255)                                                                                                    ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGcc                                                                                           >  1:691415/1‑56 (MQ=255)                                                                                                    ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGcc                                                                                           <  2:691415/56‑1 (MQ=255)                                                                                                    ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                         >  2:1072686/1‑90 (MQ=255)                                                                                                    ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                         <  1:1072686/90‑1 (MQ=255)                                                                                                    ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt                               >  2:43009/1‑116 (MQ=255)                                                                                                    ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt                               <  1:43009/116‑1 (MQ=255)                                                                                                       tttttGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGt     <  2:1356571/139‑1 (MQ=255)                                                                                                          ttGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGTccc  <  2:764218/139‑1 (MQ=255)                                                                                                             |                                                                                                                                        TGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGTCCC  >  NC_000913/2059856‑2060099

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