Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,662,588	A→G	D678G (GAT→GGT)	ynfF →	putative selenate reductase YnfF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,662,588	0	A	G	100.0%	179.0 / NA	51	D678G (GAT→GGT)	ynfF	putative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/28);  total (23/28)

CCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCATTGATGCTCA  >  NC_000913/1662451‑1662723                                                                                                                                          |                                                                                                                                        cccGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTc                                                                                                                                       >  2:1424856/1‑140 (MQ=255)   cGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTccc                                                                                                                                     <  1:760672/140‑1 (MQ=255)      tGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTg                                                                                                                                  >  1:1178921/1‑140 (MQ=255)      tGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTg                                                                                                                                  >  1:963295/1‑140 (MQ=255)        aaaaCACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCg                                                                                                                                <  1:1279584/140‑1 (MQ=255)        aaaaCACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAA‑TGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCgg                                                                                                                               <  2:1039355/140‑1 (MQ=255)              ccTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAAc                                                                                                                          <  2:1178921/140‑1 (MQ=255)                  cGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTAt                                                                                                                      >  2:100635/1‑140 (MQ=255)                       aaaaTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACt                                                                                                                 >  2:208836/1‑140 (MQ=255)                        aaaTCGAAATTTATTATGAACGTCTAGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTg                                                                                                                <  1:725096/140‑1 (MQ=255)                         aaTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGaa                                                                                                                             <  1:1399138/126‑1 (MQ=255)                         aaTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGaa                                                                                                                             >  2:1399138/1‑126 (MQ=255)                               aaTTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTaa                                                                                                         <  2:74089/140‑1 (MQ=255)                                    aTTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGc                                                                                                    <  1:1424856/140‑1 (MQ=255)                                          aaCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCtt                                                                                                  <  2:1176830/136‑1 (MQ=255)                                          aaCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCtt                                                                                                  >  1:1176830/1‑136 (MQ=255)                                           aCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCAtt                                                                                             >  2:1244683/1‑140 (MQ=255)                                               cTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACaa                                                                                         >  2:1235002/1‑140 (MQ=255)                                               cTGGCGAAGATTGCAAATACCTTGGAATTGAAAAAAGCTGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACaa                                                                                         <  1:548742/140‑1 (MQ=255)                                                  gCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCGGTTAACGGGCTTCCATTACAAAgc                                                                                      >  1:980195/1‑140 (MQ=255)                                                       gATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCGGTTAACGGGCTTCCATTACAAAGCGCGTa                                                                                 <  2:980195/140‑1 (MQ=255)                                                       gATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTa                                                                                 <  1:100635/140‑1 (MQ=255)                                                            cAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGCCTACCATTAAAAAGCGCGTACCCAc                                                                            >  1:620262/1‑140 (MQ=255)                                                             aGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACt                                                                           <  1:1244683/140‑1 (MQ=255)                                                                  ccTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGc                                                                      <  2:179857/140‑1 (MQ=255)                                                                           tGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCACTACAAAGCGCATACCCACTCCAGCTACGGCAat                                                             >  2:672583/1‑140 (MQ=255)                                                                             aaaaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAatat                                                           <  2:1344423/140‑1 (MQ=255)                                                                               aaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                         <  2:267740/140‑1 (MQ=255)                                                                               aaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                         <  2:940489/140‑1 (MQ=255)                                                                               aaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                         <  2:963295/140‑1 (MQ=255)                                                                               aaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                         <  2:990888/140‑1 (MQ=255)                                                                               aaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                         <  1:1235002/140‑1 (MQ=255)                                                                               aaaaGATGAAATCATTCCTCCCCTTCCTGCGCATTTCCCAGGTTTTGAGGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                         <  2:517686/140‑1 (MQ=255)                                                                                     tGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATgtgt                                                   <  1:1522391/140‑1 (MQ=255)                                                                                          ttattCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTAcag                                              <  1:1416777/140‑1 (MQ=255)                                                                                            attCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTAcagca                                            >  2:179476/1‑140 (MQ=255)                                                                                                  ccccTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTg                                      >  1:844589/1‑140 (MQ=255)                                                                                                       tCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTg                                      <  1:120027/135‑1 (MQ=255)                                                                                                       tCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTg                                      >  2:120027/1‑135 (MQ=255)                                                                                                       tCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCcac                                 >  2:266240/1‑140 (MQ=255)                                                                                                       tCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCcac                                 >  2:882183/1‑140 (MQ=255)                                                                                                       tCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCcac                                 >  2:1211394/1‑140 (MQ=255)                                                                                                         cTGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACaa                               <  2:112970/140‑1 (MQ=255)                                                                                                          tGCGTATACCCCAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAag                              <  1:208836/140‑1 (MQ=255)                                                                                                                    ccAGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGAt                    >  2:801064/1‑140 (MQ=255)                                                                                                                      aGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCa                  >  1:810649/1‑140 (MQ=255)                                                                                                                      aGGTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCa                  >  2:90172/1‑140 (MQ=255)                                                                                                                        gTTTTGATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAAc                <  2:1136949/140‑1 (MQ=255)                                                                                                                             gATGGCTGGGACGGTCCCCTGCGGAAAGCCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCAt           <  1:987702/140‑1 (MQ=255)                                                                                                                             gATGGCTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCAt           <  1:266240/140‑1 (MQ=255)                                                                                                                                  cTGGGACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAAGCCTGCCCACATGATGTGTGGCTCAACCCCATTGATg      >  1:654203/1‑140 (MQ=255)                                                                                                                                      gACGGTCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCATTGATGCTCa  >  2:96096/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        CCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCATTGATGCTCA  >  NC_000913/1662451‑1662723

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