Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,943,707	T→C	V124A (GTC→GCC)	csdA →	cysteine sulfinate desulfinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,943,707	0	T	C	100.0%	141.4 / NA	43	V124A (GTC→GCC)	csdA	cysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/23);  total (20/23)

CCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAG  >  NC_000913/2943571‑2943840                                                                                                                                         |                                                                                                                                      ccGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGcccc                                                                                                                                    <  1:733086/140‑1 (MQ=255)       gatAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCt                                                                                                                              >  1:156434/1‑140 (MQ=255)       gatAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCt                                                                                                                              >  2:437031/1‑140 (MQ=255)              cTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTc                                                                                                                       <  1:134669/140‑1 (MQ=255)              cTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTc                                                                                                                       >  1:314757/1‑140 (MQ=255)               tATCGTCTGGACGCGCGGCACCCCTTAATCCATCAACGTGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCg                                                                                                                      <  1:263703/140‑1 (MQ=255)               tATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTACGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCg                                                                                                                      <  2:679260/140‑1 (MQ=255)                aTCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGc                                                                                                                     <  2:314757/140‑1 (MQ=255)                   gTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAATACCACGCTAACTTCTCCCCCTGGCTGATGGTCGCCca                                                                                                                  >  1:237251/1‑140 (MQ=255)                    tCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCcaa                                                                                                                 >  2:500639/1‑140 (MQ=255)                    tCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCcaa                                                                                                                 <  1:526843/140‑1 (MQ=255)                                  accacTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCaa                                                                                                   <  2:44302/140‑1 (MQ=255)                                              aTCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAAtt                                                                                       >  1:699828/1‑140 (MQ=255)                                              aTCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAAtt                                                                                       >  1:543465/1‑140 (MQ=255)                                              aTCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAAtt                                                                                       >  1:50123/1‑140 (MQ=255)                                               tCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTg                                                                                      >  2:89932/1‑140 (MQ=255)                                                       gtggCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTaa                                                                              <  2:543465/140‑1 (MQ=255)                                                        tggCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCTCAAAAAAATGGAGCCAAAGTGATGAAATTGGCGCTTAAt                                                                             >  2:337178/1‑140 (MQ=255)                                                        tggCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAAt                                                                             >  2:638589/1‑140 (MQ=255)                                                        tggCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAAt                                                                             >  1:642213/1‑140 (MQ=255)                                                           cacaATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATgcg                                                                          >  1:511915/1‑140 (MQ=255)                                                                 gCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGa                                                                    <  1:638589/140‑1 (MQ=255)                                                                 gCTATGCGCGTCCGCATCTGCATCCTGGCTATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGa                                                                    <  1:337178/140‑1 (MQ=255)                                                                     tGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGc                                                                >  1:364746/1‑140 (MQ=255)                                                                          gtccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGCGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGAt                                                           >  1:613886/1‑140 (MQ=255)                                                                           tccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATg                                                          >  1:98386/1‑140 (MQ=255)                                                                           tccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATg                                                          >  1:493501/1‑140 (MQ=255)                                                                               cgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACACCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGGCAAAGTGGGGAAATTGGCGCTTAATACGCAGGGACTGCCGGATGTCGa                                                      >  2:316426/1‑140 (MQ=255)                                                                                gtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGAt                                                     <  1:500639/140‑1 (MQ=255)                                                                                  cTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGAttt                                                   <  2:699828/140‑1 (MQ=255)                                                                                  cTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAAAACCCCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGAttt                                                   <  2:665195/140‑1 (MQ=255)                                                                                       aCCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGc                                              <  2:642213/140‑1 (MQ=255)                                                                                       aCCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGc                                              <  2:156434/140‑1 (MQ=255)                                                                                         cGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCa                                            <  2:364746/140‑1 (MQ=255)                                                                                         cGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCa                                            <  2:511915/140‑1 (MQ=255)                                                                                           ggCGATGAGATTATTGTTAGCGCGGCAGAACCCCAAGCCGACCTCGCCGCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGa                                          <  2:237251/140‑1 (MQ=255)                                                                                            gCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGaa                                         <  1:89932/140‑1 (MQ=255)                                                                                                    attattGTCAGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTAc                                 <  1:150089/140‑1 (MQ=255)                                                                                                             aGCGTGGCAGAACACCACGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGCATGTCGACTTGCTGCCAGAACTGATTACTCCCCGTAg                        >  1:310248/1‑140 (MQ=255)                                                                                                                         caccacGCCAACCTCGCGCCCTGGCTGATGGGCGCGCAAAAAACTGGATCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGc            <  2:310248/140‑1 (MQ=255)                                                                                                                          accacGCCTACCTCGCTCCCTGTCTGATGGCCGCCCAACGAACTGGAGGCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCg           <  2:357258/140‑1 (MQ=255)                                                                                                                          accacGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCg           <  2:493501/140‑1 (MQ=255)                                                                                                                          accacGCCAACCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCg           <  2:98386/140‑1 (MQ=255)                                                                                                                                   aCCTCGCCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAg  >  1:410085/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      CCGGATGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAG  >  NC_000913/2943571‑2943840

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