Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,158,586	A→G	S33G (AGC→GGC)	oxyR →	oxidative and nitrosative stress transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,158,586	0	A	G	100.0%	184.3 / NA	50	S33G (AGC→GGC)	oxyR	oxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/23);  total (27/23)

CCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGT  >  NC_000913/4158451‑4158708                                                                                                                                        |                                                                                                                           ccTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCg                                                                                                                         >  2:790552/1‑139 (MQ=255)        ccATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAAt                                                                                                                  >  1:258991/1‑139 (MQ=255)                       gCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGa                                                                                                   >  1:832054/1‑139 (MQ=255)                            ggaggaTGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTg                                                                                                        >  1:836088/1‑129 (MQ=255)                            ggaggaTGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTg                                                                                                        <  2:836088/129‑1 (MQ=255)                            ggaggaTGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGc                                                                                              >  2:80884/1‑139 (MQ=255)                                     atAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGa                                                                                     >  2:603073/1‑139 (MQ=255)                                              atTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTgg                                                                            >  1:26433/1‑139 (MQ=255)                                                     aTCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     <  2:258991/139‑1 (MQ=255)                                                               ccTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     <  2:156025/129‑1 (MQ=255)                                                               ccTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     >  1:156025/1‑129 (MQ=255)                                                                cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCt                                                                                             >  2:664401/1‑104 (MQ=255)                                                                cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCt                                                                                             <  1:664401/104‑1 (MQ=255)                                                                cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGAt                                                                                    >  1:513171/1‑113 (MQ=255)                                                                cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGAt                                                                                    <  2:513171/113‑1 (MQ=255)                                                                cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     <  1:508355/128‑1 (MQ=255)                                                                cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     >  2:508355/1‑128 (MQ=255)                                                                    tggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGc                                                                               <  2:312644/114‑1 (MQ=255)                                                                    tggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGc                                                                               >  1:312644/1‑114 (MQ=255)                                                                         ttGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCa                                                 >  1:441495/1‑139 (MQ=255)                                                                          tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGaa                                                                                                     >  2:449453/1‑86 (MQ=255)                                                                          tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGaa                                                                                                     <  1:449453/86‑1 (MQ=255)                                                                          tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGATCAGCCGTAAAGTGTTGTTCAc                                                >  1:493277/1‑139 (MQ=255)                                                                          tGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCAc                                                >  2:308368/1‑139 (MQ=255)                                                                           ggCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCAcc                                               <  1:794709/139‑1 (MQ=255)                                                                            gCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGt                                                                                       <  1:628811/98‑1 (MQ=255)                                                                            gCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGt                                                                                       >  2:628811/1‑98 (MQ=255)                                                                                aaCACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGc                                          >  1:155038/1‑139 (MQ=255)                                                                                    ccGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGa                                      <  2:637899/139‑1 (MQ=255)                                                                                      gCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAAt                                    >  1:820812/1‑139 (MQ=255)                                                                                      gCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAAt                                    >  1:706046/1‑139 (MQ=255)                                                                                          ttttCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAAtgctg                                >  2:756844/1‑139 (MQ=255)                                                                                              cGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCtgg                            <  1:756844/139‑1 (MQ=255)                                                                                                   tGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGAt                       <  1:790552/139‑1 (MQ=255)                                                                                                   tGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGAt                       <  2:441495/139‑1 (MQ=255)                                                                                                        cAGATTCCTGCCACGTTAGCCAGCCGACGCTTGGCGCGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGgc                  <  1:526237/139‑1 (MQ=255)                                                                                                             tCCTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTa             <  2:820812/139‑1 (MQ=255)                                                                                                              ccTGCCACGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTAc            <  2:310280/139‑1 (MQ=255)                                                                                                                     cGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCt                                                                                             <  2:553096/51‑1 (MQ=255)                                                                                                                     cGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     >  2:519827/1‑75 (MQ=255)                                                                                                                     cGTTAGCCAGCCGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     <  1:519827/75‑1 (MQ=255)                                                                                                                     cGTTAGCCAGCCGACGCTAGGCGGGCAAATTCGTAAGCTGGAAGATGAGCt                                                                                             >  1:553096/1‑51 (MQ=255)                                                                                                                            cagccGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGt  <  1:533887/135‑1 (MQ=255)                                                                                                                            cagccGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGt  >  2:533887/1‑135 (MQ=255)                                                                                                                               ccGACGCTTGGTGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAg                                                                  >  1:438296/1‑68 (MQ=255)                                                                                                                               ccGACGCTTGGTGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAg                                                                  <  2:438296/68‑1 (MQ=255)                                                                                                                                cGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     <  2:187069/64‑1 (MQ=255)                                                                                                                                cGACGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGAc                                                                     >  1:187069/1‑64 (MQ=255)                                                                                                                                  aCGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGtgc       <  2:31447/124‑1 (MQ=255)                                                                                                                                  aCGCTTGGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGtgc       >  1:31447/1‑124 (MQ=255)                                                                                                                                        |                                                                                                                           CCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGCTGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAAGTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTGCGT  >  NC_000913/4158451‑4158708

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