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

TGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGT  >  NC_000913/2059836‑2060096                                                                                                                                 |                                                                                                                                     tgCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGAt                                                                                                                            >  2:573652/1‑139 (MQ=255) tgCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAATCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGAt                                                                                                                            <  1:63063/139‑1 (MQ=255)        acaacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAAt                                                                                                                     >  1:1350875/1‑139 (MQ=255)            cGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTcc                                                                                                                 >  1:1177948/1‑139 (MQ=255)            cGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTcc                                                                                                                 >  2:215701/1‑139 (MQ=255)             gATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCct                                                                                                                <  2:1350875/139‑1 (MQ=255)                tCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCtt                                                                                                             >  1:1657196/1‑139 (MQ=255)                               tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                              >  1:851447/1‑139 (MQ=255)                               tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                              >  1:552173/1‑139 (MQ=255)                                cGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAc                                                                                             <  1:573652/139‑1 (MQ=255)                                 ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAcc                                                                                            >  2:1549519/1‑139 (MQ=255)                                 ggTAGAACGGCGGACTGTTAATCCGTATGTCACGGGTTGGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGACCGCTGTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAca                                                                                            >  2:576624/1‑138 (MQ=37)                                    aGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGc                                                                                         >  1:32697/1‑139 (MQ=255)                                          gcggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTg                                                                                   <  2:1177948/139‑1 (MQ=255)                                            ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                                 <  1:215701/139‑1 (MQ=255)                                             gACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGata                                                                                >  2:636568/1‑139 (MQ=255)                                             gACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGata                                                                                >  2:334316/1‑139 (MQ=255)                                              aCTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGatat                                                                               <  1:1549519/139‑1 (MQ=255)                                                tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATcc                                                                             >  2:841311/1‑139 (MQ=255)                                                       ccGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTg                                                                      >  2:60253/1‑139 (MQ=255)                                                          tATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACa                                                                   >  2:1491180/1‑139 (MQ=255)                                                             gTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                                <  2:851447/139‑1 (MQ=255)                                                                 cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCa                                                                                                           >  1:854976/1‑92 (MQ=255)                                                                 cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCa                                                                                                           <  2:854976/92‑1 (MQ=255)                                                                    gTTCGAGTCCAGTCAGAGGAGCCAAATTTCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATCTCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCt                                                                                     >  2:607639/1‑111 (MQ=255)                                                                    gTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCt                                                                                     <  1:607639/111‑1 (MQ=255)                                                                    gTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                                >  1:1218392/1‑132 (MQ=255)                                                                    gTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                                <  2:1218392/132‑1 (MQ=255)                                                                     ttCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                         <  1:1189309/74‑1 (MQ=255)                                                                     ttCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAAt                                                                                                                         >  2:1189309/1‑74 (MQ=255)                                                                     ttCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTa                                                        <  1:636568/139‑1 (MQ=255)                                                                       cGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                      <  1:60253/139‑1 (MQ=255)                                                                        gAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACaa                                                     >  1:1458937/1‑139 (MQ=255)                                                                         aGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTg                                                                                                                              <  1:596794/65‑1 (MQ=255)                                                                         aGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTg                                                                                                                              >  2:596794/1‑65 (MQ=255)                                                                         aGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACaaa                                                    >  1:1552654/1‑139 (MQ=255)                                                                                      gagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcg                                       <  2:1552654/139‑1 (MQ=255)                                                                                         ccAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCa                                    >  1:383291/1‑139 (MQ=255)                                                                                          cAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                     >  2:1237528/1‑105 (MQ=255)                                                                                          cAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                     <  1:1237528/105‑1 (MQ=255)                                                                                             aTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACa                                >  2:1256710/1‑139 (MQ=255)                                                                                                         cccGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGc                                                                                                       <  2:355971/56‑1 (MQ=255)                                                                                                         cccGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGc                                                                                                       >  1:355971/1‑56 (MQ=255)                                                                                                           cGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGc                                                                                         <  1:1435697/68‑1 (MQ=255)                                                                                                           cGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGc                                                                                         >  2:1435697/1‑68 (MQ=255)                                                                                                           cGCTGTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAg                  <  1:334316/139‑1 (MQ=255)                                                                                                                  ataGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCt           <  2:1657196/139‑1 (MQ=255)                                                                                                                      cGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGc                                                                                         <  1:1535164/57‑1 (MQ=255)                                                                                                                      cGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGc                                                                                         >  2:1535164/1‑57 (MQ=255)                                                                                                                      cGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCt       >  1:620225/1‑139 (MQ=255)                                                                                                                        ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATc                                                                              >  2:22016/1‑66 (MQ=255)                                                                                                                        ggATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATc                                                                              <  1:22016/66‑1 (MQ=255)                                                                                                                           tttttGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGt  <  1:1491180/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                     TGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGT  >  NC_000913/2059836‑2060096

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