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%	148.7 / NA	38	D319E (GAT→GAA)	ydeN	sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (27/11);  total (27/11)

CACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCA  >  NZ_CP009273/1575740‑1575906                                                                                      |                                                                                  cACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCaa                                                                               <  2:202502/90‑1 (MQ=255)        aTCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATa                                                                        >  1:836580/1‑90 (MQ=255)             aGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACgg                                                                   >  1:918843/1‑90 (MQ=255)                              gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                                  >  2:361246/1‑90 (MQ=255)                              gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                                  >  1:86030/1‑90 (MQ=255)                              gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                                  >  2:119632/1‑90 (MQ=255)                              gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                                  <  1:523603/90‑1 (MQ=255)                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctg                                                 >  2:857185/1‑86 (MQ=255)                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctg                                                 <  1:857185/86‑1 (MQ=255)                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGt                                             >  1:414931/1‑90 (MQ=255)                                   gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGt                                             >  2:914880/1‑90 (MQ=255)                                     cATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGttt                                           >  1:429521/1‑90 (MQ=255)                                          tGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAAt                                                                         <  2:663744/55‑1 (MQ=255)                                          tGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAAt                                                                         >  1:663744/1‑55 (MQ=255)                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  2:839762/1‑90 (MQ=255)                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  1:126474/1‑90 (MQ=255)                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  2:330244/1‑90 (MQ=255)                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  2:908144/1‑90 (MQ=255)                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  1:670271/1‑90 (MQ=255)                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  1:641687/1‑90 (MQ=255)                                           gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  1:545860/1‑90 (MQ=255)                                            tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                    <  2:545860/90‑1 (MQ=255)                                            tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                    <  2:86030/90‑1 (MQ=255)                                             ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                   >  2:44444/1‑90 (MQ=255)                                               gTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgct                                               >  1:732913/1‑76 (MQ=255)                                               gTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgct                                               <  2:732913/76‑1 (MQ=255)                                                tttttCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTa                                >  1:554943/1‑90 (MQ=255)                                                          gttgttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCt                      >  2:471504/1‑90 (MQ=255)                                                            tgttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTc                    >  2:207811/1‑90 (MQ=255)                                                             gttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCt                   <  2:918843/90‑1 (MQ=255)                                                             gttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCt                   >  1:423511/1‑90 (MQ=255)                                                               tCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGa                 >  1:300000/1‑90 (MQ=255)                                                                       gCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATc         >  1:548662/1‑90 (MQ=255)                                                                       gCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATc         <  2:554943/90‑1 (MQ=255)                                                                        cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc        <  1:119632/90‑1 (MQ=255)                                                                            ttACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTg    >  2:677698/1‑90 (MQ=255)                                                                             tACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGc   <  2:481973/90‑1 (MQ=255)                                                                              acacCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCa  >  2:78418/1‑90 (MQ=255)                                                                                      |                                                                                  CACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCA  >  NZ_CP009273/1575740‑1575906

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