BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,634,527	C→T	E262K (GAG→AAG)	rob ←	DNA‑binding transcriptional dual regulator Rob

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,634,527	0	C	T	100.0%	113.4 / NA	31	E262K (GAG→AAG)	rob	DNA‑binding transcriptional dual regulator Rob
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (13/18); total (13/18)

ACAGCGTAGATTAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTCAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTCCGTATACC  >  NC_000913/4634468‑4634585

caggCGTAGATTAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAatat                                                       <  1:1812273/62‑1 (MQ=255)

acaGCGTAGATTAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAatat                                                       <  1:1619227/65‑1 (MQ=255)

acaGCGTAGATTAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAatat                                                       >  1:1638821/1‑65 (MQ=255)

  aGCGTAGATTAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATcc                                                     >  1:3225528/1‑65 (MQ=255)

  aGCGTAGATTAACTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATcc                                                     <  1:413195/65‑1 (MQ=255)

   gCGTAGATTAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCt                                                    >  1:1489905/1‑65 (MQ=255)

       aGATTAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGAcc                                                <  1:2411403/65‑1 (MQ=255)

          ttAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCttt                                             >  1:1635954/1‑65 (MQ=255)

          ttAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCttt                                             >  1:215294/1‑65 (MQ=255)

          ttAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCttt                                             >  1:1878557/1‑65 (MQ=255)

               tGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGgc                                        >  1:1165703/1‑65 (MQ=255)

                    gATCTCCCGCTTTGGCATCTTCTGCCGGGCAGTATCGCTTAATATCCTGACCTTTACGGcgcgcc                                   >  1:332702/1‑63 (MQ=255)

                      tctcCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAg                                 <  1:1720784/65‑1 (MQ=255)

                      tctcCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAg                                 <  1:2817957/65‑1 (MQ=255)

                            gCTTTGGCATCTTCTGCCGGGTAGTATCCGTTAATATCCTGACCTTTACGGCGCGTCCGGTTGAg                           >  1:2459942/1‑65 (MQ=255)

                             cTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGc                          <  1:1947572/65‑1 (MQ=255)

                             cTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGc                          <  1:3286803/65‑1 (MQ=255)

                             cTTTGGCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGc                          <  1:424032/65‑1 (MQ=255)

                              tttGGCATCTTCTGCCTGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCCGGTTGAGCa                         >  1:6893/1‑65 (MQ=255)

                                  gCATCTTCTGCCGGGTCGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTgg                     <  1:2956180/65‑1 (MQ=255)

                                  gCATCTTCTGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTgg                     <  1:1855040/65‑1 (MQ=255)

                                      cttctGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATg                 <  1:1528692/65‑1 (MQ=255)

                                      cttctGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATg                 >  1:151234/1‑65 (MQ=255)

                                        tctGCCGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCa               >  1:374127/1‑65 (MQ=255)

                                             cGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTc          <  1:2356157/65‑1 (MQ=255)

                                             cGGGTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTc          <  1:932537/65‑1 (MQ=255)

                                              gggTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTcc         <  1:1054184/65‑1 (MQ=255)

                                              gggTAGTATCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTcc         <  1:2596287/65‑1 (MQ=255)

                                                gtagtaTCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTCCGt       <  1:2344834/65‑1 (MQ=255)

                                                gtagtaTCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTCCGt       >  1:2302116/1‑65 (MQ=255)

                                                gtagtaTCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTCCGt       <  1:1163084/65‑1 (MQ=255)

                                                     aTCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTCCGTATAcc  >  1:3052953/1‑65 (MQ=255)

                                                     aTCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTCCGTATAcc  >  1:276940/1‑65 (MQ=255)

                                                     aTCGCTTAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTCCGTATAcc  >  1:1260550/1‑65 (MQ=255)

ACAGCGTAGATTAATTGGGCGATCTCCCGCTTTGGCATCTTCTGCCGGGTAGTATCGCTCAATATCCTGACCTTTACGGCGCGTCAGGTTGAGCATTGGCATGCACGTTCCGTATACC  >  NC_000913/4634468‑4634585

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 23 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`