Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,256,157	C→T	Q165Q (CAG→CAA)	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,256,157	0	C	T	100.0%	110.1 / NA	33	Q165Q (CAG→CAA)	cyuA	putative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/20);  total (13/20)

GGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTG  >  NC_000913/3255911‑3256385                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     ggATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTgg                                                                                                                                                                                                                                   <  2:145812/250‑1 (MQ=255)      aaGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAac                                                                                                                                                                                                                             >  1:239867/1‑251 (MQ=255)                  tCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATc                                                                                                                                                                                                                  <  1:51075/250‑1 (MQ=255)                     ccAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGc                                                                                                                                                                                                              <  1:192711/251‑1 (MQ=255)                     ccAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACGTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGGCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGc                                                                                                                                                                                                              <  1:14146/251‑1 (MQ=255)                                      cACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGTATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAAt                                                                                                                                                                                             <  1:216058/251‑1 (MQ=255)                                                                atGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTg                                                                                                                                                                    <  1:162298/250‑1 (MQ=255)                                                                               aCCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCAGCCGACGATGGTGACACACGCCCACTtc                                                                                                                                                     <  1:281628/250‑1 (MQ=255)                                                                                              cTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACtt                                                                                                                                     >  1:204949/1‑251 (MQ=255)                                                                                               tGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACttt                                                                                                                                    >  2:256387/1‑251 (MQ=255)                                                                                               tGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACttt                                                                                                                                    >  1:114962/1‑251 (MQ=255)                                                                                                  gATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGgc                                                                                                                                 >  1:113399/1‑251 (MQ=255)                                                                                                      aCGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGt                                                                                                                             >  1:298442/1‑251 (MQ=255)                                                                                                                    ttCGCGGAATCGAGAATAAAGCGGCTAGCCGCAAACGGGAATTCATTGACGAACTTCAGGATCTAAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGCGGATTTc                                                                                                               >  2:170577/1‑251 (MQ=255)                                                                                                                      cgcgGAATCGAGAATAAAGCGGAGCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCCGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCCCACACGCCTGTTGGGTAAACCCCACACCATCGGGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCAt                                                                                                             <  1:110016/251‑1 (MQ=255)                                                                                                                         ggAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGc                                                                                                          >  2:232478/1‑251 (MQ=255)                                                                                                                               gagaATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGtt                                                                                                    <  1:60780/251‑1 (MQ=255)                                                                                                                                                    aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACtt                                                                               >  2:89929/1‑251 (MQ=255)                                                                                                                                                     aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACttt                                                                              <  2:120341/251‑1 (MQ=255)                                                                                                                                                       cGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTcc                                                                            <  2:204949/251‑1 (MQ=255)                                                                                                                                                        gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTcc                                                                            >  2:179057/1‑250 (MQ=255)                                                                                                                                                        gggACTTCATGGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTcc                                                                            <  1:179057/250‑1 (MQ=255)                                                                                                                                                                 ttGACGAAATTCAGGATCTCAGCCAGAGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGAATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCgccgcc                                                                    >  1:208703/1‑249 (MQ=255)                                                                                                                                                                            cAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCCCACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTg                                                        <  1:96818/250‑1 (MQ=255)                                                                                                                                                                                    cagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGc                                                <  1:232478/250‑1 (MQ=255)                                                                                                                                                                                     agccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCCCACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCCCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCaa                                              <  1:256351/251‑1 (MQ=255)                                                                                                                                                                                           cgtcgtTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTGGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGGGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCt                                        <  1:118587/251‑1 (MQ=255)                                                                                                                                                                                                   tGGAAAGCACCGTCAGCGGAGACTCTTCCTCGCCCTCTGCCACACACGCCTGTTGGGCAAACACCACACCATCGTGCGTCTCGAAATGCACAACGTTGGTATGCCCGCCGACGATGGCGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGAGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGAGCTTTGGCATCGGCAATTGCCTGCGCTGtt                                >  2:305976/1‑251 (MQ=255)                                                                                                                                                                                                      aaaGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGcc                                                                                                                                                           >  2:103699/1‑125 (MQ=255)                                                                                                                                                                                                      aaaGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGcc                                                                                                                                                           <  1:103699/125‑1 (MQ=255)                                                                                                                                                                                                             cGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTTTTGCGTCTTTCa                      <  1:1031/251‑1 (MQ=255)                                                                                                                                                                                                                    ggAGCCTCTTGCTCCCCCTTGGCCCCACCCGCCTGTGGGGTAAACACCACCCCACCGTGCTTCCCGATCGGCACAATGTTGGAAGGCCCGCCGAGGATGGTGACACCCGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATGTCAGCGCAAGGTTCCTGGATCTTAACGTAGTCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGACTGCGCTGTTGCGTCTTTCAGCACTTc               <  2:37346/251‑1 (MQ=255)                                                                                                                                                                                                                                 cGCCCTCGCCCACCCACCCCTGTTGGGAAAACACCACCCCATCGGGGGTCTGGTTAGGCACATTTTTGGTATCCCCGCCGACGATGGTGACACCCGCCCCTTTCTCACCGTTCCAGACTTGGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTg  <  2:87931/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     GGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTG  >  NC_000913/3255911‑3256385

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