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

AGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGT  >  NC_000913/1998380‑1998619                                                                                                                   |                                                                                                                              aGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAAc                                                                                                       <  1:755919/139‑1 (MQ=255)                 gctgATACTGAGGTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAATGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                       <  2:233911/136‑1 (MQ=255)                   cgATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgccg                                                                                     >  2:943645/1‑139 (MQ=255)                    gATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                    <  2:118203/139‑1 (MQ=255)                      tACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                  >  2:287817/1‑139 (MQ=255)                           aGTTGCAGCGTATACCCAGCCCCTTTCAGTAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAg                                                                                                     >  2:1326760/1‑115 (MQ=255)                           aGTTGCAGCGTATACCCAGCCCCTTTCAGTAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAg                                                                                                     <  1:1326760/115‑1 (MQ=255)                                    gTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAAt                                                                    <  1:943645/139‑1 (MQ=255)                                     tataCCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGAAGTTAAACGTGGGgat                                                                                                                  >  2:1689891/1‑92 (MQ=255)                                     tataCCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGAAGTTAAACGTGGGgat                                                                                                                  <  1:1689891/92‑1 (MQ=255)                                          ccAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGccc                                                              <  2:1156108/139‑1 (MQ=255)                                             gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGa                                                                                                         <  2:15244/93‑1 (MQ=255)                                             gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGa                                                                                                         >  1:15244/1‑93 (MQ=255)                                             gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                    >  1:1110396/1‑114 (MQ=255)                                             gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                    <  2:1110396/114‑1 (MQ=255)                                              ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGTGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                    >  2:1119443/1‑113 (MQ=255)                                              ccccTTTCAGCAGGAACGGCAAGGAAGCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                    <  1:1119443/113‑1 (MQ=255)                                               cccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCa                                                                                           >  1:1218592/1‑105 (MQ=255)                                               cccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCa                                                                                           <  2:1218592/105‑1 (MQ=255)                                               cccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTGGTATAATTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCAt                                                         >  1:1133050/1‑139 (MQ=255)                                                      agcagGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGaa                                                  <  1:287817/139‑1 (MQ=255)                                                                 cAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                       >  1:433717/1‑91 (MQ=255)                                                                 cAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                       <  2:433717/91‑1 (MQ=255)                                                                  aaGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAAt                                                                    >  2:1216390/1‑109 (MQ=255)                                                                  aaGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAAt                                                                    <  1:1216390/109‑1 (MQ=255)                                                                        tCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATAcc                                <  1:1627943/139‑1 (MQ=255)                                                                         cAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCg                               >  2:487663/1‑139 (MQ=255)                                                                         cAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGAAGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCg                               >  1:429936/1‑139 (MQ=255)                                                                                cAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAg                                                 <  1:1339929/114‑1 (MQ=255)                                                                                cAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAg                                                 >  2:1339929/1‑114 (MQ=255)                                                                                  gTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                       >  1:820604/1‑74 (MQ=255)                                                                                  gTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccgc                                                                                       <  2:820604/74‑1 (MQ=255)                                                                                  gTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAAt                          <  1:1414/135‑1 (MQ=255)                                                                                  gTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAAt                          >  2:1414/1‑135 (MQ=255)                                                                                  gTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGc                      >  2:51937/1‑139 (MQ=255)                                                                                  gTTGTATACTTTCTTGCATGAGCGTTTACGGAGTTAAACGTg                                                                                                                       >  1:556544/1‑42 (MQ=255)                                                                                  gTTGTATACTTTCTTGCATGAGCGTTTACGGAGTTAAACGTg                                                                                                                       <  2:556544/42‑1 (MQ=255)                                                                                   ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAAc                                                                                                                          >  1:1274542/1‑38 (MQ=255)                                                                                   ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAAc                                                                                                                          <  2:1274542/38‑1 (MQ=38)                                                                                      tataCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATc                                                                   <  2:690174/90‑1 (MQ=255)                                                                                      tataCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATc                                                                   >  1:690174/1‑90 (MQ=255)                                                                                        taCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAg                                                                                                     >  2:1346760/1‑54 (MQ=255)                                                                                        taCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAg                                                                                                     <  1:1346760/54‑1 (MQ=255)                                                                                        taCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGcc                <  1:487663/139‑1 (MQ=255)                                                                                             tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                        <  2:1438952/46‑1 (MQ=255)                                                                                             tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                        >  1:1438952/1‑46 (MQ=255)                                                                                             tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGgcg                                                                                             >  2:753428/1‑57 (MQ=255)                                                                                             tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGgcg                                                                                             <  1:753428/57‑1 (MQ=255)                                                                                             tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCtt                                                            <  2:1182306/90‑1 (MQ=255)                                                                                             tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCtt                                                            >  1:1182306/1‑90 (MQ=255)                                                                                               ttGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCCGCCCCGCCATCGCtt         <  2:1133050/139‑1 (MQ=255)                                                                                                tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCttt        >  2:1046847/1‑139 (MQ=255)                                                                                                     gAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGgcg                                                                                             <  2:803089/49‑1 (MQ=255)                                                                                                     gAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGgcg                                                                                             >  1:803089/1‑49 (MQ=255)                                                                                                     gAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt       >  1:1559422/1‑135 (MQ=255)                                                                                                     gAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt       <  2:1559422/135‑1 (MQ=255)                                                                                                      aGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt  >  1:291694/1‑139 (MQ=255)                                                                                                      aGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt  >  1:1003963/1‑139 (MQ=255)                                                                                                      aGCGTTTCCGGAGTTAAACGTGGGGAAGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt  <  2:429936/139‑1 (MQ=255)                                                                                                                   |                                                                                                                              AGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGT  >  NC_000913/1998380‑1998619

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: ‑