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%	271.6 / NA	65	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 (38/27);  total (38/27)

CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCG  >  NC_000913/257747‑257840                                                |                                               cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                              >  1:2748557/1‑50 (MQ=255) cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                              >  1:3095273/1‑50 (MQ=255) cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                              <  1:2832577/50‑1 (MQ=255)  tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                             >  1:2217142/1‑50 (MQ=255)   ggATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTgg                                             >  1:3021729/1‑49 (MQ=255)    gATCGAAAAACGCTTGTGTTAAAAATTTGCTAACTTTTGCCAATCTggc                                            >  1:1069405/1‑48 (MQ=38)    gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                            >  1:2305507/1‑49 (MQ=255)    gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                           <  1:4276088/50‑1 (MQ=255)     aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                           >  1:2670536/1‑49 (MQ=255)       cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:1246734/1‑50 (MQ=255)        gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:1087100/1‑49 (MQ=255)        gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:388749/1‑49 (MQ=255)        gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:1858245/1‑49 (MQ=255)        gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAAc                                       <  1:1130197/50‑1 (MQ=255)         aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:2807145/1‑48 (MQ=255)          aaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                     <  1:1998677/50‑1 (MQ=255)           aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                    <  1:2613077/50‑1 (MQ=255)             aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                  <  1:3037735/50‑1 (MQ=255)             aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                  <  1:948612/50‑1 (MQ=255)               gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                                <  1:1710692/50‑1 (MQ=255)                  tgtgTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                             >  1:3972244/1‑50 (MQ=255)                  tgtgTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                             >  1:624298/1‑50 (MQ=255)                   gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                             <  1:1833727/49‑1 (MQ=255)                    tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcac                            >  1:3343968/1‑49 (MQ=255)                     gtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:278497/1‑50 (MQ=255)                      ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaac                         <  1:256519/50‑1 (MQ=255)                      ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:2780136/1‑49 (MQ=255)                      ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:3595489/1‑49 (MQ=255)                       tAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:649449/1‑48 (MQ=255)                        aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:2249423/1‑50 (MQ=255)                         aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAgg                      <  1:2332765/50‑1 (MQ=255)                         aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:3156671/1‑49 (MQ=255)                          aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:3772495/1‑48 (MQ=255)                          aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGga                     <  1:524194/50‑1 (MQ=255)                            aTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                   <  1:2890544/50‑1 (MQ=255)                             tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                   >  1:954592/1‑49 (MQ=255)                             tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                   >  1:2043496/1‑49 (MQ=255)                             tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGAt                  <  1:1295147/50‑1 (MQ=255)                             tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGAt                  <  1:3187061/50‑1 (MQ=255)                              ttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                   >  1:1647546/1‑48 (MQ=255)                              ttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                   >  1:3868399/1‑48 (MQ=255)                               tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATa                 <  1:2226526/49‑1 (MQ=255)                               tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg                >  1:2239575/1‑50 (MQ=255)                                gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg                <  1:4118156/49‑1 (MQ=255)                                 cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGc               >  1:2681425/1‑49 (MQ=255)                                  tAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa             >  1:4270678/1‑50 (MQ=255)                                    aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            <  1:2980617/49‑1 (MQ=255)                                    aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATg           <  1:2586746/50‑1 (MQ=255)                                     aTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGa          <  1:1478870/50‑1 (MQ=255)                                        ttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg        <  1:4088115/49‑1 (MQ=255)                                        ttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt       >  1:460109/1‑50 (MQ=255)                                          gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGtt      <  1:2869728/49‑1 (MQ=255)                                          gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa     >  1:2340831/1‑50 (MQ=255)                                           ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa     <  1:354609/49‑1 (MQ=255)                                            cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  1:4428172/1‑50 (MQ=255)                                            cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  1:1884596/1‑50 (MQ=255)                                            cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   <  1:2787999/50‑1 (MQ=255)                                            cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  1:3276759/1‑50 (MQ=255)                                            cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  1:3489742/1‑50 (MQ=255)                                             aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  1:2699086/1‑49 (MQ=255)                                             aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  1:2803329/1‑49 (MQ=255)                                             aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  1:888287/1‑49 (MQ=255)                                             aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCg  <  1:574448/50‑1 (MQ=255)                                             aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCg  <  1:1558621/50‑1 (MQ=255)                                             aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCg  <  1:1205399/50‑1 (MQ=255)                                                |                                               CTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCG  >  NC_000913/257747‑257840

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: ‑