Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,707,905	+C	coding (987/2223 nt)	rsxC →	SoxR [2Fe‑2S] reducing system protein RsxC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,707,901	1	.	C	100.0%	90.7 / NA	26	C328S (TGC→TCC)	rsxC	SoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (12/14);  total (12/14)

GAAACGTGCCGTTATTGATGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTG‑CCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAACTGTCTGTTGGCTCCCTCTGCCAATGAACTTGGCGAACC  >  NC_000913/1707764‑1708030                                                                                                                                           |                                                                                                                                  gAAACGTGCCGTTATTGATGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGTCGCTAGTGCGTCATTTATTGAATGATGCCGGATTCTG‑cc                                                                                                                                 >  2:218060/1‑140 (MQ=255) gAAACGTGCCGTTATTGATGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTG‑cc                                                                                                                                 >  1:37566/1‑140 (MQ=255) gAAACGTGCCGTTATTGATGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTG‑cc                                                                                                                                 >  1:38605/1‑140 (MQ=255) gAAACGTGCCGTTATTGATGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCAAGGCTGGGGACGCTAGTGCGTCATTTATTGAATGATGCCGGATTCTG‑cc                                                                                                                                 >  2:240648/1‑140 (MQ=255)  aaaCGTGCCGTTATTGATGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTG‑ccc                                                                                                                                >  1:9514/1‑140 (MQ=255)              aTTGATGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATc                                                                                                                     >  1:12666/1‑140 (MQ=255)                   tGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATg                                                                                                                <  2:38605/140‑1 (MQ=255)                   tGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATg                                                                                                                <  2:367945/140‑1 (MQ=255)                    ggCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATgg                                                                                                               >  1:336234/1‑140 (MQ=255)                        aGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGAt                                                                                                           <  1:335086/140‑1 (MQ=255)                         gCCGATTTCCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGAtt                                                                                                          >  2:387601/1‑140 (MQ=255)                                   gAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGcc                                                                                                >  2:500602/1‑140 (MQ=255)                                   gAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGcc                                                                                                <  2:336234/140‑1 (MQ=255)                                         gttgtAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTaa                                                                                          <  1:387601/140‑1 (MQ=255)                                          ttgtAACCCTGACTGCCGAAGCCATAGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGACTATGGGTGGCCCGCTAAt                                                                                         <  1:177415/140‑1 (MQ=255)                                                   tGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTAc                                                                                >  1:570339/1‑140 (MQ=255)                                                             aGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATgg                                                                      <  2:570339/140‑1 (MQ=255)                                                                   cgctcgcCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGAt                                                                >  1:546973/1‑140 (MQ=255)                                                                       cgcCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTg                                                                   >  1:605545/1‑133 (MQ=255)                                                                       cgcCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTg                                                                   <  2:605545/133‑1 (MQ=255)                                                                            gggCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGTTCAAATGGTGATTATTGGTGGCCCGCTATTGGCCTTTACCTCGCCATGGCTGGATGTCCCGGTc                                                       >  2:293787/1‑140 (MQ=255)                                                                               cAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTa                                                    <  2:9514/140‑1 (MQ=255)                                                                                aaCGTCTGGGCACGGGTGGGGACGCGAGTCCGTCATTGATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTaa                                                   <  1:293787/140‑1 (MQ=255)                                                                                aaCGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGAGTATGGGTGGCCCGCTAATGGGCTTTCCCTTGCCATGGCTGGATGTCCCGGTCGTaa                                                   <  2:37566/140‑1 (MQ=255)                                                                                       gggCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCTTGACTGAATGTCCCGGTCGTAAAGATTAc                                            >  1:264856/1‑140 (MQ=255)                                                                                        ggCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGGTTCTGCCCCCTCTGCCGCTCAACTGGTGATTATGGGTGGCCCGCTTATGGGCTTTACCTTACCGTGGCTGGCTGTACCGGTCGGAACGATTAcc                                           >  2:509067/1‑140 (MQ=255)                                                                                                 ggggACGCCAGTGCGCCATTTAGTGAATGTTGCCGGAGTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAActgtctg                                  <  1:218060/140‑1 (MQ=255)                                                                                                  gggACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAActgtctgt                                 >  1:44070/1‑140 (MQ=255)                                                                                                     aCGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAACTGTCTGTTgg                              <  1:500602/140‑1 (MQ=255)                                                                                                                          aaTGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAACTGTCTGTTGGCTcc                          <  1:96712/123‑1 (MQ=255)                                                                                                                          aaTGATGCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAACTGTCTGTTGGCTcc                          >  2:96712/1‑123 (MQ=255)                                                                                                                               tgCCGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAACTGTCTGTTGGCTCCCTCTGCCAATGAACTTGGCGaa    >  1:666606/1‑140 (MQ=255)                                                                                                                                 ccGGATTCTGCCCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAACTGTCTGTTGGCTCCCTCTGCCAATGAACTTGGCGAAcc  <  2:546973/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                  GAAACGTGCCGTTATTGATGGCGAGCCGATTACCGAGCGTGTTGTAACCCTGACTGGCGAAGCAATCGCTCGCCCGGGCAACGTCTGGGCACGGCTGGGGACGCCAGTGCGTCATTTATTGAATGATGCCGGATTCTG‑CCCCTCTGCCGATCAAATGGTGATTATGGGTGGCCCGCTAATGGGCTTTACCTTGCCATGGCTGGATGTCCCGGTCGTAAAGATTACCAACTGTCTGTTGGCTCCCTCTGCCAATGAACTTGGCGAACC  >  NC_000913/1707764‑1708030

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