Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,709,163	T→C	S8P (TCC→CCC)	rsxD →	SoxR iron‑sulfur cluster reduction factor component, putative membrane protein of electron transport complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,709,163	0	T	C	100.0%	144.9 / NA	41	S8P (TCC→CCC)	rsxD	SoxR iron‑sulfur cluster reduction factor component, putative membrane protein of electron transport complex
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/17);  total (24/17)

GGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCGGTT  >  NC_000913/1709027‑1709297                                                                                                                                         |                                                                                                                                       ggCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGcccc                                                                                                                                      >  1:1076724/1‑139 (MQ=255)      aTGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                 >  1:554362/1‑139 (MQ=255)          ggTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATAcc                                                                                                                             >  2:258767/1‑139 (MQ=255)                 gaagaaCAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAcc                                                                                                                      >  2:1153156/1‑139 (MQ=255)                         ggTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAg                                                                                                              <  2:797102/139‑1 (MQ=255)                                  gcgcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                       >  2:1309874/1‑121 (MQ=255)                                  gcgcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                       <  1:1309874/121‑1 (MQ=255)                                  gcgcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCgcgc                                                                                                     >  1:549903/1‑139 (MQ=255)                                   cgcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCa                                                                                                    <  2:835723/139‑1 (MQ=255)                                       aaGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTAt                                                                                                <  1:1153156/139‑1 (MQ=255)                                         gCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGc                                                                                              >  1:575471/1‑139 (MQ=255)                                                  gCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATAcc                                                                                                                             <  1:363164/99‑1 (MQ=255)                                                  gCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATAcc                                                                                                                             >  2:363164/1‑99 (MQ=255)                                                                 cgcgCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGc                                                                      >  2:1042967/1‑139 (MQ=255)                                                                              aaaaaGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACa                                                                                                           >  1:253941/1‑89 (MQ=255)                                                                              aaaaaGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACa                                                                                                           <  2:253941/89‑1 (MQ=255)                                                                               aaaaGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcg                                                        <  1:258767/139‑1 (MQ=255)                                                                                    ccgccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACt                                                   >  1:213614/1‑139 (MQ=255)                                                                                    ccgccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACt                                                   >  1:1331902/1‑139 (MQ=255)                                                                                    ccgccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACt                                                   >  1:737882/1‑139 (MQ=255)                                                                                      gccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGgtg                                                                                      <  1:1074047/102‑1 (MQ=255)                                                                                      gccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGgtg                                                                                      >  2:1074047/1‑102 (MQ=255)                                                                                                 ggTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTttggt                                      <  2:1331902/139‑1 (MQ=255)                                                                                                 ggTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTttggt                                      <  2:737882/139‑1 (MQ=255)                                                                                                  gTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcg                                                        <  2:346871/120‑1 (MQ=255)                                                                                                  gTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcg                                                        >  1:346871/1‑120 (MQ=255)                                                                                                      tAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTgggg                                 <  2:549903/139‑1 (MQ=255)                                                                                                             ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACtctc                          >  1:723503/1‑139 (MQ=255)                                                                                                             ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACtctc                          >  1:527692/1‑139 (MQ=255)                                                                                                                 tAAATGGTATTCAGTATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTActc                            >  1:869665/1‑133 (MQ=255)                                                                                                                 tAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTActc                            <  2:869665/133‑1 (MQ=255)                                                                                                                   aaTGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAg                    >  1:144841/1‑139 (MQ=255)                                                                                                                         aTTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTg              >  1:273261/1‑139 (MQ=255)                                                                                                                          ttCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGt             <  2:144841/139‑1 (MQ=255)                                                                                                                           tCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTg                                                                                            >  1:8106/1‑59 (MQ=255)                                                                                                                           tCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTg                                                                                            <  2:8106/59‑1 (MQ=255)                                                                                                                           tCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAg                                                          >  2:100265/1‑93 (MQ=255)                                                                                                                           tCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAg                                                          <  1:100265/93‑1 (MQ=255)                                                                                                                             aGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTa                                                                                                 >  1:570896/1‑52 (MQ=255)                                                                                                                             aGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTa                                                                                                 <  2:570896/52‑1 (MQ=255)                                                                                                                             aGAATACCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCCCATTATGCTGTTGGTGTTGCTAGCCGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGTGTACTCTCGTCCAGATCCTGTTgc          >  2:1028289/1‑138 (MQ=255)                                                                                                                              gAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGc         >  2:1106612/1‑139 (MQ=255)                                                                                                                                    ctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGt             >  2:1317184/1‑129 (MQ=255)                                                                                                                                    ctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGt             <  1:1317184/129‑1 (MQ=255)                                                                                                                                    ctagccCCCCTTATACCCATAACCA‑CGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCGGtt  >  2:388464/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       GGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCGGTT  >  NC_000913/1709027‑1709297

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