BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,961,588	C→G	R86T (AGG→ACG)	gppA ←	guanosine pentaphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,961,588	0	C	G	100.0%	421.2 / NA	158	R86T (AGG→ACG)	gppA	guanosine pentaphosphatase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (80/78); total (80/78)

Only 100 of 158 total aligned reads displayed.

CCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAT  >  NC_000913/3961517‑3961659

 cTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGcc                                                                    <  1:1363913/76‑1 (MQ=255)

 cTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGcc                                                                    >  1:3045304/1‑76 (MQ=255)

   gTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGt                                                                  <  1:1228373/76‑1 (MQ=255)

   gTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGt                                                                  >  1:1411992/1‑76 (MQ=255)

    tACCGGACAACCGGGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                 <  2:4719630/76‑1 (MQ=255)

    tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                 <  2:1592079/76‑1 (MQ=255)

     aCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGTCGTAACGTCGCAGTAg                                                                >  1:4354604/1‑76 (MQ=255)

      ccGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                               <  1:1080503/76‑1 (MQ=255)

      ccGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                               >  1:4667707/1‑76 (MQ=255)

       cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCg                                                              <  2:374027/76‑1 (MQ=255)

          aCAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACa                                                           <  2:3769980/76‑1 (MQ=255)

           cAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                          >  1:577405/1‑76 (MQ=255)

           cAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                          >  2:2965920/1‑76 (MQ=255)

             aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcg                                                        >  1:3951807/1‑76 (MQ=255)

             aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcg                                                        <  1:528083/76‑1 (MQ=255)

             aCCGAGAATGTCCTGCGCTTTGGCAATAAAATCACCCGCATTGAAGGCAACGCGTAACGTCGCCGTAGCGACAAcg                                                        >  1:745178/1‑76 (MQ=255)

              ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                       <  1:2615654/76‑1 (MQ=255)

               cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgcg                                                      >  2:1812571/1‑76 (MQ=255)

               cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgcg                                                      >  1:3140684/1‑76 (MQ=255)

                gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                     >  1:2275008/1‑76 (MQ=255)

                gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                     >  1:1936876/1‑76 (MQ=255)

                gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                     <  2:773039/76‑1 (MQ=255)

                gagaATTTCATGCGCTTTGGCAATAAAATCACACGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                     >  2:3596283/1‑76 (MQ=255)

                 agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                    >  1:2819840/1‑76 (MQ=255)

                  ggaTTTCCTGCGGTTTGGCAATAAAATCACCCGCATGGGCGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAt                                                   <  2:801665/74‑1 (MQ=255)

                  gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAt                                                   <  2:718309/76‑1 (MQ=255)

                   aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                  >  1:2328979/1‑76 (MQ=255)

                   aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACACGAAtt                                                  >  2:418394/1‑76 (MQ=255)

                    aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAttt                                                 <  1:1526019/76‑1 (MQ=255)

                    aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAttt                                                 >  1:4008245/1‑76 (MQ=255)

                    aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAttt                                                 <  2:4447352/76‑1 (MQ=255)

                    aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAttt                                                 <  2:2472734/76‑1 (MQ=255)

                     tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTg                                                >  2:2956161/1‑76 (MQ=255)

                      ttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGc                                               <  2:1154842/76‑1 (MQ=255)

                      ttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGc                                               <  2:3679140/76‑1 (MQ=255)

                       tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACACAGCGAATTTGCg                                              >  2:4380282/1‑76 (MQ=255)

                       tCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCg                                              >  2:4341543/1‑76 (MQ=255)

                        ccTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCga                                             <  2:4020423/76‑1 (MQ=255)

                         cTGCGCTTTGGCACTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgag                                            <  2:3888036/76‑1 (MQ=255)

                          tGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagg                                           <  2:2602536/76‑1 (MQ=255)

                            cgcTGTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGCCAACGCGAATTTGCgggggg                                         >  2:266171/1‑75 (MQ=255)

                             gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggag                                        >  1:3987405/1‑76 (MQ=255)

                               tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGggagga                                      <  2:767728/76‑1 (MQ=255)

                               tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGggagga                                      >  1:470701/1‑76 (MQ=255)

                               tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGggagga                                      >  1:4262950/1‑76 (MQ=255)

                                 tGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAat                                    >  2:3716622/1‑76 (MQ=255)

                                  ggCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAata                                   <  1:956311/76‑1 (MQ=255)

                                   gCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAatat                                  >  1:3850616/1‑76 (MQ=255)

                                    cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATc                                 >  2:1674237/1‑76 (MQ=255)

                                    cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATc                                 >  1:620210/1‑76 (MQ=255)

                                    cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATc                                 >  1:4818378/1‑76 (MQ=255)

                                     aaTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCAACAACGCGAATTTGCGAGGGAGGAATATCt                                <  1:2318334/76‑1 (MQ=255)

                                     aaTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCAACAACGCGAATTTGCGAGGGAGGAATATCt                                <  1:2285916/76‑1 (MQ=255)

                                      aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                               <  1:3407407/76‑1 (MQ=255)

                                      aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                               <  2:1973160/76‑1 (MQ=255)

                                      aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGGAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                               <  1:4533074/76‑1 (MQ=255)

                                        aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGc                             >  2:4827606/1‑76 (MQ=255)

                                        aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGc                             >  1:1951650/1‑76 (MQ=255)

                                          aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAg                           >  2:3325335/1‑76 (MQ=255)

                                          aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAg                           >  1:2776121/1‑76 (MQ=255)

                                          aaTCACCCGCATTGACGGCAAGGCGCAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAg                           >  1:3318115/1‑76 (MQ=255)

                                           aTCCCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGa                          <  2:1858854/76‑1 (MQ=255)

                                            gcACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGAc                         <  2:1165167/75‑1 (MQ=255)

                                             cACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACg                        <  2:3606801/76‑1 (MQ=255)

                                                ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                     <  2:2545280/76‑1 (MQ=255)

                                                ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                     <  1:830652/76‑1 (MQ=255)

                                                 cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                    >  1:1683755/1‑76 (MQ=255)

                                                  gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAGTTGCGAGGGAGGAATATCTTGCAGACGCTCAg                   >  1:2616519/1‑76 (MQ=255)

                                                   cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGc                  <  1:3062036/76‑1 (MQ=255)

                                                   cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGc                  >  1:3336659/1‑76 (MQ=255)

                                                   cATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGc                  >  1:4120120/1‑76 (MQ=255)

                                                    aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                 <  1:957676/76‑1 (MQ=255)

                                                    aTTGACGGAAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGGTCAGCa                 >  1:1190303/1‑76 (MQ=255)

                                                     ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                <  2:4120120/76‑1 (MQ=255)

                                                     ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                <  2:4441661/76‑1 (MQ=255)

                                                     ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                <  2:866032/76‑1 (MQ=255)

                                                     ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaa                <  1:4645141/76‑1 (MQ=255)

                                                      tGACGGCGAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaa               <  2:4763358/76‑1 (MQ=255)

                                                      gggcGGCAAGGGGTAGCGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaaa               <  2:874366/73‑1 (MQ=255)

                                                       gACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc              >  1:1276782/1‑76 (MQ=255)

                                                        gcGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa             <  2:1397619/75‑1 (MQ=255)

                                                        aCGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa             >  1:1091142/1‑76 (MQ=255)

                                                        aCGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa             <  1:3175104/76‑1 (MQ=255)

                                                         cGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg            <  2:430673/76‑1 (MQ=255)

                                                         cGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg            >  1:84551/1‑76 (MQ=255)

                                                         cGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg            >  2:4527911/1‑76 (MQ=255)

                                                         cGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg            <  2:613771/76‑1 (MQ=255)

                                                          ggCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg           <  2:2284607/76‑1 (MQ=255)

                                                          ggCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg           <  1:1693827/76‑1 (MQ=255)

                                                            cAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcg         >  2:4056288/1‑76 (MQ=255)

                                                            cAAGGCGTAACGCCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcg         >  1:2408352/1‑76 (MQ=255)

                                                             aaGGCGTAACGTCGCTGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc        >  1:3411882/1‑76 (MQ=255)

                                                             aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc        >  2:3130447/1‑76 (MQ=255)

                                                             aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc        <  2:3338865/76‑1 (MQ=255)

                                                              aGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCa       <  2:620210/76‑1 (MQ=255)

                                                                gCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa     <  2:4700223/76‑1 (MQ=255)

                                                                gCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa     <  1:1995979/76‑1 (MQ=255)

                                                                  gTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa   >  1:3171839/1‑76 (MQ=255)

                                                                   tAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt  <  2:1411991/76‑1 (MQ=255)

                                                                   tAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt  <  1:1516557/76‑1 (MQ=255)

CCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAT  >  NC_000913/3961517‑3961659

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