BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	803,662	C→A	L54I (CTT→ATT)	ybhJ →	putative hydratase YbhJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,662	0	C	A	100.0%	401.7 / NA	104	L54I (CTT→ATT)	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (34/70); total (34/70)

Only 100 of 105 total aligned reads displayed.

CCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGC  >  NC_000913/803597‑803726

ccAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAAtt                                                                >  1:3026031/1‑68 (MQ=255)

   aaaaaGGCACTATTACCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaaa                                                             <  1:2407937/68‑1 (MQ=255)

         gCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAg                                                       <  1:2929111/68‑1 (MQ=255)

              aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGa                                                  >  1:1747055/1‑68 (MQ=255)

                tGCCTGGTCTATTCTCTCTTCGCATANTACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                                <  1:2718502/68‑1 (MQ=255)

                tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                                >  1:532640/1‑68 (MQ=255)

                tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                                <  1:398968/68‑1 (MQ=255)

                 gCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                               <  1:1353861/68‑1 (MQ=255)

                 gCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                               <  1:827171/68‑1 (MQ=255)

                  ccTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                              >  1:3277327/1‑68 (MQ=255)

                  ccTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                              <  1:547740/68‑1 (MQ=255)

                   cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                             >  1:468041/1‑68 (MQ=255)

                   cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                             >  1:1852737/1‑68 (MQ=255)

                   cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                             >  1:120382/1‑68 (MQ=255)

                   cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                             >  1:571425/1‑68 (MQ=255)

                   cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATGCAt                                             >  1:2463780/1‑68 (MQ=255)

                      gTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAg                                          <  1:861035/68‑1 (MQ=255)

                       tCTATTCTCTCTTCGCATAATACGTCGGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                         <  1:2727158/68‑1 (MQ=255)

                       tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                         <  1:1276831/68‑1 (MQ=255)

                       tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                         >  1:2787032/1‑68 (MQ=255)

                       tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                         >  1:1966833/1‑68 (MQ=255)

                       tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                         >  1:3689/1‑68 (MQ=255)

                       tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                         >  1:1296195/1‑68 (MQ=255)

                        cTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                                        <  1:2357359/68‑1 (MQ=255)

                        cTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                                        >  1:2269695/1‑68 (MQ=255)

                          aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                      <  1:1046885/68‑1 (MQ=255)

                          aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                      <  1:1734800/68‑1 (MQ=255)

                           ttCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAACTAAGTTTGATTCATTAGCctct                                     <  1:2259588/68‑1 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAACTTAAGTTTGATTCATTAGCCtctc                                    >  1:1995384/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:3144910/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:160408/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1204244/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:3150305/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1341128/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:3286065/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1345260/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1667531/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1476614/1‑68 (MQ=255)

                            tctctcTTCGCATAATACGTCCGGAAATATGCATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:3263087/1‑68 (MQ=255)

                               ctgttCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                                 <  1:3001833/65‑1 (MQ=255)

                               ctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                                 <  1:1038620/68‑1 (MQ=255)

                               ctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                                 <  1:1927362/68‑1 (MQ=255)

                                tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                                <  1:2895521/68‑1 (MQ=255)

                                tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                                <  1:3257967/68‑1 (MQ=255)

                                tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                                <  1:3200491/68‑1 (MQ=255)

                                tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                                <  1:2510454/68‑1 (MQ=255)

                                 ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:277867/68‑1 (MQ=255)

                                 ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:27346/68‑1 (MQ=255)

                                 ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:3257744/68‑1 (MQ=255)

                                 ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               >  1:1222540/1‑68 (MQ=255)

                                 ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGGTTGATTCATTAGCCTCTCACGat                               <  1:2916582/68‑1 (MQ=255)

                                 ctTCGCACAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:1673526/68‑1 (MQ=255)

                                       ataataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                         <  1:953442/68‑1 (MQ=255)

                                       ataataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                         <  1:2683120/68‑1 (MQ=255)

                                           taCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg                     <  1:1221105/68‑1 (MQ=255)

                                           taCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg                     <  1:2163037/68‑1 (MQ=255)

                                            aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACATttgt                    >  1:908282/1‑68 (MQ=255)

                                             cGTCCGGAAATATGGATANAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt                   >  1:2580062/1‑68 (MQ=255)

                                               tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTgg                 >  1:500811/1‑68 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:638620/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:1165507/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:664259/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:555934/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:2080452/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:3403603/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:820994/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:1128102/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:2147445/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:2443110/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:2708341/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:2994607/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:1399136/68‑1 (MQ=255)

                                                  ggAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt              <  1:1544655/68‑1 (MQ=255)

                                                  ggAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt              <  1:440691/68‑1 (MQ=255)

                                                  ggAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt              <  1:2086176/68‑1 (MQ=255)

                                                  ggAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt              <  1:1311055/68‑1 (MQ=255)

                                                    aaaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTg            <  1:1673065/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:3228638/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:2801526/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           >  1:1575872/1‑68 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:2054627/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:2999956/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:3377138/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:3375274/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:1563413/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:90696/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:1727688/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:2192383/68‑1 (MQ=255)

                                                     aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:2437696/68‑1 (MQ=255)

                                                      atatGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTa          >  1:2551455/1‑68 (MQ=255)

                                                      atatGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTa          >  1:622084/1‑68 (MQ=255)

                                                       tatGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTAc         <  1:3334364/68‑1 (MQ=255)

                                                       tatGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTAc         <  1:2826524/68‑1 (MQ=255)

                                                          ggATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGa      >  1:882723/1‑68 (MQ=255)

                                                            aTAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:558919/68‑1 (MQ=255)

                                                            aTAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:3300684/68‑1 (MQ=255)

                                                            aTAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:2119450/68‑1 (MQ=255)

                                                            aTAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:3149225/68‑1 (MQ=255)

                                                            aTAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:312224/68‑1 (MQ=255)

                                                            aTAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:1807940/68‑1 (MQ=255)

CCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGC  >  NC_000913/803597‑803726

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