BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,622,811	A→G	intergenic (+70/+6)	mgtS → / ← mgrR	small protein MgtS/small regulatory RNA MgrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,622,811	0	A	G	100.0%	49.4 / NA	15	intergenic (+70/+6)	mgtS/mgrR	small protein MgtS/small regulatory RNA MgrR
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (13/2); total (13/2)

AATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGC‑GGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTT  >  NC_000913/1622687‑1622936

aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGgaaaaaaaac                                                                                                                 >  2:121039/1‑138 (MQ=255)

       ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGAAAAAAAAACGCCAGTa                                                                                                          >  2:630491/1‑140 (MQ=255)

       ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATG‑ACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGAAAAAAAAACGCCAGTaa                                                                                                         >  1:526278/1‑140 (MQ=255)

          tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGAAAAAAAAACCCCAGTaaaa                                                                                                       >  2:838696/1‑139 (MQ=255)

                      cccGCTTTTTCCACCCAAAAAGGGGATAACAAAGGAAGGGAAATACCCCCCCACCAACCCGCCCCTTTGTTACATTTGTGAACAGGGGGCT‑TTTTTTTTAGGGGGGAAAAAAAAAACCGCCAGTAAACCGGC‑GGTGAATg                                                                                          <  2:797254/139‑1 (MQ=255)

                                      acaAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGAAAAAAAACCGCCAGTAAACCGGC‑GGGGAATGCTTGCATGGGTAGAtt                                                                          >  1:116753/1‑140 (MQ=255)

                                                              agaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGCAAAAAAACCGCCAGTAAACCGGG‑GGTGAATTCTTTCATGGATAGGATTTTGTTTTTGTTTTTCGCTTACaa                                                  >  2:413750/1‑139 (MQ=255)

                                                                        tCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGAAAAAAAACCCCCAGTAAACCGGC‑GGGGGAATCTTGCGTGGATAGATTTGTGTTTTTTTTTTTCGCTAACAGGCATTTTTCt                                        >  2:82004/1‑140 (MQ=255)

                                                                            caccccGCCCCTTTTTTACATTTGTG‑AAAAGGGGCCTTTTTTTTTGGGGCGGAAAAAAAAACCCGCCAGTAAACCGGC‑GGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                    <  2:672816/136‑1 (MQ=255)

                                                                                    ccccTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGGAATGTTTGC‑TGGGTAGATTTTTGGTTTTCTTTTTCGCTAACAAGGATTTTTCTTCACTGATaaaa                            >  2:564194/1‑138 (MQ=255)

                                                                                    ccccTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGAAAAAAAACCGCCAGTAAACCGGC‑GGGGGATGCTTGCATGGATAGATTTTTGTTTTTTTTTTTCGCTAACAGGCATTTTTCTGCACTGATaaaa                            >  2:288075/1‑138 (MQ=255)

                                                                                      ccTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGGCGGCGAAAAAAAACCGCCAGTAAACCGGC‑GGCGGAATTTTGGATGGGGAAATTTTTGTTTTTCTTTTTAGCTAACAGGCAATTTTCTGCAATGATAACGaa                          >  2:792156/1‑140 (MQ=255)

                                                                                                   gtgaACAGGGGGCT‑TGTTTTTTAGGGGGGAAAAAAAAAACCGCCAGTAAACCGGC‑GGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt             <  1:17078/140‑1 (MQ=255)

                                                                                                      tACAAGGGGCC‑TGATTTTTATGGCGGCCAAAAAAAACCGCCACCAAACCCGC‑GGGGAAAGGCTTCCCGGGTAGGATTTTGTTTTTCTTTTTCCCCAAAAGGCCTTTTTCCCCCCTGCGAACCAAACGTTGGcacaccagg          >  2:194253/1‑135 (MQ=255)

                                                                                                              gCC‑TGATTTTTATGGCGGCGAAAAAAAACCGCCAGTAAACCGGC‑GGGGGAAGCTTTCATGGGTAGAATTGTGTTTTTCTTTTTCGCTAACAGGCATTTTTCTGCACTGATAAAGAATCGTTGACACAATAGCATCAGttt  >  2:17078/1‑140 (MQ=255)

                                                                                                              gCC‑TGATTTTTATGGCGGCGAAAAAAAACCGCCAGTAAACCGGC‑GGGGGAAGCTTGGATGGGTAGATTTGTGTTTTTCTTTTTCGCCAACAGGCATTTTTCTGCCCCGGTAACGAAATCTTTACACAGGAACATCAGgtt  >  2:292744/1‑140 (MQ=255)

AATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGC‑GGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTT  >  NC_000913/1622687‑1622936

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