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

GCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCC  >  NZ_CP009273/1575739‑1575900                                                                                       |                                                                            gCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCa                                                                          >  1:749358/1‑90 (MQ=255)       ttATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGaa                                                                    <  2:410701/90‑1 (MQ=255)         aTCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATa                                                                  >  1:312729/1‑90 (MQ=255)         aTCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATa                                                                  >  1:700673/1‑90 (MQ=255)         aTCGGAGGGAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATa                                                                  >  1:513219/1‑90 (MQ=255)              aGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACgg                                                             >  2:449837/1‑90 (MQ=255)               ggTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                            <  2:537578/90‑1 (MQ=255)                   aaGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCg                                                        >  1:95826/1‑90 (MQ=255)                       gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACCCCCTGTTCAACAGAATAAACGGAAGCgtagt                                                    <  1:194761/90‑1 (MQ=255)                       gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATaa                                                                 <  1:510276/77‑1 (MQ=255)                       gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATaa                                                                 >  2:510276/1‑77 (MQ=255)                       gaATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgtagt                                                    >  1:390354/1‑90 (MQ=255)                         ataatTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGtagtag                                                  >  1:133601/1‑90 (MQ=255)                               gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                            <  2:325694/90‑1 (MQ=255)                               gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                            >  1:72468/1‑90 (MQ=255)                                  ttgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgct                                         <  2:232765/90‑1 (MQ=255)                                    gtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGt                                       >  2:99783/1‑90 (MQ=255)                                           tGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTa                                <  1:248234/90‑1 (MQ=255)                                           tGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTa                                >  1:57023/1‑90 (MQ=255)                                           tGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTa                                >  2:521445/1‑90 (MQ=255)                                            gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                               >  1:2420/1‑90 (MQ=255)                                             tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                              <  2:247034/90‑1 (MQ=255)                                              ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                             >  1:45141/1‑90 (MQ=255)                                              ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                             >  1:67431/1‑90 (MQ=255)                                              ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                             <  1:238843/90‑1 (MQ=255)                                              ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                             >  1:433357/1‑90 (MQ=255)                                                 tttttCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTa                          >  1:684078/1‑90 (MQ=255)                                                               ttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTg            <  1:353213/90‑1 (MQ=255)                                                                        gCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATc   <  2:467304/90‑1 (MQ=255)                                                                         cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc  <  2:45141/90‑1 (MQ=255)                                                                                       |                                                                            GCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCC  >  NZ_CP009273/1575739‑1575900

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: ‑