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%	199.4 / NA	62	*329* (TAG→TAA)	dcyD	D‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/34);  total (28/34)

GAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACAC  >  NC_000913/1998362‑1998625                                                                                                                                     |                                                                                                                                    gaATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACg                                                                                                                               <  1:1302039/139‑1 (MQ=255)    tAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTgg                                                                                                                            <  2:620503/139‑1 (MQ=255)     aaaCCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTggt                                                                                                                           <  2:819734/139‑2 (MQ=255)                gTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCg                                                                                                                >  1:1219874/1‑139 (MQ=255)                               tGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCa                                                                                                 >  1:1018874/1‑139 (MQ=255)                               tGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCa                                                                                                 <  2:1220305/139‑1 (MQ=255)                                  cgccgATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccg                                                                                              <  1:1103821/139‑1 (MQ=255)                                   gccgATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                             <  1:1296069/139‑1 (MQ=255)                                        tACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                        <  1:843855/139‑1 (MQ=255)                                        tACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                        <  1:718812/139‑1 (MQ=255)                                               ttGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaa                                                                                 >  2:150296/1‑139 (MQ=255)                                                tGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                <  2:1219874/139‑1 (MQ=255)                                                tGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                <  1:150296/139‑1 (MQ=255)                                                  cAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACa                                                                              >  1:1062959/1‑139 (MQ=255)                                                               gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACg                                                                                                                               <  1:24976/77‑1 (MQ=255)                                                               gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACg                                                                                                                               >  2:24976/1‑77 (MQ=255)                                                               gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                <  1:1009038/124‑1 (MQ=255)                                                               gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                                >  2:1009038/1‑124 (MQ=255)                                                               gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTc                                                                 >  1:233828/1‑139 (MQ=255)                                                                   cTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCt                                                             >  2:352318/1‑139 (MQ=255)                                                                     ttCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCttt                                                           >  1:1266460/1‑139 (MQ=255)                                                                          cagGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                              >  1:605481/1‑83 (MQ=255)                                                                          cagGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                              <  2:605481/83‑1 (MQ=255)                                                                          cagGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAg                                                       >  2:568156/1‑138 (MQ=255)                                                                          cagGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAg                                                       <  1:568156/138‑1 (MQ=255)                                                                             gAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                             >  1:627005/1‑97 (MQ=255)                                                                             gAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                             <  2:627005/97‑1 (MQ=255)                                                                                  gCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTctga                                              <  2:301982/139‑1 (MQ=255)                                                                                      ggAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTctgactga                                          <  2:233828/139‑1 (MQ=255)                                                                                                aCCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAAt                                <  1:244245/139‑1 (MQ=255)                                                                                                     ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATATGCGAACAGCGCAGGCGCACCGCCgg                                                                                          >  2:1260283/1‑76 (MQ=255)                                                                                                     ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATATGCGAACAGCGCAGGCGCACCGCCgg                                                                                          <  1:1260283/76‑1 (MQ=255)                                                                                                     ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                              <  2:578016/56‑1 (MQ=255)                                                                                                     ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                              >  1:578016/1‑56 (MQ=255)                                                                                                      tGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                          >  2:395529/1‑75 (MQ=255)                                                                                                      tGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                          <  1:395529/75‑1 (MQ=255)                                                                                                               tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                        >  2:200390/1‑68 (MQ=255)                                                                                                               tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                        <  1:200390/68‑1 (MQ=255)                                                                                                               tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTctga                                              >  2:97080/1‑110 (MQ=255)                                                                                                               tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTctga                                              <  1:97080/110‑1 (MQ=255)                                                                                                               tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGc                 >  1:838403/1‑139 (MQ=255)                                                                                                               tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGc                 >  2:781164/1‑139 (MQ=255)                                                                                                               tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGc                 >  1:285001/1‑139 (MQ=255)                                                                                                               tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGc                 >  2:1018029/1‑139 (MQ=255)                                                                                                                  tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                             >  1:801878/1‑60 (MQ=255)                                                                                                                  tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                             <  2:801878/60‑1 (MQ=255)                                                                                                                  tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGGATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCttt              >  1:819610/1‑139 (MQ=255)                                                                                                                  tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGttt                                                  <  1:929321/103‑1 (MQ=255)                                                                                                                  tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGttt                                                  >  2:929321/1‑103 (MQ=255)                                                                                                                  tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTct                                                <  2:600581/105‑1 (MQ=255)                                                                                                                  tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTct                                                >  1:600581/1‑105 (MQ=255)                                                                                                                   gCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt             <  2:1018874/139‑1 (MQ=255)                                                                                                                   gCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt             <  1:352318/139‑1 (MQ=255)                                                                                                                      tGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCg          <  2:1062959/139‑1 (MQ=255)                                                                                                                        aGCGTTTCTGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGGCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt        <  1:781164/139‑1 (MQ=255)                                                                                                                        aGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt        >  1:1300863/1‑139 (MQ=255)                                                                                                                        aGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt        <  2:285001/139‑1 (MQ=255)                                                                                                                        aGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt        <  2:1300863/139‑1 (MQ=255)                                                                                                                         gCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGtc       <  2:819610/139‑2 (MQ=255)                                                                                                                              tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGaa                                                        >  1:869276/1‑85 (MQ=255)                                                                                                                              tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGaa                                                        <  2:869276/85‑1 (MQ=255)                                                                                                                              tCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATacac  >  1:380507/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                    GAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACAC  >  NC_000913/1998362‑1998625

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