Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,657,963	G→T	intergenic (+93/‑106)	ynfD → / → ynfE	DUF1161 domain‑containing protein YnfD/putative selenate reductase YnfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,657,963	0	G	T	100.0%	166.5 / NA	47	intergenic (+93/‑106)	ynfD/ynfE	DUF1161 domain‑containing protein YnfD/putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (19/28);  total (19/28)

GCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATC  >  NC_000913/1657831‑1658098                                                                                                                                     |                                                                                                                                        gcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATtaat                                                                                                                                  <  2:600591/140‑1 (MQ=255)  cAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATtaata                                                                                                                                 <  2:702459/140‑1 (MQ=255)  cAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATtaata                                                                                                                                 <  1:504837/140‑1 (MQ=255)     tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAAtt                                                                                                                              >  1:650493/1‑140 (MQ=255)      cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTc                                                                                                                             <  1:665238/140‑1 (MQ=255)            aGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAAcc                                                                                                                       <  1:922797/140‑1 (MQ=255)               aTGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCt                                                                                                                            <  1:449012/132‑1 (MQ=255)               aTGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCt                                                                                                                            >  2:449012/1‑132 (MQ=255)               aTGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCaa                                                                                                                    >  2:706519/1‑140 (MQ=255)                tGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCaaa                                                                                                                   <  2:767185/140‑1 (MQ=255)                   cgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATg                                                                                                                >  1:98578/1‑140 (MQ=255)                     cgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATggg                                                                                                              <  1:889614/140‑1 (MQ=255)                            cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTc                                                                                                       <  1:914565/140‑1 (MQ=255)                            cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTc                                                                                                       <  2:706943/140‑1 (MQ=255)                                          acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAgc                                                                                         <  1:706519/140‑1 (MQ=255)                                          acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAgc                                                                                         <  1:720849/140‑1 (MQ=255)                                           caTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAgcg                                                                                        <  1:116327/140‑1 (MQ=255)                                            aTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAgcgc                                                                                       <  2:12580/140‑1 (MQ=255)                                              tATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCC‑CCTTATTACTCATACCCTTCTATTGATATTGATTAGTAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCt                                                                                    >  2:816821/1‑140 (MQ=255)                                               aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCt                                                                                    <  2:584142/140‑1 (MQ=255)                                                tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTa                                                                                   >  2:418420/1‑140 (MQ=255)                                                      gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATAtt                                                                             >  1:1025968/1‑140 (MQ=255)                                                      gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCC‑TCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTc                                                                            >  1:105040/1‑140 (MQ=255)                                                       aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTc                                                                            >  2:552146/1‑140 (MQ=255)                                                       aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATCCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTc                                                                            <  2:76807/140‑1 (MQ=255)                                                        aCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCa                                                                           <  1:434589/140‑1 (MQ=255)                                                             ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAatata                                                                      <  2:101556/140‑1 (MQ=255)                                                                        aTCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTtat                                                           >  2:285516/1‑140 (MQ=255)                                                                           cGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTtatata                                                        <  2:98578/140‑1 (MQ=255)                                                                              cGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACg                                                     <  1:774363/140‑1 (MQ=255)                                                                                  ttGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATaa                                                 <  2:157338/140‑1 (MQ=255)                                                                                  ttGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATaa                                                 <  2:1025968/140‑1 (MQ=255)                                                                                  ttGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATaa                                                 <  1:418420/140‑1 (MQ=255)                                                                                      aCATACAACATAATCCC‑CCTTATTACTCATACCCTTCTATTGATATTGATTAGTAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAAc                                            <  1:816821/140‑1 (MQ=255)                                                                                        aTACAACATAATCCCACCTTATTACTCATACCC‑TCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAg                                          <  2:105040/140‑1 (MQ=255)                                                                                           caacaTAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAggg                                        <  2:650493/140‑1 (MQ=255)                                                                                               aTAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGga                                       >  1:1071842/1‑137 (MQ=255)                                                                                               aTAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGga                                       <  2:1071842/137‑1 (MQ=255)                                                                                                tAATCCCACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGgagtga                                   >  1:1079558/1‑140 (MQ=255)                                                                                                      cACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGgag                                      <  1:167830/131‑1 (MQ=255)                                                                                                      cACCTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGgag                                      >  2:167830/1‑131 (MQ=255)                                                                                                        ccTTATTACTCATACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTc                           >  1:126566/1‑140 (MQ=255)                                                                                                                   aTACCCTTCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAAc                >  1:812996/1‑140 (MQ=255)                                                                                                                         ttCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAAtgg          >  2:746603/1‑140 (MQ=255)                                                                                                                          tCTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAAtggt         >  1:1022/1‑140 (MQ=255)                                                                                                                           cTATTGATATTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAATTTTTATATAACGTTAAAGAGCAGGGAGTAAGTTATGTCCAAAAATGAACGATaggtg        >  2:569019/1‑140 (MQ=255)                                                                                                                             attgatattgatTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTggg      >  1:1091918/1‑140 (MQ=255)                                                                                                                                 atatTGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATc  >  1:19026/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                        GCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATC  >  NC_000913/1657831‑1658098

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