Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,998,555	G→A	P309L (CCG→CTG)	dcyD ←	D‑cysteine desulfhydrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,998,555	0	G	A	100.0%	164.6 / NA	49	P309L (CCG→CTG)	dcyD	D‑cysteine desulfhydrase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/19);  total (30/19)

TACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTCCATGCCTTCGTCGTTCG  >  NC_000913/1998418‑1998690                                                                                                                                          |                                                                                                                                        taCCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGGTAGGCGAAAAGCCCAGGCGCACCGCGGGTATGAATAAACAGAATCAGc                                                                                                                                       >  2:321298/1‑140 (MQ=255) taCCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGc                                                                                                                                       >  2:613581/1‑140 (MQ=255)       gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCa                                                                                                                                 >  2:252360/1‑140 (MQ=255)       gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCa                                                                                                                                 >  2:240862/1‑140 (MQ=255)       gCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCa                                                                                                                                 >  2:119717/1‑140 (MQ=255)         cccTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCGCCGTCGGTATGAGTAAACAGAATCAGCCCTTCATc                                                                                                                               >  2:294959/1‑140 (MQ=255)               cagcagGAACGGCAAAGAATCACTACCCAATTGTATCCTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGaa                                                                                                                         <  1:347585/140‑1 (MQ=255)                agcggGACCGGCAAGGAATCAATAACCAGTTGTATTCTTTCTGGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAg                                                                                                                        <  1:709375/140‑1 (MQ=255)                agcagGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAg                                                                                                                        <  2:1572/140‑1 (MQ=255)                     gAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAg                                                                                                                        >  2:873951/1‑135 (MQ=255)                     gAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAg                                                                                                                        <  1:873951/135‑1 (MQ=255)                     gAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGttt                                                                                                                   <  1:252360/140‑1 (MQ=255)                                           aGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCCCCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGc                                                                                             <  2:841136/140‑1 (MQ=255)                                             ttGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGccc                                                                                           >  1:552682/1‑140 (MQ=255)                                                tataCTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGc                                                                                        >  2:234574/1‑140 (MQ=255)                                                tataCTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGc                                                                                        >  2:153024/1‑140 (MQ=255)                                                  taCTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCa                                                                                      >  2:610336/1‑140 (MQ=255)                                                    cTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATc                                                                                    <  1:592686/140‑1 (MQ=255)                                                       tCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCt                                                                                 >  1:714437/1‑140 (MQ=255)                                                       tCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCt                                                                                 >  2:694544/1‑140 (MQ=255)                                                          tGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGTGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt                                                                              >  2:604932/1‑140 (MQ=255)                                                          tGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt                                                                              <  1:180201/140‑1 (MQ=255)                                                          tGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt                                                                              >  1:187510/1‑140 (MQ=255)                                                          tGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt                                                                              <  2:187510/140‑1 (MQ=255)                                                          tGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt                                                                              <  1:694544/140‑1 (MQ=255)                                                          tGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt                                                                              >  2:712191/1‑140 (MQ=255)                                                          tGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTCCTGACTGTTACTGTGAATCAGGCCCGCCATCGCtttt                                                                              >  1:310388/1‑140 (MQ=255)                                                              tGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCgg                                                                          >  2:205211/1‑140 (MQ=255)                                                               gggCGTTTCGGGAGCTAGACGGTGGGAGGATGGGCTAACCGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt                                                                         <  2:310388/138‑1 (MQ=255)                                                                ggCGTTTCGGGGGCTAAACGGGGGGATGATAGGCGAAAAGCGGAGGCCCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGta                                                                        <  1:9784/139‑1 (MQ=255)                                                                aGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGta                                                                        >  2:111447/1‑140 (MQ=255)                                                                  ggTTTCCGGAGCTAATCGGGGGCATCATACGCGAACACCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGtata                                                                      <  1:294959/139‑1 (MQ=255)                                                                    tttCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATaca                                                                    <  1:872680/140‑1 (MQ=255)                                                                    tttCCGGAGCTAAACCTGGGGATGAGAGGGGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATaca                                                                    <  1:400796/140‑1 (MQ=255)                                                                        cGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAgg                                                                <  2:552682/140‑1 (MQ=255)                                                                           aGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATc                                                             >  1:256698/1‑140 (MQ=255)                                                                                aaCGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGACCCGCCAACGCGTTGCCGGTATACACAGGATCAAGCa                                                        >  1:465928/1‑140 (MQ=255)                                                                                     ggggATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATg                                                   >  2:14111/1‑140 (MQ=255)                                                                                      gggATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGc                                                  >  1:243201/1‑140 (MQ=255)                                                                                       ggATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAAGCAGCCCCGGCATCGCTATTCCGGTAAACACAGGATCAAGCAGAATGcc                                                 >  2:274734/1‑140 (MQ=255)                                                                                           gatAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTcc                                             >  1:705387/1‑140 (MQ=255)                                                                                            atAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCa                                            >  1:1000424/1‑140 (MQ=255)                                                                                                     cAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGGAGAATGCCTTCCAGCCGCGCCa                                   >  2:338165/1‑140 (MQ=255)                                                                                                     cAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCa                                   >  2:753329/1‑140 (MQ=255)                                                                                                      aGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAg                                  <  1:14111/140‑1 (MQ=255)                                                                                                            ggCGCACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAAttt                            <  2:1000424/140‑1 (MQ=255)                                                                                                               gcACCGCCGGTATGAATAAACAGAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCAc                         <  2:705387/140‑1 (MQ=255)                                                                                                                                     gAATCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTCCATGCCTTCGTCGTTc   >  1:572106/1‑140 (MQ=255)                                                                                                                                      aaTCAGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGACTTCCAGCGGCGGCAGTAATTTCACAGCTTGCATGCCTTCGTCGTTCg  >  2:259588/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        TACCCAGCCCCTTTCAGCAGGAACGGCAAGGAATCAATAACCAGTTGTATACTTTCTTGCATGAGCGTTTCCGGAGCTAAACGTGGGGATGATAGGCGAACAGCGCAGGCGCACCGCCGGTATGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTCCATGCCTTCGTCGTTCG  >  NC_000913/1998418‑1998690

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