BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,161,702	C→G	A162G (GCA→GGA)	btuB →	BtuB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,702	0	C	G	99.2%	317.7 / NA	118	A162G (GCA→GGA)	btuB	BtuB
Reads supporting (aligned to +/- strand): ref base C (0/0); major base G (64/53); minor base T (0/1); total (64/54)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 121 total aligned reads displayed.

ATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGG  >  NC_000913/4161630‑4161773

aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                      <  1:1827599/76‑1 (MQ=255)

aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                      <  1:3759769/76‑1 (MQ=255)

aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                      <  1:2958284/76‑1 (MQ=255)

  ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAgggg                                                                    >  1:2276677/1‑75 (MQ=255)

  ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAgggg                                                                    >  2:3585107/1‑75 (MQ=255)

  ggTGCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGt                                                                    <  2:663901/76‑1 (MQ=255)

   gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                   >  2:741181/1‑76 (MQ=255)

   gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                   >  2:5651579/1‑76 (MQ=255)

   gTTCCGATGCAATACGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                   >  1:3278235/1‑76 (MQ=255)

    ttCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgg                                                                  <  2:190144/76‑1 (MQ=255)

    ttCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgg                                                                  <  1:5379708/76‑1 (MQ=255)

     tCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGACCGGAAATTTCAGGAGGGTggg                                                                 <  2:2762020/76‑1 (MQ=255)

     tCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTggg                                                                 <  2:2910941/76‑1 (MQ=255)

     tCCGATGCAAGAGGCGGGGAGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATGTCAGGAGGGCggg                                                                 >  1:4277197/1‑76 (MQ=255)

       cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                               <  1:4978205/76‑1 (MQ=255)

       cGATGCAAGAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                               <  2:562418/76‑1 (MQ=255)

        gATGCACGAGGCGGGGTGTGGAATATCATCACGACGCGCGATCAAACCGGACCGGAAATTTCAGGAGGGTGGGGaa                                                              <  2:1655454/76‑1 (MQ=255)

        gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                              <  1:561014/76‑1 (MQ=255)

        gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                              <  1:4275282/76‑1 (MQ=255)

        gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                              >  2:4344541/1‑76 (MQ=255)

        gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                              <  2:4297884/76‑1 (MQ=255)

          tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGc                                                            >  1:2597359/1‑76 (MQ=255)

          tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGc                                                            >  1:2057322/1‑76 (MQ=255)

           gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGCAATTTCCGGAGGGTGGGGAAGCa                                                           >  1:5295599/1‑76 (MQ=255)

           gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                           >  2:3036772/1‑76 (MQ=255)

           gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                           <  2:2472073/76‑1 (MQ=255)

            cAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGAGGAAGCaa                                                          >  1:4920499/1‑76 (MQ=255)

             aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAAt                                                         >  2:4689779/1‑76 (MQ=255)

              aTAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGGGGGGGGGGGAGGaaaa                                                        >  2:509991/1‑74 (MQ=255)

              aTAGGCGGGGGGGGGAATATCATCACGACGCGCGCTGAACCCGGGACGGGAATTTTAGGGGGGGGGGGGAGGAATa                                                        >  2:4626618/1‑76 (MQ=255)

               tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                       <  1:233490/76‑1 (MQ=255)

               tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                       <  1:1718793/76‑1 (MQ=255)

                aGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                      >  1:2147659/1‑76 (MQ=255)

                aGGCGGGGGGGTGAATATCATCACGACGCGCGCTGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAATAGt                                                      >  2:3785317/1‑76 (MQ=255)

                 ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                     <  1:741181/76‑1 (MQ=255)

                  gCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                    <  2:901163/76‑1 (MQ=255)

                   cGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAt                                                   >  2:4166374/1‑76 (MQ=255)

                    ggggggGTGAATATCATCACGACGAGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGAGAAGCAATAGTTATc                                                  >  1:1944443/1‑76 (MQ=255)

                    ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                  <  2:1431228/76‑1 (MQ=255)

                    ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTAta                                                  >  2:5641047/1‑75 (MQ=255)

                     gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCa                                                 >  2:196013/1‑76 (MQ=255)

                     gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCa                                                 <  2:3163061/76‑1 (MQ=255)

                      ggtggtGTATATCATCACGACGCGCGCTGATCCCGGTGCGGAGATTTCGGGAGGGTGGGGACGCAATAGTTATCAg                                                <  2:3337103/76‑1 (MQ=255)

                       gtggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGCTATCAGa                                               >  2:686596/1‑76 (MQ=255)

                         ggtgAGTATCCTCACGACGCGCGGTGACCCCGGAAGGGAAATTGCAGGGGGGTGGGGAATCAATAGTTATCAGAAc                                             <  2:4718897/76‑1 (MQ=255)

                          gttaaTATCATCACGACGCGCGATGGACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                            <  1:4233288/73‑1 (MQ=255)

                          gtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                            <  1:4195744/76‑1 (MQ=255)

                           tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                           >  1:4215525/1‑76 (MQ=255)

                            gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                          <  1:1744706/76‑1 (MQ=255)

                              atatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                        >  2:3497187/1‑76 (MQ=255)

                              atatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                        <  1:3636872/76‑1 (MQ=255)

                              atatCATCAAGACGCGCGATGAACCCGAAAAGGAAATTTCAGGGGGGTGGGGGAGCAATCGTTATCCGAACTatga                                        >  1:3900910/1‑76 (MQ=255)

                               tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                       <  2:2730956/76‑1 (MQ=255)

                                atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                      >  1:3972010/1‑76 (MQ=255)

                                atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                      <  1:3479770/76‑1 (MQ=255)

                                atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                      <  2:802903/76‑1 (MQ=255)

                                 tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                     <  2:1617558/76‑1 (MQ=255)

                                 tcatcaCGACGCGCGATGAACCCGGAACGGAAATGTCAGGAGGGTGGGGAGGCATTTGTTATCAGAACTATGATGt                                     >  1:1765971/1‑76 (MQ=255)

                                  catcaCGACGCGCGATGACCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                    <  2:3136351/76‑1 (MQ=255)

                                  catcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                    >  2:2045554/1‑76 (MQ=255)

                                  catcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                    >  2:4457367/1‑76 (MQ=255)

                                   atcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtct                                   >  1:4669575/1‑76 (MQ=255)

                                    tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctc                                  <  2:3285272/76‑1 (MQ=255)

                                     caCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                 <  1:4539185/76‑1 (MQ=255)

                                      acgacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                >  1:3564710/1‑76 (MQ=255)

                                         acgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                             >  2:2117279/1‑76 (MQ=255)

                                           gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcag                           <  1:4626619/76‑1 (MQ=255)

                                            cgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                          >  2:1723781/1‑76 (MQ=255)

                                             gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                         >  1:839882/1‑76 (MQ=255)

                                              cgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                        >  1:4343600/1‑76 (MQ=255)

                                               gATGACCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaac                       <  1:3668452/76‑1 (MQ=255)

                                                aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                      <  2:1714421/76‑1 (MQ=255)

                                                  gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                    >  2:2100547/1‑76 (MQ=255)

                                                  gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                    >  1:5645443/1‑76 (MQ=255)

                                                      ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggg                <  2:2112909/76‑1 (MQ=255)

                                                        tgaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg              >  2:1378483/2‑76 (MQ=255)

                                                        ggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg              >  1:480373/1‑76 (MQ=255)

                                                         gaacggaaATTGCAGGAGGGTGGGGAAGCAATAGTTATCAGAAGTATGATGTCTCTACGCAGCCACAACTGGGGGa             >  1:3099394/1‑76 (MQ=255)

                                                          aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt            >  1:2714151/1‑76 (MQ=255)

                                                           acggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa           >  2:4770129/1‑76 (MQ=255)

                                                            cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACt                   >  2:2981445/1‑67 (MQ=255)

                                                            cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa          <  2:3001261/76‑1 (MQ=255)

                                                             ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAg         <  1:406955/76‑1 (MQ=255)

                                                             ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAg         >  1:3454810/1‑76 (MQ=255)

                                                             ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAg         >  2:4334465/1‑76 (MQ=255)

                                                             ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAg         >  1:2864343/1‑76 (MQ=255)

                                                               aaaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac       <  2:1903976/76‑1 (MQ=255)

                                                               aaaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac       >  2:5815530/1‑76 (MQ=255)

                                                               aaaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGACAAGac       >  2:445560/1‑76 (MQ=255)

                                                                aaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGaca      >  2:4009376/1‑76 (MQ=255)

                                                                 aTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac     <  1:4427323/76‑1 (MQ=255)

                                                                 aTTTCAGGAGGGTGGGGAAGCAATAGGTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac     >  2:1516747/1‑76 (MQ=255)

                                                                  tttCAGGAGGGAGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg    >  2:4736234/1‑76 (MQ=255)

                                                                   ttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg   >  1:1109227/1‑76 (MQ=255)

                                                                   ttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg   >  1:5608412/1‑76 (MQ=255)

                                                                   ttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg   >  1:4790590/1‑76 (MQ=255)

                                                                   ttCAGGAGGGTGGGGAAGCAATAGTCATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg   >  1:2941856/1‑76 (MQ=255)

                                                                    tCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg  >  1:3402569/1‑76 (MQ=255)

                                                                    tCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg  >  2:2927265/1‑76 (MQ=255)

                                                                    tCAGGAGGGTGGGAAAGCCATAGTTCACAGAACTATGATGTCTCTACGCAGCAGCAACTGGGGGCTCAGACACggg  <  2:3900909/76‑1 (MQ=255)

ATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGG  >  NC_000913/4161630‑4161773

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