Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,575,730	T→C	D351G (GAT→GGT)	ydeN ←	sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,575,730	0	T	C	100.0%	114.0 / NA	33	D351G (GAT→GGT)	ydeN	sulfatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/20);  total (13/20)

CGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGC  >  NZ_CP009273/1575599‑1575860                                                                                                                                    |                                                                                                                                   cGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAAc                                                                                                                             >  2:746587/1‑139 (MQ=255)     aaTCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCa                                                                                                                         <  2:819349/139‑1 (MQ=255)              gTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTAt                                                                                                                  <  1:206529/137‑1 (MQ=255)                    aTTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTa                                                                                                          <  2:736834/139‑1 (MQ=255)                         cGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTaaa                                                                                                        >  1:812054/1‑136 (MQ=255)                         cGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAagag                                                                                                     <  1:746587/139‑1 (MQ=255)                                 aGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattg                                                                                             <  1:475631/139‑1 (MQ=255)                                 aGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattg                                                                                             <  2:480917/139‑1 (MQ=255)                                 aGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattg                                                                                             <  2:529544/139‑1 (MQ=255)                                   ttttCCTTTCCACCACTTAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCACGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattgta                                                                                           <  2:178309/139‑1 (MQ=255)                                                       aCATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCa                                                                                                                      >  2:346093/1‑92 (MQ=255)                                                       aCATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCa                                                                                                                      <  1:346093/92‑1 (MQ=255)                                                             gggTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAg                                                                 >  2:670172/1‑139 (MQ=255)                                                                  tgAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAgttgtt                                                            >  1:191876/1‑139 (MQ=255)                                                                        ccGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCgaga                                                       <  2:539336/138‑1 (MQ=255)                                                                         cGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAAt                                                     <  1:798378/139‑1 (MQ=255)                                                                          gCCGGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTACGAGAATg                                                    <  2:62630/139‑1 (MQ=255)                                                                                 tAGGTCTGACTCTTATAGCCGTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTa                                             <  1:670172/139‑1 (MQ=255)                                                                                    gTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattgta                                                                                           <  2:506277/90‑1 (MQ=255)                                                                                    gTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattgta                                                                                           >  1:506277/1‑90 (MQ=255)                                                                                    gTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGc                                                   >  1:45934/1‑130 (MQ=255)                                                                                    gTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGc                                                   <  2:45934/130‑1 (MQ=255)                                                                                         aCTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAg                                                                                                             <  1:138558/67‑1 (MQ=255)                                                                                         aCTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAg                                                                                                             >  2:138558/1‑67 (MQ=255)                                                                                         aCTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattgta                                                                                           >  2:256365/1‑85 (MQ=255)                                                                                         aCTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattgta                                                                                           <  1:256365/85‑1 (MQ=255)                                                                                         aCTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGa                                     >  1:105583/1‑139 (MQ=255)                                                                                               ataGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCTCTAACCGCACCATTATCGGCTGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCaa                                 <  2:118124/137‑1 (MQ=255)                                                                                               ataGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCaa                                 <  1:798594/137‑1 (MQ=255)                                                                                                              cccGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCg                <  2:105583/139‑1 (MQ=255)                                                                                                                          cAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATct    >  2:13985/1‑139 (MQ=255)                                                                                                                            gagGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATctgc  >  2:521238/1‑139 (MQ=255)                                                                                                                            gagGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATctgc  >  1:688501/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                   CGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGC  >  NZ_CP009273/1575599‑1575860

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