Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,984,524	G→T	T495K (ACA→AAA)	aslA ←	putative sulfatase AslA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,984,524	0	G	T	95.8%	74.9 / NA	24	T495K (ACA→AAA)	aslA	putative sulfatase AslA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (11/12);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGC  >  NC_000913/3984278‑3984769                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      gCAGCTTTCTGTAGGTTGGATAATATGCTTCAGCTTCGCATCCGGCAAAGGCAGATCTCAGCGATAGCCCCGGCTTATGCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCAAATACGCGTCCATTTCGGTCTGTAGGGGTCCACCCATCGGAATATGTCGCACGCCGATGGAGTCTCTTTCCTGCGGATCGGTGTAGAGGTTAAACCCCGACGATCCCGCCGTTTGCATTACTTTGcc                                                                                                                                                                                                                                                   <  1:78045/251‑1 (MQ=255)           gTAGGCTGGAAAAGTTGCGTCAGCAGCGAATCCGGAAAAGGCAGAGCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATTGGCGCGCGTGTTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTGCGGTCTGTAGCGGTCCACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCa                                                                                                                                                                                                                                         <  1:130006/251‑1 (MQ=255)                  ggATAAGATGCGTAAGCATCGCATACGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATGTAATCTGCGCGCGTGGGGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGat                                                                                                                                                                                                                                  >  2:172001/1‑251 (MQ=255)                                      gcatcCGGAAAAGGCAGATCTCGGCGATAGCGCGGGCTTAGTCAGATTTAATCTGCGCGGGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCTTTTCGGTCGGTAGCGGTACCCCCATCGGAATAGGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCg                                                                                                                                                                                                              <  2:159266/251‑1 (MQ=255)                                       catcCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGACGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGt                                                                                                                                                                                                             >  1:124470/1‑251 (MQ=255)                                                                  aGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACt                                                                                                                                                                                  >  1:90506/1‑251 (MQ=255)                                                                          cTTAGTCAGATTTAATCTGCGAGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTc                                                                                                                                                                           >  1:173915/1‑250 (MQ=255)                                                                                cAGATTTAATCTGCGCGGGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTTTAGCGGTACCCCCATCGAAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCAt                                                                                                                                                                     <  2:62560/250‑1 (MQ=255)                                                                                 aGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACCCCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATAc                                                                                                                                                                   <  1:210922/251‑1 (MQ=255)                                                                                                     gtggATATTTTTTCAGGATCTCCATATACGCGCGCATTTCGGTCTGACGCGGTACAACCATCGGAATATGGCGCCCGCAGATGGCGTCGCTTTCCTGCGGATAGGTGTAGAGGTTAAACAACGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCg                                                                                                                                               >  1:145031/1‑251 (MQ=255)                                                                                                                              ataCGCGGGCATTTCGGTTTGTAGCGGTACCCCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTGGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTAc                                                                                                                      <  2:96700/251‑1 (MQ=255)                                                                                                                              aacggtggtcatttgggttgggtggggTACCCCCAGCGGAATTGGGCGCACGACGAGGGAGTCGCTTGCCTGTGGACTGTTGTAGAGGTTAAACCCCGACGATCCCGCCGGTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTAc                                                                                                                      <  2:30924/226‑1 (MQ=255)                                                                                                                               aagggtggattttCGGTTGGTAGGGGTACACCCATCGGAATAGGGCGCACGCCGAGGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTAc                                                                                                                      <  2:206340/240‑1 (MQ=255)                                                                                                                                                  gaagccgtaaccccATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCATACGGCTCGGTGTAGAGGTTAAAAACCGACGCTCCAGACGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGAACGGTTGCTGAATCAGGACGTGATCCTTGAACTAATCCATACGCACAGCCGAGAGTTTACCGTTGAGGATGTATAGCTCGGCCTTACGGTTAGACTGCCCATTTGtt                                                                                                  >  2:159551/12‑251 (MQ=255)                                                                                                                                                                                    gATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCg                                                                 >  1:189549/1‑250 (MQ=255)                                                                                                                                                                                               tttCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtt                                                     >  1:138652/1‑251 (MQ=255)                                                                                                                                                                                                  ccTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtttt                                                   <  1:35583/250‑1 (MQ=255)                                                                                                                                                                                                            ggTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGTTGGTTTTCGTCACTAAAt                                        <  1:59039/251‑1 (MQ=255)                                                                                                                                                                                                                                 cgaTCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGAt                   >  2:143180/1‑251 (MQ=255)                                                                                                                                                                                                                                  gaTCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCAGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGCGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATACAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATg                  >  1:35674/1‑251 (MQ=255)                                                                                                                                                                                                                                   aTCCCGCCGTTTGCATTACGTTGCCGGTGAATCCACCCTGATATCCGCGCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTGACGGTTAGACTGACCATTTGGGCCTAGGAGGAAGGATGTCTGGTCCACGCCATCGATAAAGGTGGGTTTCGGCGCGAAATTCGCCACTGTCGCACCAGGAt                   <  2:118348/249‑1 (MQ=255)                                                                                                                                                                                                                                    tCCCGCCGTTTGCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCCGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGcc                >  1:216966/1‑251 (MQ=255)                                                                                                                                                                                                                                               gCATTACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCCGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATc      <  2:216966/250‑1 (MQ=255)                                                                                                                                                                                                                                                   tACTTTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGTTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGTATGCCCTGCCAGATCCAgc  <  2:173915/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      GCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATCGGAATATGGCGCACGCCGATGGAGTCGCTTTCCTGCGGATCGGTGTAGAGGTTAAACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGC  >  NC_000913/3984278‑3984769

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