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%	137.6 / NA	40	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 (20/20);  total (20/20)

TAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCGGT  >  NC_000913/1709030‑1709296                                                                                                                                      |                                                                                                                                      tAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCt                                                                                                                                  >  1:591848/1‑139 (MQ=255)             agaagaACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                      <  2:586175/139‑1 (MQ=255)                    cAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGcc                                                                                                               <  1:1574739/139‑1 (MQ=255)                     aGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                      >  2:662775/1‑131 (MQ=255)                     aGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                      <  1:662775/131‑1 (MQ=255)                     aGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCa                                                                                                              >  2:995690/1‑139 (MQ=255)                               gcgcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCgcgc                                                                                                    >  1:1636148/1‑139 (MQ=255)                                cgcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                <  1:1119714/110‑1 (MQ=255)                                cgcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                >  2:1119714/1‑110 (MQ=255)                                cgcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCa                                                                                                   <  1:894081/139‑1 (MQ=255)                                    aaGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTAt                                                                                               <  2:1636148/139‑1 (MQ=255)                                          cAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                <  1:1140906/100‑1 (MQ=255)                                          cAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                >  2:1140906/1‑100 (MQ=255)                                                        aTTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAg                                                                           <  1:995690/139‑1 (MQ=255)                                                              cgcgCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGc                                                                     >  1:879365/1‑139 (MQ=255)                                                               gcgcTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGcc                                                                    <  1:560134/139‑1 (MQ=255)                                                                     aGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAg                                                                           <  1:26357/126‑1 (MQ=255)                                                                     aGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAg                                                                           >  2:26357/1‑126 (MQ=255)                                                                      ggCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATc                                                             >  1:210074/1‑139 (MQ=255)                                                                       gCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCa                                                                                                                    <  2:909316/83‑1 (MQ=255)                                                                       gCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCa                                                                                                                    >  1:909316/1‑83 (MQ=255)                                                                       gCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGcc                                                                         >  2:873869/1‑126 (MQ=255)                                                                       gCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGcc                                                                         <  1:873869/126‑1 (MQ=255)                                                                            aaaaGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcg                                                       <  2:879365/139‑1 (MQ=255)                                                                                 ccgccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                <  1:1018840/61‑1 (MQ=255)                                                                                 ccgccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                >  2:1018840/1‑61 (MQ=255)                                                                                                   tAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTttggt                                     <  2:1483337/134‑1 (MQ=255)                                                                                                   tAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTttggt                                     >  1:1483337/1‑134 (MQ=255)                                                                                                         aGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTActct                          >  1:1605452/1‑139 (MQ=255)                                                                                                          ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                <  1:509259/36‑1 (MQ=255)                                                                                                          ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                                >  2:509259/1‑36 (MQ=255)                                                                                                          ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTttgg                                      >  2:1630209/1‑126 (MQ=255)                                                                                                          ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTttgg                                      <  1:1630209/126‑1 (MQ=255)                                                                                                          ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACtctc                         >  2:660691/1‑139 (MQ=255)                                                                                                             cTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcgc                                                      <  1:596677/107‑1 (MQ=255)                                                                                                             cTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATCATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcgc                                                      >  2:596677/1‑107 (MQ=255)                                                                                                                aaTGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcg                                                       <  2:986911/103‑1 (MQ=255)                                                                                                                aaTGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcg                                                       >  1:986911/1‑103 (MQ=255)                                                                                                                  tGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTgggg                                <  2:566496/124‑1 (MQ=255)                                                                                                                  tGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTgggg                                >  1:566496/1‑124 (MQ=255)                                                                                                                                 ctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAAt                                                              <  1:1427085/79‑1 (MQ=255)                                                                                                                                 ctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAAt                                                              >  2:1427085/1‑79 (MQ=255)                                                                                                                                 ctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCGGt  >  2:1529001/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                      TAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCGGT  >  NC_000913/1709030‑1709296

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