Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,943,436	G→T	A34S (GCC→TCC)	csdA →	cysteine sulfinate desulfinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,943,436	0	G	T	100.0%	144.2 / NA	41	A34S (GCC→TCC)	csdA	cysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (18/23);  total (18/23)

TTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCA  >  NC_000913/2943299‑2943570                                                                                                                                          |                                                                                                                                       ttttGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATcc                                                                                                                                      <  2:338650/140‑1 (MQ=255) ttttGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATcc                                                                                                                                      <  2:338489/140‑1 (MQ=255)  tttGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCg                                                                                                                                     <  2:942333/140‑1 (MQ=255)       gTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGtt                                                                                                                                <  2:639969/140‑1 (MQ=255)           cGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAg                                                                                                                            <  2:907389/140‑1 (MQ=255)                 gTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCAccc                                                                                                                      <  2:182020/140‑1 (MQ=255)                   gCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCaa                                                                                                                    <  2:332249/140‑1 (MQ=255)                   gCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTACCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCaa                                                                                                                    <  1:109608/140‑1 (MQ=255)                     aTCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACa                                                                                                                  <  2:43521/140‑1 (MQ=255)                          gCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGtttt                                                                                                             >  1:349813/1‑140 (MQ=255)                            cGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCCCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTAc                                                                                                           >  2:79711/1‑140 (MQ=255)                              aggagTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAg                                                                                                         >  1:115428/1‑140 (MQ=255)                              aggagTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAg                                                                                                         >  2:23406/1‑140 (MQ=255)                               ggagTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGt                                                                                                        >  2:775329/1‑140 (MQ=255)                                  gTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTg                                                                                                     <  2:103554/140‑1 (MQ=255)                                     ccATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAgc                                                                                                  >  2:119441/1‑140 (MQ=255)                                        tGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGcc                                                                                               >  2:684950/1‑140 (MQ=255)                                                         gcgcAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAg                                                                              <  1:844677/140‑1 (MQ=255)                                                         gcgcAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAg                                                                              >  1:766497/1‑140 (MQ=255)                                                          cgcAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGc                                                                             >  2:912658/1‑140 (MQ=255)                                                                     gcCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGa                                                                  <  2:240383/140‑1 (MQ=255)                                                                       ccAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAg                                                                >  1:496454/1‑140 (MQ=255)                                                                         aGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCTACCGTGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGGTTTACTGTCTGAGCGCCGGAAACGTCCATCGCAGCCATTTTGCCGAAGcc                                                              <  2:236924/140‑1 (MQ=255)                                                                                        aGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACcgcg                                               <  2:62975/140‑1 (MQ=255)                                                                                        aGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACcgcg                                               <  2:756352/140‑1 (MQ=255)                                                                                        aGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACcgcg                                               <  1:955208/140‑1 (MQ=255)                                                                                               gggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                        >  2:893836/1‑140 (MQ=255)                                                                                               gggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                        >  1:608682/1‑140 (MQ=255)                                                                                               gggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                        >  2:100262/1‑140 (MQ=255)                                                                                               gggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTAAACGCCGGAAACGTCCATCGCAGGCTGTTTGCCGGAGCCTAACGCCTGATCGCGCGGTATg                                        >  2:950659/1‑140 (MQ=255)                                                                                               cggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATg                                        >  1:496925/2‑140 (MQ=255)                                                                                                ggCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGa                                       <  2:137892/140‑1 (MQ=255)                                                                                                      tATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGc                                 >  1:126439/1‑140 (MQ=255)                                                                                                       aTCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCa                                <  1:549547/140‑1 (MQ=255)                                                                                                       aTCTCGACAGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCa                                <  1:23406/140‑1 (MQ=255)                                                                                                               aGCGCCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaa                        <  1:684950/140‑1 (MQ=255)                                                                                                                gcgcCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaaa                       <  1:633381/140‑1 (MQ=255)                                                                                                                gcgcCGCGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGaaa                       >  2:890002/1‑140 (MQ=255)                                                                                                                 cgccgcGACCGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAg                      >  2:567303/1‑140 (MQ=255)                                                                                                                         ccGCGCTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGTGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCaa              <  2:115428/140‑1 (MQ=255)                                                                                                                              cTTAAACCTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACt         >  2:302156/1‑140 (MQ=255)                                                                                                                                    ccTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGc   <  1:119441/140‑1 (MQ=255)                                                                                                                                     cTGAATCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCa  <  1:518246/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       TTTTGCGTCACGATCGGTGCATCAAGCCGAGGAGTACCATGAACGTTTTTAATCCCGCGCAGTTTCGCGCCCAGTTTCCCGCACTACAGGATGCGGGCGTCTATCTCGACAGCGCCGCGACCGCGCTTAAACCTGAAGCCGTGGTTGAAGCCACCCAACAGTTTTACAGTCTGAGCGCCGGAAACGTCCATCGCAGCCAGTTTGCCGAAGCCCAACGCCTGACCGCGCGTTATGAAGCTGCACGAGAGAAAGTGGCGCAATTACTGAATGCA  >  NC_000913/2943299‑2943570

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