BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,538,270	A→G	F21S (TTT→TCT)	narV ←	nitrate reductase 2 (NRZ), gamma subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,538,270	0	A	G	100.0%	232.4 / NA	61	F21S (TTT→TCT)	narV	nitrate reductase 2 (NRZ), gamma subunit
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (0/61); total (0/61)

GGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTC  >  W3110S.gb/1538227‑1538317

ggggtACTGCCCTTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:669512/62‑1 (MQ=255)

ggTGTACTGCCCGTAGTCGTAGCGCAGCTAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:973916/65‑1 (MQ=38)

ggTGTACTGCCCGTAGTCGTAGCGCAGCTAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:970328/65‑1 (MQ=38)

ggTGTACTGCCCGTAGTCGTAGCGCAGCCTGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:81388/65‑1 (MQ=38)

ggTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:1178292/65‑1 (MQ=255)

ggTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:2780142/65‑1 (MQ=255)

ggTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:543530/65‑1 (MQ=255)

 gtgtACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:1861043/64‑1 (MQ=255)

  tgtACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:919267/63‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:719910/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:2698539/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:275731/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:2101489/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:2856830/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:189751/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:1829442/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:815526/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:826803/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:827162/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:852289/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:920165/61‑1 (MQ=255)

    tACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:1197984/61‑1 (MQ=255)

             tAGTCGTAGGGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1212505/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGtc               <  1:282897/65‑2 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:2592331/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:2275542/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:943081/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:891980/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:3005816/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:880397/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:635767/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:639670/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:655884/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1927482/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1393005/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1413128/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1434937/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1255339/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1245892/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1794997/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:116089/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:2041890/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:2185548/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1010997/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:2386486/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1062361/65‑1 (MQ=255)

             tAGTCGTAGCGCAGCCACCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTa               <  1:1446334/65‑1 (MQ=255)

                       gcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACg     <  1:1961069/65‑1 (MQ=255)

                       gcAGCCAGCTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACg     <  1:1434375/65‑1 (MQ=255)

                         agccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGtt   <  1:1079138/65‑1 (MQ=255)

                         agccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGtt   <  1:2072660/65‑1 (MQ=255)

                         agccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGtt   <  1:1482368/65‑1 (MQ=255)

                          gccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTc  <  1:685468/65‑1 (MQ=255)

                          gccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTc  <  1:1070567/65‑1 (MQ=255)

                          gccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTc  <  1:1375894/65‑1 (MQ=255)

                          gccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTc  <  1:1654682/65‑1 (MQ=255)

                          gccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTc  <  1:1878266/65‑1 (MQ=255)

                          gccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTc  <  1:351232/65‑1 (MQ=255)

                          gccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTc  <  1:221268/65‑1 (MQ=255)

                          gccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTc  <  1:2435399/65‑1 (MQ=255)

                           ccagcTGCCGAGGAAAGACACCGTCGCACAAATGTACg                            <  1:2518214/38‑1 (MQ=255)

GGTGTACTGCCCGTAGTCGTAGCGCAGCCAGCTGCCGAGGAAAAACACCGTCGCACAAATGTACGGGTAGATATCGTAAAAAAAGACGTTC  >  W3110S.gb/1538227‑1538317

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