Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,710,659	T→C	G152G (GGT→GGC)	rsxG →	SoxR [2Fe‑2S] reducing system protein RsxG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,710,659	0	T	C	100.0%	95.4 / NA	27	G152G (GGT→GGC)	rsxG	SoxR [2Fe‑2S] reducing system protein RsxG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/15);  total (12/15)

GGTGGGAGCCGATTTTAACGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGTGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCCGCGCGGTGGTTAATGCGGTAAAACGCGCCGGATTGTACGCT  >  NC_000913/1710524‑1710776                                                                                                                                        |                                                                                                                      ggtggGAGCCGATTTTAACGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAg                                                                                                                   <  1:270104/140‑1 (MQ=255) ggtggGAGCCGATTTTAACGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAg                                                                                                                   >  1:554782/1‑140 (MQ=255)      gAGCCGATTTTAACGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAgatgat                                                                                                              <  1:373492/140‑1 (MQ=255)         ccGATTTTAACGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATgcg                                                                                                           >  1:470057/1‑140 (MQ=255)           gATTTTAACGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCa                                                                                                         <  2:554782/140‑1 (MQ=255)            aTTTTAACGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCAc                                                                                                        >  1:753147/1‑140 (MQ=255)                   cGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCgg                                                                                                 <  2:753147/140‑1 (MQ=255)                                ggCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGAtgg                                                                                    >  2:247545/1‑140 (MQ=255)                                   aCGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGAtggtgg                                                                                 >  2:913708/1‑140 (MQ=255)                                   aCGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGAtggtgg                                                                                 >  2:875332/1‑140 (MQ=255)                                           gACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCg                                                                         <  1:913708/140‑1 (MQ=255)                                             cAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGAc                                                                       <  2:90092/140‑1 (MQ=255)                                                     cacGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCAc                                                               >  1:520387/1‑140 (MQ=255)                                                       cGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGAc                                                                       <  1:450626/130‑1 (MQ=255)                                                       cGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGAc                                                                       >  2:450626/1‑130 (MQ=255)                                                          aaCGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGgcg                                                          >  1:883225/1‑140 (MQ=255)                                                           aCGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGgcgc                                                         <  2:23021/140‑1 (MQ=255)                                                                       ggCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTcc                                             >  1:106520/1‑140 (MQ=255)                                                                         cGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTcccc                                           <  2:470057/140‑1 (MQ=255)                                                                            tAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCcgcg                                        <  2:520387/140‑1 (MQ=255)                                                                            tAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCcgcg                                        <  2:890788/140‑1 (MQ=255)                                                                            tAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCcgcg                                        <  1:875332/140‑1 (MQ=255)                                                                              aaaTCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCcgcgcg                                      <  2:488545/140‑1 (MQ=255)                                                                                   gAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCCCCGCggtggt                                 <  2:674791/140‑1 (MQ=255)                                                                                                        aTCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCCGCGCGGTGGTTAATGCGGTAAAACGCGCCgg            >  1:551253/1‑140 (MQ=255)                                                                                                        aTCACCCATTTTGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCAACCAGTTCACCGGCGCGACGATTACTCCCCGCGCGGTGGTTAATGCGGTAAAACGCGCCgg            >  2:118645/1‑140 (MQ=255)                                                                                                                  ttGCGGGTAAAAAAATCAGTGGCGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCCGCGCGGTGGTTAATGCGGTAAAACGCGCCGGATTGTACGCt  <  1:33686/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                      GGTGGGAGCCGATTTTAACGGCACGGTACTTGGCACGCGCGTGACAGAGCACCACGAAACGCCAGGGCTTGGCGATAAAATCGAACTGCGCCTTTCTGACTGGATCACCCATTTTGCGGGTAAAAAAATCAGTGGTGCAGATGATGCGCACTGGGCGGTGAAGAAAGATGGTGGTGATTTCGACCAGTTCACCGGCGCGACGATTACTCCCCGCGCGGTGGTTAATGCGGTAAAACGCGCCGGATTGTACGCT  >  NC_000913/1710524‑1710776

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