Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,479,358	C→A	H53N (CAT→AAT)	dsdA →	D‑serine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,479,358	0	C	A	100.0%	168.9 / NA	48	H53N (CAT→AAT)	dsdA	D‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/25);  total (23/25)

CTCGCTCATCGCCCAGTATCCGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCCATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCAAAAACGGCTGGAAAAAGAATATCAG  >  NC_000913/2479222‑2479489                                                                                                                                         |                                                                                                                                    cTCGCTCATCGCCCAGTATCCGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATg                                                                                                                                  <  1:510910/140‑1 (MQ=255)   cGCTCATCGCCCAGTATCCGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATgcg                                                                                                                                <  1:145914/140‑1 (MQ=255)   cGCTCATCGCCCAGTATCCGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATgcg                                                                                                                                <  1:510513/140‑1 (MQ=255)    gCTCATCGCCCAGTATCCGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATgcgc                                                                                                                               >  1:1067020/1‑140 (MQ=255)    gCTCATCGCCCAGTATCCGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATgcgc                                                                                                                               >  1:1002545/1‑140 (MQ=255)              cAGTATCCGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCg                                                                                                                     <  2:894318/140‑1 (MQ=255)                     cGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCa                                                                                                              <  2:629692/140‑1 (MQ=255)                           tAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTAt                                                                                                        >  1:813673/1‑140 (MQ=255)                           tAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTAt                                                                                                        >  1:550447/1‑140 (MQ=255)                                   cTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCaaa                                                                                                <  2:801656/140‑1 (MQ=255)                                    tGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCaaaa                                                                                               >  1:871608/1‑140 (MQ=255)                                     ggTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAg                                                                                              >  2:378384/1‑140 (MQ=255)                                         gCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCAtt                                                                                          <  2:1002545/140‑1 (MQ=255)                                         gCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCAtt                                                                                          <  2:1067020/140‑1 (MQ=255)                                          ctctTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCAttt                                                                                         >  1:517725/1‑140 (MQ=255)                                                     accaccTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAActgc                                                                              <  2:550447/140‑1 (MQ=255)                                                     accaccTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAActgc                                                                              <  2:517725/140‑1 (MQ=255)                                                     accaccTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAActgc                                                                              >  2:382616/1‑140 (MQ=255)                                                     accaccTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAActgc                                                                              <  2:850090/140‑1 (MQ=255)                                                         ccTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGcc                                                                          >  2:917152/1‑140 (MQ=255)                                                         ccTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGcc                                                                          >  1:747762/1‑140 (MQ=255)                                                          cTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAActgct                                                                             <  2:240574/136‑1 (MQ=255)                                                          cTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAActgct                                                                             >  1:240574/1‑136 (MQ=255)                                                          cTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCa                                                                         >  1:952321/1‑140 (MQ=255)                                                                  aTCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACt                                                                       >  1:956764/1‑134 (MQ=255)                                                                  aTCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACt                                                                       <  2:956764/134‑1 (MQ=255)                                                                    ccTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGat                                                               >  2:308713/1‑140 (MQ=255)                                                                       ggCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGattat                                                            >  2:72818/1‑140 (MQ=255)                                                                       ggCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGattat                                                            <  2:998769/140‑1 (MQ=255)                                                                            gACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAAt                                                       <  2:1082424/140‑1 (MQ=255)                                                                             aCCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATc                                                      >  1:249349/1‑140 (MQ=255)                                                                               cTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAg                                                    <  2:871608/140‑1 (MQ=255)                                                                                 cATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGaa                                                  <  1:917152/140‑1 (MQ=255)                                                                                        tGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTg                                           >  2:403617/1‑140 (MQ=255)                                                                                           aGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGt                                             >  1:858447/1‑135 (MQ=255)                                                                                           aGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGt                                             <  2:858447/135‑1 (MQ=255)                                                                                              tttACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTc                                     <  1:72818/140‑1 (MQ=255)                                                                                                 aCCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAg                                  <  2:318109/140‑1 (MQ=255)                                                                                                      aTGTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGATGCCATTCCAGCTATg                             >  1:1160123/1‑140 (MQ=255)                                                                                                        gTTGGCCTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCa                           >  2:357125/1‑140 (MQ=255)                                                                                                             ccTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCaaaaa                       >  2:284350/1‑139 (MQ=255)                                                                                                             ccTGACCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCaaaaa                       <  1:284350/139‑1 (MQ=255)                                                                                                                 aCCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCAAAAACGGCt                  <  1:120286/140‑1 (MQ=255)                                                                                                                 aCCGAACAGGATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCAAAAACGGCt                  >  1:866561/1‑140 (MQ=255)                                                                                                                          gATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCAAAAACGGCTGGAAAAAGa         <  1:382616/140‑1 (MQ=255)                                                                                                                          gATGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCAAAAACGGCTGGAAAAAGa         <  1:403617/140‑1 (MQ=255)                                                                                                                           aTGTTCAGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCAAAAACGGCTGGAAAAAGaa        >  2:692938/1‑140 (MQ=255)                                                                                                                                 aGGACGCCAATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCAAAAACGGCTGGAAAAAGAATATCAg  <  1:308713/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    CTCGCTCATCGCCCAGTATCCGTTGGTAAAGGATCTGGTTGCTCTTAAAGAAACCACCTGGTTTAATCCTGGCACGACCTCATTGGCTGAAGGTTTACCTTATGTTGGCCTGACCGAACAGGATGTTCAGGACGCCCATGCGCGCTTATCCCGTTTTGCACCCTATCTGGCAAAAGCATTTCCTGAAACTGCTGCCACTGGGGGGATTATTGAATCAGAACTGGTTGCCATTCCAGCTATGCAAAAACGGCTGGAAAAAGAATATCAG  >  NC_000913/2479222‑2479489

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