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%	376.2 / NA	98	L54I (CTT→ATT)	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (31/67); total (31/67)

Only 100 of 104 total aligned reads displayed.

CCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGC  >  NC_000913/803597‑803726

     aaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGAtaaaattaaaat                                                           >  1:1442594/1‑68 (MQ=255)

      aaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATaaaattaaaatt                                                          >  1:3720885/1‑68 (MQ=255)

         gCACTATTGCCTGGTCTATCCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAg                                                       <  1:240871/68‑1 (MQ=255)

          cACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGt                                                      <  1:2289995/68‑1 (MQ=255)

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

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

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

                 gCCTGGTCTTTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                               <  1:2666672/68‑1 (MQ=255)

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

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

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

                   cTGGTCTATTCTTTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                             >  1:2539316/1‑68 (MQ=255)

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

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

                     ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTa                                           >  1:3681667/1‑68 (MQ=255)

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

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

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

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

                       tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTCTGATTCATTAGc                                         >  1:684120/1‑68 (MQ=255)

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

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

                        cTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTCTGATTCATTAGcc                                        <  1:1410182/68‑1 (MQ=255)

                         tATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAAGTAAAATTAAGTTTGATTCATTAGCct                                       <  1:3000364/68‑1 (MQ=255)

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

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

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

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

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

                           ttCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctct                                     <  1:4184054/68‑1 (MQ=255)

                           ttCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctct                                     <  1:776251/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                            tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCACTAGCCtctc                                    >  1:1284529/1‑68 (MQ=255)

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

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

                                tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                                >  1:1718132/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

                                 ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCACTAGCCTCTCACGat                               <  1:3700804/68‑1 (MQ=255)

                                 ctTCGCATAATACGCCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:3805711/68‑1 (MQ=255)

                                    cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATAtt                            <  1:3929526/68‑1 (MQ=255)

                                    cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATAtt                            <  1:3493892/68‑1 (MQ=255)

                                      cATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGCTATTAc                          >  1:2808488/1‑68 (MQ=255)

                                       ataataCGTCCGGAAATATGGATAAAATTAAAAGTAAGTTTGATTCATTAGCCTCTCACGATATTGcc                         <  1:746332/68‑1 (MQ=255)

                                         aataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt                       <  1:1031286/68‑1 (MQ=255)

                                          ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt                      >  1:23236/1‑68 (MQ=255)

                                            aCGTCCGGAAATATGGATAAAATTAAATTTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt                    >  1:3808453/1‑68 (MQ=255)

                                             cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt                   >  1:1916098/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTCTTGGt                <  1:2420050/68‑1 (MQ=255)

                                                ccGGAAATATGGATAAAATTAAAATTAAGTTTGAGTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:7802/68‑1 (MQ=255)

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

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

                                                    aaaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTg            >  1:3655214/1‑68 (MQ=255)

                                                    aaaTATGGATAAAATTAAAATTAAGTTTGATTAATTAGCCTNTCACGATATTACCTTGGTGGTTAtgg            >  1:3256281/1‑66 (MQ=255)

                                                     anTATGGATAAANNTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:3192273/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                              aaaattaaaattAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGc  <  1:572912/68‑1 (MQ=255)

CCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGC  >  NC_000913/803597‑803726

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