BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,472,408	A→G	intergenic (+117/‑28)	b4251 → / → b4252	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,472,408	0	A	G	100.0%	302.2 / NA	110	intergenic (+117/‑28)	b4251/b4252	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (64/46); total (64/46)

Only 100 of 110 total aligned reads displayed.

TGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCC  >  NC_000913/4472336‑4472478

tgtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgt                                                                     <  1:3547520/76‑1 (MQ=255)

tgtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgt                                                                     >  1:4350234/1‑76 (MQ=255)

tgtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgt                                                                     >  2:3748736/1‑76 (MQ=255)

 gtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgtt                                                                    <  2:3930332/76‑1 (MQ=255)

 gtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgtt                                                                    <  1:5347186/76‑1 (MQ=255)

   gATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                  >  1:830895/1‑76 (MQ=255)

   gATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                  <  2:598211/76‑1 (MQ=255)

     tAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGttt                                                                >  2:5589717/1‑76 (MQ=255)

     tAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGttt                                                                >  2:307587/1‑76 (MQ=255)

      aaGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                               <  1:3149696/76‑1 (MQ=255)

       aGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTtta                                                              <  2:4119596/76‑1 (MQ=255)

       aGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTtta                                                              <  2:4087341/76‑1 (MQ=255)

          cAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatt                                                           <  2:5833059/76‑1 (MQ=255)

          cAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatt                                                           >  2:3734557/1‑76 (MQ=255)

           aataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                          >  1:827273/1‑76 (MQ=255)

            ataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaattaa                                                         >  1:3453108/1‑76 (MQ=255)

              aatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAgg                                                       >  1:3337204/1‑76 (MQ=255)

               atTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                      >  1:147478/1‑76 (MQ=255)

                 tCACGCGTTAAGGATAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGt                                                    <  2:765723/76‑1 (MQ=255)

                  cACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTg                                                   >  2:4429677/1‑76 (MQ=255)

                  cACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTg                                                   >  1:4163614/1‑76 (MQ=255)

                  cACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTg                                                   >  1:45016/1‑76 (MQ=255)

                  cACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTg                                                   <  2:1254027/76‑1 (MQ=255)

                   aCGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTgg                                                  <  1:3476091/76‑1 (MQ=255)

                   aCGCGTGAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTgg                                                  <  2:152430/76‑1 (MQ=255)

                     gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTCAGGAGTggca                                                >  2:2326399/1‑74 (MQ=255)

                     gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaa                                                >  1:5569371/1‑76 (MQ=255)

                      cgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaa                                               >  1:875385/1‑76 (MQ=255)

                      cgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaa                                               >  2:1353606/1‑76 (MQ=255)

                       gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                              >  1:1682856/1‑76 (MQ=255)

                       gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                              >  1:5687101/1‑76 (MQ=255)

                       gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                              <  2:4706429/76‑1 (MQ=255)

                       gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                              >  1:3947984/1‑76 (MQ=255)

                       gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                              <  1:2996512/76‑1 (MQ=255)

                        ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAg                                             >  2:2624557/1‑76 (MQ=255)

                        ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAg                                             <  1:5439657/76‑1 (MQ=255)

                         tAAGGCTAGCGGAATTGATTATCTTTTCTTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGCGTGGAACAga                                            <  1:1703564/76‑1 (MQ=255)

                         tAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAga                                            >  2:3763880/1‑76 (MQ=255)

                          aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgat                                           <  2:5407453/76‑1 (MQ=255)

                          aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgat                                           <  1:1567162/76‑1 (MQ=255)

                          aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgat                                           >  1:2506870/1‑76 (MQ=255)

                           aGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatg                                          >  1:2156644/1‑76 (MQ=255)

                            ggCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatga                                         >  1:4932510/1‑76 (MQ=255)

                             gCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatgat                                        >  1:5119892/1‑76 (MQ=255)

                             gCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatgat                                        <  1:4031509/76‑1 (MQ=255)

                               tAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatca                                      >  1:746518/1‑76 (MQ=255)

                                aGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcat                                     >  2:5465694/1‑76 (MQ=255)

                                aGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcat                                     >  2:2819657/1‑76 (MQ=255)

                                 gCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcatc                                    >  1:5576316/1‑76 (MQ=255)

                                   ggAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCgg                                  >  1:5168012/1‑76 (MQ=255)

                                   ggAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCgg                                  <  2:1995997/76‑1 (MQ=255)

                                    gAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGa                                 <  2:1343780/76‑1 (MQ=255)

                                      aTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaat                               >  1:5764046/1‑75 (MQ=255)

                                      aTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                               >  1:4027071/1‑76 (MQ=255)

                                       ttGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAat                              <  2:4932509/76‑1 (MQ=255)

                                       ttGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAat                              <  2:827273/76‑1 (MQ=255)

                                         gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                            <  1:719869/76‑1 (MQ=255)

                                         gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                            >  2:2456692/1‑76 (MQ=255)

                                         gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                            <  1:4196546/76‑1 (MQ=255)

                                         gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                            <  2:1441631/76‑1 (MQ=255)

                                          aTTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTCATTAAGGAGTGTAAAAGATGATCATCGGAAatagt                           >  1:3158516/1‑74 (MQ=255)

                                          aTTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATAtt                           >  1:3452149/1‑76 (MQ=255)

                                           ttATCTTTTCGTCTAACGATAGAAATGACGCGTTGTTTTAATTAATGAGTGGAAAAGATTATCATCGGAAATATTc                          <  1:1310256/76‑1 (MQ=255)

                                           ttATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTc                          >  2:4992077/1‑76 (MQ=255)

                                           ttATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTc                          >  2:4934578/1‑76 (MQ=255)

                                            tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                         <  1:1645191/76‑1 (MQ=255)

                                             aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCAt                        >  1:5067265/1‑76 (MQ=255)

                                             aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCAt                        <  1:2631018/76‑1 (MQ=255)

                                              tCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATa                       <  2:786451/76‑1 (MQ=255)

                                              tCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATa                       >  1:4613356/1‑76 (MQ=255)

                                              tCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATa                       <  1:2287526/76‑1 (MQ=255)

                                              tCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATa                       >  1:2708345/1‑76 (MQ=255)

                                               cTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTGAGGAGTGGAAAAGATGATCATCGGAAATATTCATaa                      >  1:2603522/1‑76 (MQ=255)

                                                ttttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                     <  1:1329175/76‑1 (MQ=255)

                                                 tttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATc                    <  1:3542040/76‑1 (MQ=255)

                                                  ttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCt                   >  2:635588/1‑76 (MQ=255)

                                                  ttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCt                   >  1:2474074/1‑76 (MQ=255)

                                                  ttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCt                   >  1:3285951/1‑76 (MQ=255)

                                                   tCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAATGATGATCATCGGAAATATTCATAATCtt                  >  1:818779/1‑76 (MQ=255)

                                                    cGTATAGCGATAGAAATGAAGCCTTGTTTTAATTAAGGAGTGGAAAAGATCATCATCGGACGTATTCATAATCTTc                 <  2:584711/76‑1 (MQ=255)

                                                    cGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTc                 >  1:1769770/1‑76 (MQ=255)

                                                     gTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCa                <  2:643546/76‑1 (MQ=255)

                                                     gTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCa                <  2:2537114/76‑1 (MQ=255)

                                                      tataACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAGAAGATGATCATCGGAAATATTCATAATCTTCAg               <  1:5550532/76‑1 (MQ=255)

                                                       ataACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGc              <  1:4320949/76‑1 (MQ=255)

                                                       ataACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTAAGc              >  2:5656217/1‑76 (MQ=255)

                                                        taACGATAGAAAGGAAGCGTTGTTTTAATTAAGGAGAGGAAAAGATGATCATCGGAAATATTAATAATCTGCAGcc             >  1:334767/1‑76 (MQ=255)

                                                         aaCGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCg            <  2:2612623/76‑1 (MQ=255)

                                                          aCGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGt           <  2:3453108/76‑1 (MQ=255)

                                                             aTAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc        >  1:1862543/1‑76 (MQ=255)

                                                             aTAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc        >  1:4317660/1‑76 (MQ=255)

                                                              tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATTATCATCGGAGATATTCATAATCTTCAGCCGTGGCt       <  1:251540/76‑1 (MQ=255)

                                                              tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt       <  1:2630875/76‑1 (MQ=255)

                                                                gAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAc     >  2:3023135/1‑76 (MQ=255)

                                                                gAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAc     <  2:5295581/76‑1 (MQ=255)

                                                                 aaaTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcc    >  1:1632407/1‑76 (MQ=255)

                                                                 aaaTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcc    >  2:1670954/1‑76 (MQ=255)

                                                                   aTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcccc  >  1:1620968/1‑76 (MQ=255)

                                                                   aTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcccc  >  1:3677828/1‑76 (MQ=255)

                                                                   aTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcccc  >  1:666836/1‑76 (MQ=255)

TGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCC  >  NC_000913/4472336‑4472478

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