Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,998,494	C→T	*329* (TAG→TAA)	dcyD ←	D‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,998,494	0	C	T	100.0%	113.4 / NA	37	*329* (TAG→TAA)	dcyD	D‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/18);  total (19/18)

CGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACA  >  NC_000913/1998359‑1998626                                                                                                                                        |                                                                                                                                     cgAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTaa                                                                                                                                   <  1:1044491/139‑1 (MQ=255)                      aGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAAc                                                                                                              >  1:803349/1‑139 (MQ=255)                         cAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAgc                                                                                                           >  2:597459/1‑139 (MQ=255)                                     cgccgATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccg                                                                                               <  1:597459/139‑1 (MQ=255)                                           tACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTGGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                         <  2:630840/139‑1 (MQ=255)                                                  ttGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTgggg                                                                                                                           >  2:805504/1‑98 (MQ=255)                                                  ttGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTgggg                                                                                                                           <  1:805504/98‑1 (MQ=255)                                                                  gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGt                                                                                          <  1:699336/115‑1 (MQ=255)                                                                  gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGt                                                                                          >  2:699336/1‑115 (MQ=255)                                                                  gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCtt                                                                   <  1:583787/138‑1 (MQ=255)                                                                  gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCtt                                                                   >  2:583787/1‑138 (MQ=255)                                                                  gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTc                                                                  >  1:757314/1‑139 (MQ=255)                                                                  gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTc                                                                  >  1:3882/1‑139 (MQ=255)                                                                   ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccg                                                                                               <  2:1068415/109‑1 (MQ=255)                                                                   ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccg                                                                                               >  1:1068415/1‑109 (MQ=255)                                                                   ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                           >  2:39098/1‑113 (MQ=255)                                                                   ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                           >  2:607959/1‑113 (MQ=255)                                                                   ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                           <  1:607959/113‑1 (MQ=255)                                                                   ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                           <  1:39098/113‑1 (MQ=255)                                                                   ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCa                                                                 <  2:803349/139‑1 (MQ=255)                                                                                 aaCGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGttt                                                   <  2:3882/139‑1 (MQ=255)                                                                                                    ccAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATc                                >  1:895641/1‑139 (MQ=255)                                                                                                     cAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCa                               >  2:709745/1‑139 (MQ=255)                                                                                                       gTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCa                                                                 <  1:942482/103‑1 (MQ=255)                                                                                                       gTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCa                                                                 >  2:942482/1‑103 (MQ=255)                                                                                                       gTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGc                             <  1:709745/139‑1 (MQ=255)                                                                                                         tGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATc                                                                          >  2:1208531/1‑92 (MQ=255)                                                                                                         tGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATc                                                                          <  1:1208531/92‑1 (MQ=255)                                                                                                                  tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGc                  >  1:908914/1‑139 (MQ=255)                                                                                                                  tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGc                  >  1:803474/1‑139 (MQ=255)                                                                                                                      gCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt              <  2:757314/139‑1 (MQ=255)                                                                                                                           aGCGTTTGCTGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt         <  2:1257379/139‑1 (MQ=255)                                                                                                                             cGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt         <  1:903023/137‑1 (MQ=255)                                                                                                                             cGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt         >  2:903023/1‑137 (MQ=255)                                                                                                                                 tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCtt                                                                   >  2:945153/1‑75 (MQ=255)                                                                                                                                 tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCtt                                                                   <  1:945153/75‑1 (MQ=255)                                                                                                                                  ccGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAACCCCGCCATCGCTTTTCCGGTATAcaca  >  2:559156/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                     CGAGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACA  >  NC_000913/1998359‑1998626

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