Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,277,878	A→G	D137G (GAC→GGC)	soxR →	redox‑sensitive transcriptional activator of soxS; autorepressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,277,878	0	A	G	100.0%	159.5 / NA	42	D137G (GAC→GGC)	soxR	redox‑sensitive transcriptional activator of soxS; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/16);  total (26/16)

ACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGC  >  NC_000913/4277818‑4277938                                                             |                                                             aCGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTa                                                       >  1:591298/1‑68 (MQ=255)   ggATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAgg                                                     <  1:782641/68‑1 (MQ=255)    gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGga                                                    >  1:90803/1‑68 (MQ=255)      tGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGgaga                                                  >  2:333193/1‑68 (MQ=255)        tATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGagaag                                                <  1:924648/68‑1 (MQ=255)        tATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCGTAGGagaag                                                <  2:775422/68‑1 (MQ=255)        tATTGGGTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGagaag                                                <  1:325621/68‑1 (MQ=255)         aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGagaaga                                               <  2:615795/68‑1 (MQ=255)         aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGagaaga                                               >  2:549984/1‑68 (MQ=255)           tGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGaagaag                                             >  1:1208/1‑68 (MQ=255)           tGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGaagaag                                             >  2:575385/1‑68 (MQ=255)           tGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGaagaag                                             >  1:769087/1‑68 (MQ=255)           tGGTTGTGGCTGACTTTCGCGCAGTGATTGCCCGTTGCGTGACCCGGGCGGCCGCTTAGGAGaagaag                                             >  1:495599/1‑68 (MQ=255)             gTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGt                                           >  1:827424/1‑68 (MQ=255)             gTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGt                                           <  1:417142/68‑1 (MQ=255)              ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTa                                          >  2:576552/1‑68 (MQ=255)                gtgGCTGCCTTTCGCGCAGTGATTGCCCGTCGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTAcc                                        <  1:695049/68‑1 (MQ=255)                   gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGgaga                                                  >  1:466064/1‑55 (MQ=255)                   gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGgaga                                                  <  2:466066/55‑1 (MQ=255)                   gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGgc                                     <  1:113352/68‑1 (MQ=255)                   gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGgc                                     >  2:282710/1‑68 (MQ=255)                   gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGgc                                     >  1:76261/1‑68 (MQ=255)                   gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGgc                                     >  2:45144/1‑68 (MQ=255)                       ccTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCAc                                 >  1:456768/1‑68 (MQ=255)                          ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCt                              <  1:576550/68‑1 (MQ=255)                          ttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCt                              <  2:547240/68‑1 (MQ=255)                                       tGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGa                 <  2:395036/68‑1 (MQ=255)                                        gCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGaa                <  1:616347/68‑1 (MQ=255)                                         cccGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAc               >  2:749699/1‑68 (MQ=255)                                         cccGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAc               >  2:886565/1‑68 (MQ=255)                                         cccGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAc               >  2:770854/1‑68 (MQ=255)                                         cccGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAc               >  1:482214/1‑68 (MQ=255)                                         cccGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAc               >  1:619391/1‑68 (MQ=255)                                         cccGATGCGTAACACGGGCGGCCGCATAGGAGAAGAAGGTACCGGCGCAAGAATGCTGGAAGATGAAc               >  1:411507/1‑68 (MQ=255)                                          ccGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACa              >  2:315778/1‑68 (MQ=255)                                          ccGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACa              <  1:364448/68‑1 (MQ=255)                                          ccGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACa              <  1:710056/68‑1 (MQ=255)                                           cGTTGCGTACAAAGGGAGGCCGCATAGGAGAAGAAGGTACCGGAGCACGCTTGCTGGAAGATGAACaa             >  2:580079/1‑68 (MQ=255)                                           cGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACaa             >  2:451933/1‑68 (MQ=255)                                            gTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACaaa            <  1:374505/68‑1 (MQ=255)                                                      ccGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGCTGAACAAAACTAAAgcgc  >  1:144715/1‑68 (MQ=255)                                                      ccGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAgcgc  >  2:457405/1‑68 (MQ=255)                                                             |                                                             ACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGC  >  NC_000913/4277818‑4277938

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