Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,480,361	C→T	A387V (GCA→GTA)	dsdA →	D‑serine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,480,361	0	C	T	96.3%	73.2 / ‑3.7	28	A387V (GCA→GTA)	dsdA	D‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (5/21);  total (5/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

AACGCACTCCCCTTGTATGTTGTTAGGCGTCCATACAGGATTACACGATCAGATTTCTGTTCAGGATATTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACGCAGCATCGCAATCCTTTCCCTGGGTGAGCGATGCTGCCGATGGCGCAGACTTAAGATCCCCG  >  NC_000913/2480119‑2480605                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         aaCGCACTCTCCTTGTATGTTGTTAGGCGTCCATAAAGGATTACACGATCAGATTTCTGTTCAGGATATTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCgg                                                                                                                                                                                                                                                  >  2:560456/1‑251 (MQ=255) aaCGCACTCCCCTTGTATGTTGTTAGGCGTCCATACAGGATTACACGATCAGATTTCTGTTCAGGATATTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGACGGTATTCGTCTTGAACCTTCGGTACTGGCgg                                                                                                                                                                                                                                                  >  1:82517/1‑251 (MQ=255)       cTCCCCTTGTATGTTGTTAGGCGTCCATTCAGTAATACTCCATCAGATTTCTGTTCAGGATTTTTGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGTCTTTGTCTGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATTCCCTTATCGATCAATCCATGTGTGACATGCTTGTCTGGTTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTggttgttat                                                                                                                                                                                                                                              <  2:598118/249‑8 (MQ=255)         ccccTTGTATGTTGTTAGGCGTCCATACAGGATTACACGATCAGATTTCTGTTCAGGATATTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATgg                                                                                                                                                                                                                                            <  1:579453/249‑1 (MQ=255)          cccTTGTATGTTGTTAGGCGTCCATACAGGATTACACGATCAGATTTCTGTTCAGGATATTGGTATCGCCACCCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGGATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCCGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGCCATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGGCTTGAACCTTCGGTACTGGCAGGTATGGCCg                                                                                                                                                                                                                                         <  1:535244/251‑1 (MQ=255)           ccTTGTATGTTGTTAGGCGTCCATACAGGATTACACGATCAGATTTCTGTTCAGGATATTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCggagatggggtg                                                                                                                                                                                                                                        >  2:168902/1‑241 (MQ=255)                          ggCGTCCATACAGGATTACCCGATCAGATTTCTGTTCAGGATATTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTGTGTCGGGCTGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGTCATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCgtgtg                                                                                                                                                                                                                         <  2:282818/251‑1 (MQ=255)                                                              cAGGATATTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCATTAAGTTACCAACAGATGCACGG‑TTTCAg                                                                                                                                                                                    <  2:394969/251‑1 (MQ=255)                                                                  atatTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTGCAGCGCa                                                                                                                                                                                <  2:539651/251‑1 (MQ=255)                                                                          aTCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGaacaac                                                                                                                                                                         <  1:137786/250‑1 (MQ=255)                                                                            cGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGc                                                                                                                                                                      <  2:131628/251‑1 (MQ=255)                                                                                            gCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGTCCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCt                                                                                                                                                     <  1:613647/251‑1 (MQ=255)                                                                                                  ggCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGgtgt                                                                                                                                                <  1:140819/250‑1 (MQ=255)                                                                                                                                gTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCCATGCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTTAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCATTAAGTTCCCAACAGATGCACTT‑TTTCATCGCAGAACAACTGCGGAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGgaaga                                                                                                                 <  1:591067/251‑1 (MQ=255)                                                                                                                                      cgggcAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGc                                                                                                                        <  1:263273/238‑1 (MQ=255)                                                                                                                                             gggAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACTTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCTGTAGAAGAGA‑TGAATCaa                                                                                                   <  1:465485/250‑1 (MQ=255)                                                                                                                                                            aTGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGGATGGGCGGACCCCCGCGCGGGGGGGCGGCAGTAAGTTAACAACAGATGCACGGGTTTC‑GCGCCGAAAAAAAGCGGAATACCACCCCAT‑TGGTGTGGGGGACGGGAGGGGGGATGGGGCCGGGAGAAG‑GAGTGAATCAATATTTGGCaaaaagg                                                                                    >  1:483153/1‑248 (MQ=255)                                                                                                                                                                     ataCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTtaataa                                                                           <  1:474507/251‑1 (MQ=255)                                                                                                                                                                             aGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑Caa                                                                  <  1:189550/251‑1 (MQ=255)                                                                                                                                                                                gATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACgc                                                               <  1:573076/251‑1 (MQ=255)                                                                                                                                                                                 aTCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACgca                                                              <  1:167229/251‑1 (MQ=255)                                                                                                                                                                                      aCCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGTACTGGCGGGGATGGACGGAGCTCACCGCCGGTGTGCATCAGTAAATTTTCAACACATGAACGG‑GTTCCGCGCCGGAAAAATGCGGCATACAACAC‑ATCCGGGGGGGGCGGCGGGAGGGGGGAAGGGGGCGGGAGGAGAGA‑AGAAACAATAACTGGCAAAAGGCCGTTAAAAACGTTTTCAA‑GCAGCATc                                                         >  1:322417/1‑251 (MQ=255)                                                                                                                                                                                                                    cAGGAAGAAGGTATTCGTCTTGAACCTTCGGAACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACGCAGCATCGCAATCCTTTCCCTGGGTGAGCGAtgctgc                           <  1:220458/251‑1 (MQ=255)                                                                                                                                                                                                                         agaagGTATTCGTCTTGAACCTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACGCAGCATCGCAATCCTTTCCCTGGGTGAGCGATGCTGCCGATg                      <  1:342213/251‑1 (MQ=255)                                                                                                                                                                                                                              gTATTCGTCTTGAACCTTCGGTACTGTCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACGCAGCATCGCAATCCTTTCCCTGGGTGAGCGATGCTGCCGATGGCGCa                 <  2:529106/251‑1 (MQ=255)                                                                                                                                                                                                                                         gagccTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACGCAGCATCGCAATCCTTTCCCTGGGTGAGCGATGCTGCCGATGGCGCAGACGTAAGAt       <  2:82517/247‑1 (MQ=255)                                                                                                                                                                                                                                            ccTTCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACGCAGCATCGCAATCCTTTCCCTGGGTGAGCGATGCTGCCGATGGCGCAGACTTAAGATcccc   <  1:473240/251‑1 (MQ=255)                                                                                                                                                                                                                                               tCGGTACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACGCAGCATCGCAATCCTTTCCCTGGGTGAGCGATGCTGCCGATGGCGCAGACTTAAGATCCCCg  <  2:283903/249‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         AACGCACTCCCCTTGTATGTTGTTAGGCGTCCATACAGGATTACACGATCAGATTTCTGTTCAGGATATTGGTATCGACAACCTTACCGCAGCGGATGGCCTTGCAGTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGG‑TTTCAGCGCAGAACAACTGCGTAATACCACTC‑ATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGA‑TGAATCAATATCTGGCAAAAGGCCGTTAATAACGTTT‑CAACGCAGCATCGCAATCCTTTCCCTGGGTGAGCGATGCTGCCGATGGCGCAGACTTAAGATCCCCG  >  NC_000913/2480119‑2480605

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