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%	229.2 / NA	57	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 (35/22);  total (35/22)

GATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCG  >  NC_000913/257750‑257840                                             |                                               gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                            >  1:2108095/1‑49 (MQ=255) gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                            >  1:2064353/1‑49 (MQ=255) gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                            <  1:2310380/49‑1 (MQ=255) gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                            <  1:1436476/49‑1 (MQ=255)  aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                           >  1:2714682/1‑49 (MQ=255)   tCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaa                                          >  1:1408362/1‑49 (MQ=255)    cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:1507901/1‑50 (MQ=255)    cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:1473713/1‑50 (MQ=255)    cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:981855/1‑50 (MQ=255)    cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:2580411/1‑50 (MQ=255)     gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:2497556/1‑49 (MQ=255)     gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:2446047/1‑49 (MQ=255)     gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                        >  1:1657436/1‑49 (MQ=255)     gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAAc                                       <  1:1527740/50‑1 (MQ=255)     gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAAc                                       <  1:1963706/50‑1 (MQ=255)         aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                   <  1:1367506/50‑1 (MQ=255)          aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                  >  1:2611954/1‑50 (MQ=255)          aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                  <  1:2532905/50‑1 (MQ=255)          aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                  >  1:2237460/1‑50 (MQ=255)           cGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTg                                  >  1:2577542/1‑49 (MQ=255)            gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                                <  1:337994/50‑1 (MQ=255)            gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCa                                >  1:2302372/1‑50 (MQ=255)             cTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                               >  1:1071408/1‑50 (MQ=255)              ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                               >  1:774602/1‑49 (MQ=255)                gtgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAATCAGTTGCATca                             <  1:2377423/49‑1 (MQ=39)                 tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATcac                            >  1:2178403/1‑49 (MQ=255)                  gtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:348992/1‑50 (MQ=255)                   ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaac                         <  1:975968/50‑1 (MQ=255)                   ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:923472/1‑49 (MQ=255)                    tAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                          >  1:1507471/1‑48 (MQ=255)                     aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaaca                        <  1:752497/49‑1 (MQ=255)                     aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaaca                        <  1:2757665/49‑1 (MQ=255)                     aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaaca                        <  1:2797050/49‑1 (MQ=255)                     aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:2336617/1‑50 (MQ=255)                     aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:2684754/1‑50 (MQ=255)                      aaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                       >  1:1176783/1‑49 (MQ=255)                          tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                   >  1:1371052/1‑49 (MQ=255)                          tttGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                   >  1:1409984/1‑49 (MQ=255)                             gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg                >  1:2496989/1‑49 (MQ=255)                             gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg                <  1:2483579/49‑1 (MQ=255)                                 aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            <  1:1600990/49‑1 (MQ=255)                                 aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt            >  1:1689761/1‑49 (MQ=255)                                 aaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATg           <  1:753145/50‑1 (MQ=255)                                   ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc         >  1:306076/1‑50 (MQ=255)                                   ttttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc         >  1:449302/1‑50 (MQ=255)                                    tttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg        >  1:2764403/1‑50 (MQ=255)                                    tttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg        >  1:2184504/1‑50 (MQ=255)                                    tttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACg        >  1:898545/1‑50 (MQ=255)                                      tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt       <  1:654806/49‑1 (MQ=255)                                      tGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt       >  1:1634451/1‑49 (MQ=255)                                       gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGtt      <  1:535733/49‑1 (MQ=255)                                       gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGtt      <  1:2604473/49‑1 (MQ=255)                                        ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc    <  1:803473/50‑1 (MQ=255)                                        ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc    <  1:1888485/50‑1 (MQ=255)                                        ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc    <  1:2303721/50‑1 (MQ=255)                                          aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc   >  1:875958/1‑49 (MQ=255)                                          aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCg  <  1:2661864/50‑1 (MQ=255)                                             |                                               GATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCG  >  NC_000913/257750‑257840

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