Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,660,581	T→C	L9P (CTC→CCC)	ynfF →	putative selenate reductase YnfF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,660,581	0	T	C	100.0%	98.5 / NA	28	L9P (CTC→CCC)	ynfF	putative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/12);  total (16/12)

GCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCG  >  NC_000913/1660445‑1660717                                                                                                                                         |                                                                                                                                         gCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCAt                                                                                                                                       >  2:586393/1‑140 (MQ=255)  cGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATg                                                                                                                                      <  1:693069/140‑1 (MQ=255)  cGAACGGAAAACCGCAGCACACTAAGCAGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTAGCATCACACTTTGGTCGAAAAGCAATAAGTGAGGAATGATGAAAACCCATCCCGCAGGGGCTCCCATg                                                                                                                                      <  2:844632/140‑1 (MQ=255)      gggAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAgg                                                                                                                                  >  2:519238/1‑140 (MQ=255)             ccGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAAt                                                                                                                           <  1:586393/140‑1 (MQ=255)                cagcaCACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAg                                                                                                                        <  2:421703/140‑1 (MQ=255)                 agcaCACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGc                                                                                                                       >  1:709576/1‑140 (MQ=255)                   cacaCTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCg                                                                                                                     <  1:1062774/140‑1 (MQ=255)                      acTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAg                                                                                                                  <  2:468118/140‑1 (MQ=255)                                  gagaTCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGaa                                                                                                      >  1:619888/1‑140 (MQ=255)                                    gaTCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGaaaa                                                                                                    >  1:846240/1‑140 (MQ=255)                                     aTCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAAc                                                                                                   >  1:329496/1‑140 (MQ=255)                                                        cGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCtt                                                                                >  2:433777/1‑140 (MQ=255)                                                           cAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTgcg                                                                             <  2:329496/140‑1 (MQ=255)                                                              ccATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTg                                                                          <  1:1037521/140‑1 (MQ=255)                                                               cATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTgg                                                                         >  1:996709/1‑140 (MQ=255)                                                                                                          gTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGggc                              >  1:715380/1‑140 (MQ=255)                                                                                                          gTAATGATGAAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGggc                              >  2:546379/1‑140 (MQ=255)                                                                                                                atgaAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTc                                   <  1:1085688/129‑1 (MQ=255)                                                                                                                atgaAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTc                                   >  2:1085688/1‑129 (MQ=255)                                                                                                                 tgaAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCggg                               <  1:1037963/132‑1 (MQ=255)                                                                                                                 tgaAAATCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCggg                               >  2:1037963/1‑132 (MQ=255)                                                                                                                   aaaaTCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCggg                               <  1:736428/130‑1 (MQ=255)                                                                                                                   aaaaTCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCggg                               >  2:736428/1‑130 (MQ=255)                                                                                                                     aaTCCATACCACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTc                   <  2:715380/140‑1 (MQ=255)                                                                                                                             ccACAGAGGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTagaa           >  1:216689/1‑140 (MQ=255)                                                                                                                                   agGCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGaaa     >  1:360796/1‑140 (MQ=255)                                                                                                                                    ggCGCCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAg    <  1:433777/140‑1 (MQ=255)                                                                                                                                      cgcCCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCg  >  2:1020856/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                         GCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCG  >  NC_000913/1660445‑1660717

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