BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,936,381	C→T	T35M (ACG→ATG)	rbsB →	D‑ribose ABC transporter periplasmic binding protein, ribose chemotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,936,381	0	C	T	100.0%	127.0 / NA	37	T35M (ACG→ATG)	rbsB	D‑ribose ABC transporter periplasmic binding protein, ribose chemotaxis receptor
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (19/18); total (19/18)

AATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTA  >  NC_000913/3936252‑3936511

aataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaa                                                                                                                           <  2:211173/139‑1 (MQ=255)

 ataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                                                          <  2:171788/139‑1 (MQ=255)

 ataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCAACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                                                          <  2:434064/139‑1 (MQ=255)

   aaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAAccc                                                                                                                        <  2:66986/139‑1 (MQ=255)

         tACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCtt                                                                                                                  >  2:561342/1‑139 (MQ=255)

         tACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCtt                                                                                                                  >  1:605220/1‑139 (MQ=255)

           cAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTg                                                                                                                >  2:100503/1‑139 (MQ=255)

             ggACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTa                                                                                                              >  2:557523/1‑139 (MQ=255)

              gACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTAt                                                                                                             <  1:417523/139‑1 (MQ=255)

                 aTCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGc                                                                                                          >  2:240933/1‑139 (MQ=255)

                   cTTGAATATGAACATGAAAAACCTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTg                                                                                                        >  2:322787/1‑139 (MQ=255)

                           tGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATgg                                                                                                <  1:240933/139‑1 (MQ=255)

                                     aaaaCTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGaaa                                                                                      <  1:345766/139‑1 (MQ=255)

                                            gCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCgg                                                                               >  1:363040/1‑139 (MQ=255)

                                             cTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGa                                                                              >  2:707786/1‑139 (MQ=255)

                                              tACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGAt                                                                             <  1:373200/139‑1 (MQ=255)

                                                  cTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                         <  1:557523/139‑1 (MQ=255)

                                                        tCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCGGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCt                                                                   <  1:451726/139‑1 (MQ=255)

                                                         ccGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCta                                                                  <  1:322787/139‑1 (MQ=255)

                                                          cgttGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCtat                                                                 >  2:235062/4‑139 (MQ=255)

                                                          cGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCtat                                                                 >  1:578348/1‑139 (MQ=255)

                                                           gCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCtata                                                                >  2:107060/1‑139 (MQ=255)

                                                               ttGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTa                                                                                                                               >  2:677702/1‑72 (MQ=255)

                                                               ttGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTa                                                                                                                               <  1:677702/72‑1 (MQ=255)

                                                                          gcCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTgg                                                 <  1:707786/139‑1 (MQ=255)

                                                                            cACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGAc                                               >  1:379513/1‑139 (MQ=255)

                                                                                   aGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                         >  1:632130/1‑106 (MQ=255)

                                                                                   aGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                         <  2:632130/106‑1 (MQ=255)

                                                                                                tGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGaa                           <  2:394797/139‑1 (MQ=255)

                                                                                                    aaaaGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                         >  1:734540/1‑89 (MQ=255)

                                                                                                    aaaaGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                         <  2:734540/89‑1 (MQ=255)

                                                                                                      aaGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCt                     >  2:580262/1‑139 (MQ=255)

                                                                                                      aaGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCt                     <  2:605220/139‑1 (MQ=255)

                                                                                                      aaGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCt                     >  2:189516/1‑139 (MQ=255)

                                                                                                                cgcgCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTg           >  2:473288/1‑139 (MQ=255)

                                                                                                                cgcgCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTg           <  1:189516/139‑1 (MQ=255)

                                                                                                                         ggtCTCCATGCATAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGCAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTa  >  2:83963/1‑139 (MQ=255)

AATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTA  >  NC_000913/3936252‑3936511

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