Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,576,388	A→G	S132P (TCC→CCC)	ydeN ←	sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,576,388	0	A	G	100.0%	171.3 / NA	46	S132P (TCC→CCC)	ydeN	sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/22);  total (24/22)

GAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCA  >  NZ_CP009273/1576301‑1576470                                                                                        |                                                                                   gAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGgaa                                                                                  >  1:855044/1‑90 (MQ=255)  aaTACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAac                                                                                 <  1:1097900/90‑1 (MQ=255)    tACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                               <  1:2472329/90‑1 (MQ=255)    tACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                               <  2:2089348/90‑1 (MQ=255)        aTCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtg                                                                           >  2:2029113/1‑90 (MQ=255)         tCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtgt                                                                          <  2:84555/90‑1 (MQ=255)         tCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtgt                                                                          <  1:863212/90‑1 (MQ=255)         tCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtgt                                                                          <  2:652397/90‑1 (MQ=255)           cTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                        >  2:2966799/1‑90 (MQ=255)           cTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                        <  2:1853587/90‑1 (MQ=255)           cTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                        <  1:2472520/90‑1 (MQ=255)             gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                               <  2:1422489/81‑1 (MQ=255)             gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                               >  1:1422489/1‑81 (MQ=255)             gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                               <  2:135028/81‑1 (MQ=255)             gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                               >  1:135028/1‑81 (MQ=255)             gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCca                                                                      >  2:1841093/1‑90 (MQ=255)             gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCca                                                                      >  2:1203191/1‑90 (MQ=255)             gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCca                                                                      >  1:2801321/1‑90 (MQ=255)               gCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCcaca                                                                    >  2:1386056/1‑90 (MQ=255)                cATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACAt                                                                   <  2:1322787/90‑1 (MQ=255)                  tCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAg                                                                 >  1:1133905/1‑90 (MQ=255)                  tCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAg                                                                 >  2:1593221/1‑90 (MQ=255)                      tATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGt                                                             <  2:2800705/90‑1 (MQ=255)                            aaTAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTaa                                                       >  2:2745143/1‑90 (MQ=255)                               aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTaaaa                                                     >  2:1625515/1‑89 (MQ=255)                               aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTaaaa                                                     <  1:1625515/89‑1 (MQ=255)                               aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                    >  2:431750/1‑90 (MQ=255)                               aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                    >  1:2752284/1‑90 (MQ=255)                               aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                    >  1:2463598/1‑90 (MQ=255)                               aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                    >  2:2188899/1‑90 (MQ=255)                               aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                    >  2:2146753/1‑90 (MQ=255)                               aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                    >  1:757854/1‑90 (MQ=255)                                   acCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGGCACATAGCCGTTAGTAAAACGTAc                                                <  2:2159050/90‑1 (MQ=255)                                   acCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTAc                                                <  2:1874779/90‑1 (MQ=255)                                   acCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTAc                                                <  1:1841093/90‑1 (MQ=255)                                    ccAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACg                                               >  1:2761615/1‑90 (MQ=255)                                     cAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGc                                              <  2:2761615/90‑1 (MQ=255)                                         gcgggcggGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAac                                                                                 <  2:2572903/51‑1 (MQ=255)                                         gcgggcggGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAac                                                                                 >  1:2572903/1‑51 (MQ=255)                                           gggcgggAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCAt                                        <  1:1203191/90‑1 (MQ=255)                                                      cGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACAt                                                                   >  1:1590059/1‑52 (MQ=255)                                                      cGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACAt                                                                   <  2:1590059/52‑1 (MQ=255)                                                                      cggcgCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTg             >  1:2775802/1‑90 (MQ=255)                                                                        gcgcGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAt           >  1:2282509/1‑90 (MQ=255)                                                                          gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTAc                                                >  2:1457337/1‑51 (MQ=255)                                                                          gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTAc                                                <  1:1457337/51‑1 (MQ=255)                                                                                 gggCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCa  <  1:2751496/90‑1 (MQ=255)                                                                                        |                                                                                   GAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCA  >  NZ_CP009273/1576301‑1576470

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