BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,177,274	G→A	R3C (CGT→TGT)	yqiB ←	DUF1249 protein YqiB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,274	0	G	A	100.0%	125.7 / NA	37	R3C (CGT→TGT)	yqiB	DUF1249 protein YqiB
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (13/24); total (13/24)

CCGATATTGTGCGTTTGCCACCTGATAGCTTACAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACGCTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTA  >  NC_000913/3177143‑3177399

ccGATATTGTGCGTTTGCCACCTGATAGCTTACAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCAtt                                                                                                                        <  1:453615/139‑1 (MQ=255)

    tatTGTGCGTTTGCCACCTGATAGCTTACAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATg                                                                                                                    <  2:97895/139‑1 (MQ=255)

     atTGTGCGTTTGCCACCTGATAGCTTACAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGc                                                                                                                   <  1:778722/139‑1 (MQ=255)

           cGTTTGCCACCTGATAGCTTACAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTc                                                                                                             >  2:72335/1‑139 (MQ=255)

                    ccTGATAGCTTACAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTaa                                                                                                    <  1:815356/139‑1 (MQ=255)

                     cTGATAGCTTACAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAAc                                                                                                   <  2:939514/139‑1 (MQ=255)

                               aCAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTgatgat                                                                                         >  1:324272/1‑139 (MQ=255)

                               aCAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTgatgat                                                                                         >  1:220717/1‑139 (MQ=255)

                                cAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGatgatg                                                                                        <  2:757070/139‑1 (MQ=255)

                                     tCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGAtgcagct                                                                                   <  2:220717/139‑1 (MQ=255)

                                     tCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGAtgcagct                                                                                   >  1:921776/1‑139 (MQ=255)

                                      cGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGatgatg                                                                                        <  2:338942/133‑1 (MQ=255)

                                      cGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGatgatg                                                                                        >  1:338942/1‑133 (MQ=255)

                                            gTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCa                                                                            <  2:356475/139‑1 (MQ=255)

                                                      gcgcGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCg                                                                  <  2:604876/139‑1 (MQ=255)

                                                       cgcgGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGcc                                                                                                                  <  1:791797/90‑1 (MQ=255)

                                                       cgcgGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGcc                                                                                                                  >  2:791797/1‑90 (MQ=255)

                                                                      cgcAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACg                                                                                                  <  2:243189/91‑1 (MQ=255)

                                                                      cgcAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACg                                                                                                  >  1:243189/1‑91 (MQ=255)

                                                                      cgcAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGt                                                  >  1:869703/1‑139 (MQ=255)

                                                                      cgcAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGt                                                  >  2:17992/1‑139 (MQ=255)

                                                                              tgAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGAtt                                          <  1:1099269/139‑1 (MQ=255)

                                                                                   aGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTccc                                     <  2:181288/139‑1 (MQ=255)

                                                                                     ttCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACt                                                                                                              <  1:606815/64‑1 (MQ=255)

                                                                                     ttCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACt                                                                                                              >  2:606815/1‑64 (MQ=255)

                                                                                     ttCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCct                                   <  1:569757/139‑1 (MQ=255)

                                                                                       cATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCctct                                 >  1:868432/1‑139 (MQ=255)

                                                                                          ctcGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCtcttct                              <  1:854194/139‑1 (MQ=255)

                                                                                          ctcGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCtcttct                              <  1:72335/139‑1 (MQ=255)

                                                                                          ctcGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCtcttct                              <  1:538774/139‑1 (MQ=255)

                                                                                           tcGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTa                             <  1:417228/139‑1 (MQ=255)

                                                                                              cacaGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTAccc                          <  1:17992/139‑1 (MQ=255)

                                                                                               acaGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGAc                                                             <  2:63869/103‑1 (MQ=255)

                                                                                               acaGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGAc                                                             >  1:63869/1‑103 (MQ=255)

                                                                                                        tcatTTCAGGAAAGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGc                <  1:312834/139‑1 (MQ=255)

                                                                                                                    aGTCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCt    <  2:324272/139‑1 (MQ=255)

                                                                                                                      tCAGGTGTGTAACACTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTa  >  1:836702/1‑139 (MQ=255)

CCGATATTGTGCGTTTGCCACCTGATAGCTTACAGTTTCGCCGGGTGCGTCATTGCGCGGTAACAAACGGCGCAATTGTGAAAAGTTCATCTCGCACAGGCGCATCATTTCAGGAAAGTCAGGTGTGTAACGCTTCATTTATGCCCACTCATTTTTTAACGCTTGATGATGCAGCTGCAGCCATTGCAAAGCGATGACCGACGCTGCGTTGTCGATTTTCCCCTCTTCTACCCACTGGTATGCCTGTTCCCGGCTTA  >  NC_000913/3177143‑3177399

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