Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	257,794	T→C	intergenic (+23/‑35)	frsA → / → crl	fermentation‑respiration switch protein/RNA polymerase holoenzyme assembly factor Crl

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	257,794	0	T	C	100.0%	118.3 / NA	31	intergenic (+23/‑35)	frsA/crl	fermentation‑respiration switch protein/RNA polymerase holoenzyme assembly factor Crl
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/15);  total (16/15)

CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257747‑257839                                                |                                              cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                             >  1:1852274/1‑50 (MQ=255)    gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           <  1:713614/49‑1 (MQ=255)    gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           <  1:1731383/49‑1 (MQ=255)    gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           <  1:2036000/49‑1 (MQ=255)    gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                          <  1:235888/50‑1 (MQ=255)        gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAAc                                      <  1:903070/50‑1 (MQ=255)          aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                    >  1:2207896/1‑50 (MQ=255)          aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                    <  1:2249150/50‑1 (MQ=255)          aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                    <  1:2202866/50‑1 (MQ=255)           aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   <  1:2149750/50‑1 (MQ=255)                cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:587881/1‑50 (MQ=255)                 ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:1741521/1‑49 (MQ=255)                   gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            >  1:866233/1‑49 (MQ=255)                   gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            >  1:423908/1‑49 (MQ=255)                    tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcaca                          <  1:2035458/50‑1 (MQ=255)                    tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcaca                          <  1:22499/50‑1 (MQ=255)                    tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcaca                          <  1:1651966/50‑1 (MQ=255)                      ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         <  1:44424/49‑1 (MQ=255)                        aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaaca                       <  1:1844469/49‑1 (MQ=255)                         aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:1392740/1‑49 (MQ=255)                          aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:552319/1‑48 (MQ=255)                             tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:782965/1‑49 (MQ=255)                             tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:1226188/1‑49 (MQ=255)                             tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGAt                 <  1:1139714/50‑1 (MQ=255)                              ttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:1772584/1‑48 (MQ=255)                                 cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc              >  1:549361/1‑49 (MQ=255)                                    aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt           <  1:2069393/49‑1 (MQ=255)                                      ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc        >  1:2186461/1‑50 (MQ=255)                                          gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa    >  1:303383/1‑50 (MQ=255)                                           ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   >  1:690523/1‑50 (MQ=255)                                             aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc  >  1:753565/1‑49 (MQ=255)                                                |                                              CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257747‑257839

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