Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,575,825	A→T	D319E (GAT→GAA)	ydeN ←	sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,575,825	0	A	T	100.0%	150.4 / NA	39	D319E (GAT→GAA)	ydeN	sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (22/17);  total (22/17)

ACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGA  >  NZ_CP009273/1575741‑1575909                                                                                     |                                                                                     aCCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAc                                                                                 <  2:1078737/90‑1 (MQ=255) aCCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAc                                                                                 <  1:183974/90‑1 (MQ=255)    aTTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGa                                                                              >  2:416103/1‑90 (MQ=255)      tATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAAt                                                                            <  1:384145/90‑1 (MQ=255)            aGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACgg                                                                      >  1:1174101/1‑90 (MQ=255)             ggTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                     <  1:1010672/90‑1 (MQ=255)                aaaGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGc                                                                  <  2:968368/90‑1 (MQ=255)                 aaGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCg                                                                 >  2:1081457/1‑90 (MQ=255)                  agagAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgt                                                                <  1:1174537/90‑1 (MQ=255)                  agagAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgt                                                                <  2:56037/90‑1 (MQ=255)                  agagAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgt                                                                <  1:254450/90‑1 (MQ=255)                     gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                             >  1:1143772/1‑90 (MQ=255)                       ataatTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGtagtag                                                           <  1:1081457/90‑1 (MQ=255)                               attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                   >  2:95152/1‑90 (MQ=255)                               attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                   >  1:38020/1‑90 (MQ=255)                                ttgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgct                                                  >  2:589266/1‑90 (MQ=255)                                    cATACTTTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTAGACCCTGTTCAACAGTAGAAACGGAATCGTAGTAGTTATCTGCTGttt                                              >  1:774622/1‑90 (MQ=255)                                    cATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGttt                                              >  2:710985/1‑90 (MQ=255)                                     aTACTGTCCGTTTTTCTTCAGTTGTTCGTGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTg                                             >  2:222989/1‑90 (MQ=255)                                          gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAc                                                                                 <  2:192217/49‑1 (MQ=255)                                          gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAAc                                                                                 >  1:192217/1‑49 (MQ=255)                                          gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                        >  1:5448/1‑90 (MQ=255)                                           tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                       <  1:376639/90‑1 (MQ=255)                                            ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                      >  1:784746/1‑90 (MQ=255)                                            ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                      >  1:358692/1‑90 (MQ=255)                                            ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                      >  1:331795/1‑90 (MQ=255)                                            ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                      >  1:314598/1‑90 (MQ=255)                                             cGTTTTTCTTCAGTTGTTTGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCgg                                     <  2:817636/90‑1 (MQ=255)                                                           tgttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTc                       >  1:1101852/1‑90 (MQ=255)                                                             ttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTg                     <  2:583657/90‑1 (MQ=255)                                                              tCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGa                    >  1:1134223/1‑90 (MQ=255)                                                                       cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc           <  1:1035937/90‑1 (MQ=255)                                                                       cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc           <  2:1048070/90‑1 (MQ=255)                                                                       cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc           <  1:589266/90‑1 (MQ=255)                                                                           ttACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTg       >  2:1163341/1‑90 (MQ=255)                                                                           ttACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTg       >  2:1155998/1‑90 (MQ=255)                                                                            tACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGc      <  2:38020/90‑1 (MQ=255)                                                                                cccTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGa  >  2:1029362/1‑90 (MQ=255)                                                                                cccTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGa  >  2:957497/1‑90 (MQ=255)                                                                                     |                                                                                     ACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCAGGA  >  NZ_CP009273/1575741‑1575909

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