BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,417,275	A→G	Y38H (TAC→CAC)	recE ←	Rac prophage; exonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,417,275	0	A	G	100.0%	148.8 / NA	43	Y38H (TAC→CAC)	recE	Rac prophage; exonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (23/20); total (23/20)

TCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCA  >  NC_000913/1417139‑1417406

tCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTc                                                                                                                                  <  2:523176/140‑1 (MQ=255)

  tttGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCa                                                                                                                                >  1:610064/1‑140 (MQ=255)

  tttGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCa                                                                                                                                >  1:806046/1‑140 (MQ=255)

        gAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTgg                                                                                                                          >  1:122880/1‑140 (MQ=255)

              aTAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAgg                                                                                                                    >  1:671893/1‑140 (MQ=255)

                    ttCACTCCAGGTAAGATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGAtt                                                                                                              >  2:294472/1‑140 (MQ=255)

                     tCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTc                                                                                                             >  1:623493/1‑140 (MQ=255)

                          ccAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAAt                                                                                                        <  2:115615/140‑1 (MQ=255)

                             ggTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGt                                                                                                     <  1:530252/140‑1 (MQ=255)

                                 aaaTCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGt                                                                                                                   >  2:637696/1‑122 (MQ=255)

                                 aaaTCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGt                                                                                                                   <  1:637696/122‑1 (MQ=255)

                                  aaTCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgctt                                                                                                <  1:454546/140‑1 (MQ=255)

                                  aaTCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgctt                                                                                                <  2:122880/140‑1 (MQ=255)

                                   aTCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgcttg                                                                                               >  2:462012/1‑140 (MQ=255)

                                   aTCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgcttg                                                                                               >  2:683451/1‑140 (MQ=255)

                                        tctcACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCcac                                                                                          <  1:381319/140‑1 (MQ=255)

                                                  aGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt                                                                                >  1:466893/1‑140 (MQ=255)

                                                  aGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt                                                                                <  2:651909/140‑1 (MQ=255)

                                                  aGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt                                                                                >  1:651909/1‑140 (MQ=255)

                                                         aGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTc                                                                         >  1:90325/1‑140 (MQ=255)

                                                          ggTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCa                                                                        <  1:179869/140‑1 (MQ=255)

                                                                tAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGgat                                                                  <  1:359738/140‑1 (MQ=255)

                                                                                cGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTcc                                                  >  1:580760/1‑140 (MQ=255)

                                                                                             ttAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAagag                                     <  2:610064/140‑1 (MQ=255)

                                                                                                    aGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGtttt                              <  2:580760/140‑1 (MQ=255)

                                                                                                     gctgctCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAagag                                     <  1:106220/132‑1 (MQ=255)

                                                                                                     gctgctCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAagag                                     >  2:106220/1‑132 (MQ=255)

                                                                                                     gctgctCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTtg                             >  1:416214/1‑140 (MQ=255)

                                                                                                     gctgctCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTtg                             >  2:68442/1‑140 (MQ=255)

                                                                                                      ctgctCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTtgt                            >  1:534481/1‑140 (MQ=255)

                                                                                                        gctCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGc                          <  2:806046/140‑1 (MQ=255)

                                                                                                         ctCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCt                         <  2:416214/140‑1 (MQ=255)

                                                                                                         ctCAGTTTTTTACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCt                         >  1:232630/1‑140 (MQ=255)

                                                                                                                   tACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGttttt               >  2:245641/1‑140 (MQ=255)

                                                                                                                   tACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGttttt               >  2:171961/1‑140 (MQ=255)

                                                                                                                   tACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGttttt               >  2:564055/1‑140 (MQ=255)

                                                                                                                   tACCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGttttt               >  1:267368/1‑140 (MQ=255)

                                                                                                                    aCCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGtttttt              <  1:462012/140‑1 (MQ=255)

                                                                                                                    aCCTGACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGtttttt              <  2:232630/140‑1 (MQ=255)

                                                                                                                        gACTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAAcc          >  1:293699/1‑140 (MQ=255)

                                                                                                                         aCTTAACGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCt         <  2:267368/140‑1 (MQ=255)

                                                                                                                              aCGATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGTAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACt    <  2:623493/140‑1 (MQ=255)

                                                                                                                                gATCAGGTGGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCa  <  2:90325/140‑1 (MQ=255)

TCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCA  >  NC_000913/1417139‑1417406

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