Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,256	A→G	G114G (GGT→GGC)	yjcS ←	alkyl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,256	0	A	G	100.0%	124.5 / NA	35	G114G (GGT→GGC)	yjcS	alkyl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/12);  total (23/12)

AGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGCAGAATACCTTCATCCAGCAGC  >  NC_000913/4306133‑4306391                                                                                                                            |                                                                                                                                        aGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCTAGAAATGCCGTTGATTTGCGACt                                                                                                                         >  2:8290/1‑140 (MQ=255)                              taatGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTCATTTGCGACTGACGCCCCAGGCTGGGGTTTACGGTTTcc                                                                                            >  2:267840/1‑140 (MQ=255)                                      ttCTCGCCCTCAACGAATGTAACGTTAGAGATATCCTGGCCGCGCACCCGATACATTTTATCGGTGACTTTGAACAGGCAAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGc                                                                                    <  2:175526/140‑1 (MQ=255)                                      ttCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGc                                                                                    >  1:406738/1‑140 (MQ=255)                                        ctcGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCgg                                                                                  >  1:396895/1‑140 (MQ=255)                                                cAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAata                                                                          >  1:562229/1‑140 (MQ=255)                                                cAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAata                                                                          >  2:59873/1‑140 (MQ=255)                                                 aacgaacgTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAatat                                                                         <  2:369893/140‑1 (MQ=255)                                                                 gagaTATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGAGTTTACGGTTTCAGGCCCTGCGGCAATAATTTCAAATTTGTAATCGTc                                                         >  2:280876/1‑140 (MQ=255)                                                                    atatCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCcgc                                                      >  1:498226/1‑140 (MQ=255)                                                                     tatCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCcgcg                                                     <  2:406738/140‑1 (MQ=255)                                                                     tatCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCcgcg                                                     <  2:659162/140‑1 (MQ=255)                                                                      atCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCcgcg                                                     <  1:155501/139‑1 (MQ=255)                                                                      atCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCcgcg                                                     >  2:155501/1‑139 (MQ=255)                                                                        ccTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCgg                                                  >  2:69371/1‑140 (MQ=255)                                                                                cgcACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAAc                                          >  2:546290/1‑140 (MQ=255)                                                                                cgcACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAAc                                          >  1:464450/1‑140 (MQ=255)                                                                                cgcACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAAc                                          >  1:365198/1‑140 (MQ=255)                                                                                    ccTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTa                                      <  2:464450/140‑1 (MQ=255)                                                                                        aTACATTTTATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGt                                  <  1:439257/140‑1 (MQ=255)                                                                                             ttttATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCa                             <  1:59873/140‑1 (MQ=255)                                                                                               ttATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGGGGAATTATTTTAAAATTTGTAATCTTCCGGGCGGTAGTAAACTTTACGGtcggcttc                           >  2:497217/1‑133 (MQ=255)                                                                                               ttATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATc                           >  2:584127/1‑140 (MQ=255)                                                                                               ttATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATc                           >  2:48328/1‑140 (MQ=255)                                                                                                tATCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCa                          <  2:396895/140‑1 (MQ=255)                                                                                                 aTCGGTGACTTTGAATAGGCCATAATTGCGGTTGATTAGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCAc                         <  1:267840/140‑1 (MQ=255)                                                                                                  tCGGTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACg                        <  2:498226/140‑1 (MQ=255)                                                                                                     gTGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGCAg                     >  2:591543/1‑140 (MQ=255)                                                                                                      tGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGAAGTAAACTTTACCTTTCGCATCACGCAGa                    >  1:201374/1‑140 (MQ=255)                                                                                                      tGACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTACGCGCGGTAGTGAAGTTTACCATTCGCATCACGCAGa                    >  1:311823/1‑140 (MQ=255)                                                                                                       gACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGCAGaa                   >  2:109322/1‑140 (MQ=255)                                                                                                       gACTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGCAGaa                   <  1:495450/140‑1 (MQ=255)                                                                                                        aCTTTGAACAGGCCAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGCAGAAt                  >  2:396831/1‑140 (MQ=255)                                                                                                                    ccAGAAATGCCGTTGATTTGCGACTGACGCCACAGGCTGTGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTAACGTTCGCAACACGCAGAATACCTTCATCcag      >  2:287654/1‑140 (MQ=255)                                                                                                                        aaaTGCCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGCAGAATACCTTCATCCagcagc  >  2:13558/1‑140 (MQ=255)                                                                                                                            |                                                                                                                                        AGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGCAGAATACCTTCATCCAGCAGC  >  NC_000913/4306133‑4306391

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