Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,706,464	A→T	K39N (AAA→AAT)	rsxB →	SoxR [2Fe‑2S] reducing system protein RsxB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,706,464	0	A	T	100.0%	128.9 / NA	34	K39N (AAA→AAT)	rsxB	SoxR [2Fe‑2S] reducing system protein RsxB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (15/19);  total (15/19)

GCTGGGCCTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAAAATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTCGCCCCTACGCGGAAGCC  >  NC_000913/1706386‑1706536                                                                               |                                                                         gctgGGCCTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAatat                                                                        >  1:307685/1‑81 (MQ=255)  ctgGGCCTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAatat                                                                        <  2:307628/80‑1 (MQ=255)     ggCCTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTg                                                                      <  2:307630/79‑1 (MQ=255)     ggCCTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTg                                                                      <  1:307629/79‑1 (MQ=255)      gCCTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACg                                                                   <  1:307631/81‑1 (MQ=255)       ccTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGa                                                                  <  1:307632/81‑1 (MQ=255)        cTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGaa                                                                 <  2:307633/81‑1 (MQ=255)          ggCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTTACGAAAt                                                               >  1:307673/1‑81 (MQ=255)          ggCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAAt                                                               >  1:307674/1‑81 (MQ=255)             gTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCtt                                                            <  2:307634/81‑1 (MQ=255)              tttGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTa                                                           <  2:307635/81‑1 (MQ=255)               ttGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTAc                                                          >  1:307667/1‑81 (MQ=255)                tGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTAcc                                                         >  1:307672/1‑81 (MQ=255)                 ggCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTAcc                                                         <  1:307637/80‑1 (MQ=255)                 ggCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCg                                                        >  2:307687/1‑81 (MQ=255)                 ggCGACATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTAcc                                                         <  2:307636/80‑1 (MQ=255)                  gcgcCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCg                                                        >  1:307678/1‑80 (MQ=255)                    gcCATTCTGGGTTACGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCg                                                        >  2:307677/1‑78 (MQ=255)                     ccATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCaga                                                    <  2:307638/81‑1 (MQ=255)                          cTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAg                                               >  2:307668/1‑81 (MQ=255)                                 aTGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGt                                         >  1:307684/1‑80 (MQ=255)                                 aTGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGt                                         <  1:307639/80‑1 (MQ=255)                                    ccTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGt                                     <  1:307641/81‑1 (MQ=255)                                        ccgccgTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGc                                   <  2:307642/79‑1 (MQ=255)                                          gccgTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGtt                               <  1:307643/81‑1 (MQ=255)                                          gccgTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGt                                >  2:307691/1‑80 (MQ=255)                                             gTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATc                            >  2:307682/1‑81 (MQ=255)                                                ttGCGGTGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATccc                          <  2:307644/80‑1 (MQ=255)                                                      tGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTc                   >  1:307681/1‑81 (MQ=255)                                                      tGGAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTc                   <  1:307645/81‑1 (MQ=255)                                                        gAAGACGATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTCGc                 <  1:742414/81‑1 (MQ=255)                                                              gATCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTCGCCCCTAc           >  2:307679/1‑81 (MQ=255)                                                                tCCGGTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTCGCCCCTAcgc         <  1:307646/81‑1 (MQ=255)                                                                    gTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTCGCCCCTACGCgg       >  2:307693/1‑79 (MQ=255)                                                                    gTCGTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTCGCCCCTACGCgg       >  1:307694/1‑79 (MQ=255)                                                                       gTTGAGAATATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTCGCCCCTACGCGGAAGcc  <  1:307647/81‑1 (MQ=255)                                                                               |                                                                         GCTGGGCCTGGCGTTTGGCGCCATTCTGGGTTATGCCTCCCGCCGTTTTGCGGTGGAAGACGATCCGGTCGTTGAGAAAATTGACGAAATCTTACCGCAGAGCCAGTGTGGTCAGTGCGGTTATCCCGGCTGTCGCCCCTACGCGGAAGCC  >  NC_000913/1706386‑1706536

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