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

AGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTA  >  NC_000913/1998361‑1998620                                                                                                                                      |                                                                                                                               agaATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAAc                                                                                                                           <  2:404736/139‑1 (MQ=255)             aGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATTCCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATa                                                                                                               >  1:10478/1‑139 (MQ=255)                          aaaaCATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCa                                                                                                  >  2:480199/1‑139 (MQ=255)                              cATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGTGAACAGCGCAGgcg                                                                                              <  1:224631/139‑1 (MQ=255)                                tGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAAc                                                                                                                           <  1:1223922/108‑1 (MQ=255)                                tGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAAc                                                                                                                           >  2:1223922/1‑108 (MQ=255)                                       gATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                     <  2:1483844/139‑1 (MQ=255)                                               gTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATa                                                                             >  2:1403932/1‑139 (MQ=255)                                                                gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                   <  1:112961/116‑1 (MQ=255)                                                                gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                   >  2:112961/1‑116 (MQ=255)                                                                 ccccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCa                                                           >  2:644300/1‑139 (MQ=255)                                                                  cccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCAt                                                          >  2:873224/1‑139 (MQ=255)                                                                     tttCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCtt                                                       >  2:682820/1‑139 (MQ=255)                                                                          gcagGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAg                                                  <  1:682820/139‑1 (MQ=255)                                                                                         aaTCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                   <  2:86771/91‑1 (MQ=255)                                                                                         aaTCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                   >  1:86771/1‑91 (MQ=255)                                                                                           tCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATAcc                                 <  1:873224/139‑1 (MQ=255)                                                                                            cAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCg                                <  1:1403932/139‑1 (MQ=255)                                                                                                      ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                         >  2:1356317/1‑56 (MQ=255)                                                                                                      ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                         <  1:1356317/56‑1 (MQ=255)                                                                                                      ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGTGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCaa                            >  2:720826/1‑133 (MQ=255)                                                                                                      ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGTGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCaa                            <  1:720826/133‑1 (MQ=255)                                                                                                      ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                           >  2:536942/1‑86 (MQ=255)                                                                                                      ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATaaa                                                                           <  1:536942/86‑1 (MQ=255)                                                                                                      ttGTATACTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGcc                      >  2:380680/1‑139 (MQ=255)                                                                                                         tataCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccg                                                                                         >  1:1150968/1‑69 (MQ=255)                                                                                                         tataCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAccg                                                                                         <  2:1150968/69‑1 (MQ=255)                                                                                                         tataCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAcc                                                                                          <  1:1138136/68‑1 (MQ=255)                                                                                                         tataCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCAcc                                                                                          >  2:1138136/1‑68 (MQ=255)                                                                                                           taCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCtt                                                             >  2:501/1‑95 (MQ=255)                                                                                                           taCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCtt                                                             <  1:501/95‑1 (MQ=255)                                                                                                           taCTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGcc                 >  2:1192775/1‑139 (MQ=255)                                                                                                             cTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACa                                                                         <  2:1265458/81‑1 (MQ=255)                                                                                                             cTTTCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACa                                                                         >  1:1265458/1‑81 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCg                                                                                                           <  1:1042371/44‑1 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCg                                                                                                           >  2:1042371/1‑44 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                         <  1:490693/46‑1 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGaa                                                                                                         >  2:490693/1‑46 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAAc                                                                                                        >  1:37604/1‑47 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAAc                                                                                                        <  2:37604/47‑1 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                     >  1:411979/1‑66 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                     <  2:411979/66‑1 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                   <  1:194509/68‑1 (MQ=255)                                                                                                                tCTTGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTa                                                                                   >  2:194509/1‑68 (MQ=255)                                                                                                                   tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                     >  2:448736/1‑63 (MQ=255)                                                                                                                   tGCATGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCgg                                                                                     <  1:448736/63‑1 (MQ=255)                                                                                                                   tGCATGAGCGTTTCCGGAGTTAAAC‑‑GGGGATGATAGGCGAACAGCGCAGgcg                                                                                              >  2:515394/1‑52 (MQ=255)                                                                                                                   tGCATGAGCGTTTCCGGAGTTAAAC‑‑GGGGATGATAGGCGAACAGCGCAGgcg                                                                                              <  1:515394/52‑1 (MQ=255)                                                                                                                       tGAGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCg     <  2:10478/139‑1 (MQ=255)                                                                                                                         aGCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt   <  1:644300/139‑1 (MQ=255)                                                                                                                          gCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTc       >  2:1345661/1‑134 (MQ=255)                                                                                                                          gCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTc       <  1:1345661/134‑1 (MQ=255)                                                                                                                          gCGTTTCCGGAGTTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGta  <  1:1192775/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                               AGAATAAACCCGAGCAGTAAGCCAAAAAACATGCCGCCGATACTGAGTTGCAGCGTATACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTA  >  NC_000913/1998361‑1998620

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