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	92.9%	145.7 / 2.2	42	D319E (GAT→GAA)	ydeN	sulfatase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (24/15);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCA  >  NZ_CP009273/1575739‑1575906                                                                                       |                                                                                  gCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCa                                                                                >  2:17901/1‑90 (MQ=255) gCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCa                                                                                >  2:251190/1‑90 (MQ=255) gCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCa                                                                                >  2:347172/1‑90 (MQ=255)      aTTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGa                                                                           >  1:101063/1‑90 (MQ=255)       ttATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGaa                                                                          <  1:727467/90‑1 (MQ=255)              aGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACgg                                                                   >  1:153283/1‑90 (MQ=255)               ggTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                  <  2:73568/90‑1 (MQ=255)               ggTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGa                                                                  <  1:251190/90‑1 (MQ=255)                gTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGaa                                                                 >  1:330577/1‑90 (MQ=255)                 tAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAg                                                                <  2:272422/90‑1 (MQ=255)                    agagAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCgt                                                             <  2:629637/90‑1 (MQ=255)                               gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                                  >  1:737503/1‑90 (MQ=255)                               gtattgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATct                                                  >  1:151863/1‑90 (MQ=255)                                 attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                >  1:398048/1‑90 (MQ=255)                                 attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                >  1:695576/1‑90 (MQ=255)                                 attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                >  1:723610/1‑90 (MQ=255)                                 attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                >  1:97285/1‑90 (MQ=255)                                 attgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATctgc                                                >  2:205675/1‑90 (MQ=255)                                   tgtCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATctgctg                                              >  2:835543/1‑90 (MQ=255)                                      cATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGttt                                           >  2:735560/1‑90 (MQ=255)                                      cATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGttt                                           >  2:347252/1‑90 (MQ=255)                                         aCTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGAc                                        >  1:237335/1‑90 (MQ=255)                                            gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  1:105173/1‑90 (MQ=255)                                            gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  2:1181/1‑90 (MQ=255)                                            gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  1:318145/1‑90 (MQ=255)                                            gTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAc                                     >  1:398871/1‑90 (MQ=255)                                             tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                    <  1:448147/90‑1 (MQ=255)                                             tCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                    <  2:330577/90‑1 (MQ=255)                                              ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                    >  1:162885/1‑89 (MQ=255)                                              ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTAcc                                    <  2:162885/89‑1 (MQ=255)                                              ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                   >  2:154514/1‑90 (MQ=255)                                              ccGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCg                                   >  2:146185/1‑90 (MQ=255)                                                    ttcttcAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTa                             >  2:591784/1‑90 (MQ=255)                                                        tcAGTTGTTCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAAtt                         <  2:559471/90‑1 (MQ=255)                                                           gttgttCGAGAATGCGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCt                      >  2:562039/1‑90 (MQ=255)                                                                         cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc        <  1:347172/90‑1 (MQ=255)                                                                         cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc        <  2:695576/90‑1 (MQ=255)                                                                         cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc        <  2:723610/90‑1 (MQ=255)                                                                         cGTTTTACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc        <  1:205675/90‑1 (MQ=255)                                                                         cGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATcc        <  1:835543/90‑1 (MQ=255)                                                                              tACACCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGc   <  2:105173/90‑1 (MQ=255)                                                                               acacCCTGTTCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCa  <  1:320200/90‑1 (MQ=255)                                                                                       |                                                                                  GCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGCTGTTTGACTACCGGTATTAAATTGCTTCTGATATTGATCCGGTGCA  >  NZ_CP009273/1575739‑1575906

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