Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,160,832	G→T	intergenic (‑42/‑235)	sthA ← / → fabR	soluble pyridine nucleotide transhydrogenase/DNA‑binding transcriptional repressor FabR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,160,832	0	G	T	100.0%	210.8 / NA	60	intergenic (‑42/‑235)	sthA/fabR	soluble pyridine nucleotide transhydrogenase/DNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (35/25);  total (35/25)

CTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCT  >  NC_000913/4160710‑4160959                                                                                                                           |                                                                                                                                cTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCggg                                                                                                                <  1:49035/140‑1 (MQ=255) cTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCggg                                                                                                                <  1:15191/140‑1 (MQ=255) cTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCggg                                                                                                                <  1:192054/140‑1 (MQ=255) cTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCggg                                                                                                                <  1:635081/140‑1 (MQ=255)  ttAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGt                                                                                                               <  2:493983/140‑1 (MQ=255)    aaCCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTaa                                                                                                             <  1:449796/140‑1 (MQ=255)    aaCCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTaa                                                                                                             <  2:181174/140‑1 (MQ=255)     aCCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAg                                                                                                            >  2:385211/1‑140 (MQ=255)     aCCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAg                                                                                                            <  1:363132/140‑1 (MQ=255)      ccAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGc                                                                                                           <  1:119479/140‑1 (MQ=255)       cAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGTTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGcc                                                                                                          >  2:312169/1‑140 (MQ=255)       cAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGcc                                                                                                          >  1:180772/1‑140 (MQ=255)         ggCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCaa                                                                                                        >  2:576117/1‑140 (MQ=255)          gCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAAt                                                                                                       >  1:685168/1‑140 (MQ=255)           cccATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAAtt                                                                                                      >  2:493093/1‑140 (MQ=255)               ttGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATc                                                                                                  >  1:682980/1‑140 (MQ=255)                tGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATcc                                                                                                 >  2:121330/1‑140 (MQ=255)                tGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATcc                                                                                                 >  1:516747/1‑140 (MQ=255)                 gCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATccc                                                                                                <  2:516747/140‑1 (MQ=255)                 gCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATccc                                                                                                <  2:105616/140‑1 (MQ=255)                    gcgcCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGt                                                                                             >  2:141395/1‑140 (MQ=255)                     cgcCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTc                                                                                            >  2:640731/1‑140 (MQ=255)                     cgcCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTc                                                                                            >  1:566619/1‑140 (MQ=255)                      gcCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCa                                                                                           <  1:386471/140‑1 (MQ=255)                       ccTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCAt                                                                                          <  2:351692/140‑1 (MQ=255)                        cTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCAtt                                                                                         >  2:712868/1‑140 (MQ=255)                          tCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTa                                                                                       >  2:366108/1‑140 (MQ=255)                           cgccgcCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTC‑TATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATc                                                                                     <  1:661087/140‑1 (MQ=255)                            gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATc                                                                                     <  1:521957/140‑1 (MQ=255)                            gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATc                                                                                     <  1:121330/140‑1 (MQ=255)                            gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATc                                                                                     <  2:509280/140‑1 (MQ=255)                            gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATc                                                                                     <  1:510769/140‑1 (MQ=255)                            gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATc                                                                                     <  1:369959/140‑1 (MQ=255)                            gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATc                                                                                     <  1:107732/140‑1 (MQ=255)                                 cgGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGTTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCt                                                                                <  1:312169/140‑1 (MQ=255)                                 cgGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCt                                                                                >  2:183771/1‑140 (MQ=255)                                                      tGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgct                                                           >  1:600664/1‑140 (MQ=255)                                                      tGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgct                                                           >  1:63890/1‑140 (MQ=255)                                                       ggCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgctt                                                          >  2:715282/1‑140 (MQ=255)                                                        gCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgcttt                                                         >  2:573385/1‑140 (MQ=255)                                                        gCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgctct                                                         >  1:140560/1‑138 (MQ=255)                                                              taatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGt                                                   >  2:453926/1‑140 (MQ=255)                                                              taatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGt                                                   >  2:4512/1‑140 (MQ=255)                                                              taatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGt                                                   >  2:659173/1‑140 (MQ=255)                                                              taatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGt                                                   >  2:708787/1‑140 (MQ=255)                                                                            tgGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGcc                                     >  1:88215/1‑140 (MQ=255)                                                                                 tGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGAc                                >  2:633367/1‑140 (MQ=255)                                                                                  ggTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACg                               <  1:715282/140‑1 (MQ=255)                                                                                   gTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGt                              >  1:218254/1‑140 (MQ=255)                                                                                   gTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGt                              >  2:298017/1‑140 (MQ=255)                                                                                      gggCTTACCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGtttt                           >  2:54917/1‑140 (MQ=255)                                                                                            aCCTGTTCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGa                     >  1:237238/1‑140 (MQ=255)                                                                                                 ttCTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAttttt                >  1:290335/1‑140 (MQ=255)                                                                                                   cTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTAt              >  2:546778/1‑140 (MQ=255)                                                                                                   cTTATACATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTAACTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTAt              <  2:685168/140‑1 (MQ=255)                                                                                                     tataCATAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATcg            <  2:88215/140‑1 (MQ=255)                                                                                                          aTAAAAGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGAt       <  2:565842/140‑1 (MQ=255)                                                                                                            aaaaGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTTATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATAc     >  2:408485/1‑140 (MQ=255)                                                                                                              aaGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATActc   >  2:76783/1‑140 (MQ=255)                                                                                                               aGCAACAGAATGTTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATActct  >  1:397778/1‑140 (MQ=255)                                                                                                                           |                                                                                                                                CTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCT  >  NC_000913/4160710‑4160959

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