Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,255,685	A→G	C323R (TGT→CGT)	cyuA ←	putative L‑cysteine desulfidase CyuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,255,685	0	A	G	100.0%	157.5 / NA	45	C323R (TGT→CGT)	cyuA	putative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/24);  total (21/24)

CGACGCACCATCGCAAATCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACACAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAGTGTTCTGCTACCACCACCACAGGCATTGTTGCGGTA  >  NC_000913/3255551‑3255818                                                                                                                                       |                                                                                                                                      cGACGCACCATCGCAAATCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAgcg                                                                                                                                  <  2:257897/140‑1 (MQ=255)   aCGCACCATCGCAAATCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGTGGCACGCAGCGCa                                                                                                                                <  2:757790/140‑1 (MQ=255)      cACCATCGCAAATCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGAc                                                                                                                             >  1:809628/1‑140 (MQ=255)            cGCAAATCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAAcgc                                                                                                                       >  1:285393/1‑140 (MQ=255)              cAAATCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGACACGCAGCGCAGACAAACGCgg                                                                                                                     >  2:151614/1‑140 (MQ=255)                 aTCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTaa                                                                                                                  <  2:776388/140‑1 (MQ=255)                 aTCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTaa                                                                                                                  <  2:676943/140‑1 (MQ=255)                        cGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTa                                                                                                           <  1:369611/140‑1 (MQ=255)                        cGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTa                                                                                                           <  2:107091/140‑1 (MQ=255)                                       tCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTa                                                                                                           <  2:694862/125‑1 (MQ=255)                                       tCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTa                                                                                                           >  1:694862/1‑125 (MQ=255)                                            ctgctgATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCaa                                                                                       <  1:125136/140‑1 (MQ=255)                                                 gATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGc                                                                                   <  2:411502/139‑1 (MQ=255)                                                 gATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGc                                                                                   >  1:411502/1‑139 (MQ=255)                                                  aTCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGTGa                                                                                 >  2:562174/1‑140 (MQ=255)                                                  aTCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGa                                                                                 >  2:780389/1‑140 (MQ=255)                                                  aTCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGa                                                                                 <  1:780389/140‑1 (MQ=255)                                                   tCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGaa                                                                                <  1:47542/140‑1 (MQ=255)                                                    cGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGaaa                                                                               <  2:285393/140‑1 (MQ=255)                                                     gCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAg                                                                              >  2:435920/1‑140 (MQ=255)                                                       cATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCa                                                                            >  1:780805/1‑140 (MQ=255)                                                        aTCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCAt                                                                           <  1:103261/140‑1 (MQ=255)                                                             gaTGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAgcg                                                                      <  1:618893/140‑1 (MQ=255)                                                              aTGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAgcgc                                                                     >  1:53787/1‑140 (MQ=255)                                                              aTGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAgcgc                                                                     >  1:368268/1‑140 (MQ=255)                                                               tGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCa                                                                    <  1:232477/140‑1 (MQ=255)                                                                    tCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACgcgc                                                               <  2:571527/140‑1 (MQ=255)                                                                       tAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAg                                                            >  2:629486/1‑140 (MQ=255)                                                                           gCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGt                                                        >  2:686284/1‑140 (MQ=255)                                                                            cccATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACgcgc                                                               >  2:33373/1‑132 (MQ=255)                                                                            cccATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACgcgc                                                               <  1:33373/132‑1 (MQ=255)                                                                                  caccaGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTcatcat                                                 <  2:53787/140‑1 (MQ=255)                                                                                     caGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCg                                              >  2:520924/1‑140 (MQ=255)                                                                                          atgccatCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCg                                         <  1:562174/140‑1 (MQ=255)                                                                                           tgccatCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGa                                        >  1:95855/1‑140 (MQ=255)                                                                                           tgccatCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGa                                        >  1:746304/1‑140 (MQ=255)                                                                                           tgccatCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGa                                        >  2:533192/1‑140 (MQ=255)                                                                                              catCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAgt                                     >  2:83092/1‑140 (MQ=255)                                                                                               atCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAgtg                                    <  1:686284/140‑1 (MQ=255)                                                                                                 cccGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAGTGtt                                  <  1:629486/140‑1 (MQ=255)                                                                                                  ccGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAGTGTTc                                 >  2:804822/1‑140 (MQ=255)                                                                                                        gcggcCCCCATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAGTGTTCTGCTac                           <  1:533192/140‑1 (MQ=255)                                                                                                              cccATTGCTGCGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAGTGTTCTGCTaccaccac                     <  2:368268/140‑1 (MQ=255)                                                                                                                        cGGTCGTTGCGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAGTGTTCTGCTACCACCACCACAGGCAtt           >  2:187607/1‑140 (MQ=255)                                                                                                                                 cGGCACGCAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAGTGTTCTGCTACCACCACCACAGGCATTGTTGCGGTa  <  2:780805/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      CGACGCACCATCGCAAATCATGCCGCTGACATCGCCGATCATACTGCTGATCGCCATCGAGATGGTTTCATAACGCCCATCCACCAGCCATGCCATCCCGGCGGCGGCCCCCATTGCTGCGGTCGTTGCGGCACACAGCGCAGACAAACGCGGTAACTGGTTATGGATGTAAATTGCGCTCAAATGCGAAAGCATCAGCGCACGCGCCAGCCGTTCATCATCCGCTCCGAAGTGTTCTGCTACCACCACCACAGGCATTGTTGCGGTA  >  NC_000913/3255551‑3255818

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