BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,115,116	+C	coding (2013/2028 nt)	b2974 ←	putative endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,115,115	1	.	C	99.4%	436.3 / ‑5.9	164	K672E (AAA→GAA)	b2974	putative endoglucanase
Reads supporting (aligned to +/- strand): ref base . (1/0); major base C (90/72); minor base A (1/0); total (92/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 179 total aligned reads displayed.

CGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCG  >  NC_000913/3115040‑3115185

cGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAttt                                                                         <  1:404208/76‑1 (MQ=255)

cGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAttt                                                                         >  1:4329434/1‑76 (MQ=255)

cGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAttt                                                                         >  2:1132836/1‑76 (MQ=255)

 gCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑c                                                                       <  1:1814330/76‑1 (MQ=255)

   aGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cct                                                                     <  1:2144448/76‑3 (MQ=255)

   aGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cct                                                                     >  1:3126119/1‑74 (MQ=255)

       ggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCgg                                                                  <  1:2410495/76‑1 (MQ=255)

        gATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGa                                                                 >  1:847598/1‑76 (MQ=255)

           tGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                              <  1:1309317/76‑1 (MQ=255)

            gCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCga                                                             <  2:1753864/76‑1 (MQ=255)

             cATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgag                                                            >  1:2797950/1‑76 (MQ=255)

             cATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgag                                                            >  2:1293645/1‑76 (MQ=255)

              aTAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgaga                                                           <  1:3913194/76‑1 (MQ=255)

              aTAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGAGCgaga                                                           >  2:2841696/1‑76 (MQ=255)

                aaGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATc                                                         >  1:3849347/1‑76 (MQ=255)

                aaGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATc                                                         <  1:2815176/76‑1 (MQ=255)

                aaGGATCGTTACCCACCTAAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATc                                                         <  1:2561770/76‑1 (MQ=255)

                 gggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCg                                                        <  2:2022867/75‑1 (MQ=255)

                  ggATCGTTCCCCCCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTCCCTGCGGATCGAGATCGc                                                       <  2:1329704/76‑1 (MQ=255)

                  ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGc                                                       <  1:1366089/76‑1 (MQ=255)

                  ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGc                                                       >  1:302585/1‑76 (MQ=255)

                  ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGAACGc                                                       >  2:262661/1‑76 (MQ=255)

                    aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCa                                                     >  1:3161377/1‑76 (MQ=255)

                    aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCa                                                     <  2:1433024/76‑1 (MQ=255)

                      cGTTGCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATa                                                   <  2:1379634/76‑1 (MQ=255)

                        ttACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc                                                 >  1:3872371/1‑76 (MQ=255)

                        ttACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc                                                 <  1:3988744/76‑1 (MQ=255)

                         tACCCACCTGAGACACATATTCAAAGCCGTGAAGGATCAGCAGCGGCATTTCCTGCGGATTGAGATCGTCATAgat                                                >  1:2104464/1‑74 (MQ=255)

                          aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa                                               >  1:3551945/1‑76 (MQ=255)

                           cccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                              <  1:4582499/76‑1 (MQ=255)

                            ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACt                                             >  1:1309209/1‑76 (MQ=255)

                            ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACt                                             >  2:4363293/1‑76 (MQ=255)

                             cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACtt                                            <  1:2718265/76‑1 (MQ=255)

                              cccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACtaa                                           <  2:1003164/75‑3 (MQ=255)

                              aCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTa                                           >  1:2911821/1‑76 (MQ=255)

                              aCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTa                                           >  2:1508732/1‑76 (MQ=255)

                               ccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaa                                          <  1:2016661/76‑1 (MQ=255)

                               ccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCCCCATCGCTACTGaa                                          >  2:108729/1‑76 (MQ=255)

                                cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                         <  2:4638361/76‑1 (MQ=255)

                                 tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAAt                                        >  1:1096450/1‑76 (MQ=255)

                                 tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAAt                                        <  1:2861618/76‑1 (MQ=255)

                                 tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAAt                                        >  1:4059224/1‑76 (MQ=255)

                                 tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGCATCGAGATCGCCATAGCTACTTAAAt                                        >  2:2734239/1‑76 (MQ=255)

                                  gcgTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATg                                       >  2:4711521/3‑76 (MQ=255)

                                   aGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                      <  1:4776376/76‑1 (MQ=255)

                                    gTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCt                                     >  1:1171127/1‑76 (MQ=255)

                                    gTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCt                                     >  2:3592814/1‑76 (MQ=255)

                                     tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGGGATCGCCATAGCTACTTAAATGCTc                                    <  1:1442613/76‑1 (MQ=255)

                                     tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                    >  2:1915430/1‑76 (MQ=255)

                                     tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                    >  2:1891177/1‑76 (MQ=255)

                                     tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                    <  1:3738731/76‑1 (MQ=255)

                                         atatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACa                                >  1:2777516/1‑76 (MQ=255)

                                         atatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCACAGCTACTTAAATGCTCCAcc                                >  2:531958/1‑75 (MQ=255)

                                           atTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                              >  1:4587482/1‑76 (MQ=255)

                                             tCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGACAGACATAGCTATCGAAATGCGCAACAaaga                            >  1:1579743/1‑76 (MQ=255)

                                              cAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                           <  1:1034615/76‑1 (MQ=255)

                                              cAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                           <  1:1906517/76‑1 (MQ=255)

                                              cAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                           <  2:626713/76‑1 (MQ=255)

                                               aaaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATg                          >  1:3416244/1‑76 (MQ=255)

                                                aaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGc                         >  1:4005207/1‑76 (MQ=255)

                                                 aGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                        <  2:4315960/76‑1 (MQ=255)

                                                 aGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                        >  1:452415/1‑76 (MQ=255)

                                                  gCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGccc                       >  1:4242008/1‑76 (MQ=255)

                                                  gCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGccc                       >  2:2864069/1‑76 (MQ=255)

                                                    cGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGc                     <  2:2098827/76‑1 (MQ=255)

                                                    cGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGc                     >  1:2027538/1‑76 (MQ=255)

                                                       tAAGGATCAGCATCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGATATGCCCGCaaa                  >  2:2070578/1‑76 (MQ=255)

                                                       tAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaca                  >  1:661586/1‑74 (MQ=255)

                                                       tAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaa                  >  1:554180/1‑76 (MQ=255)

                                                       tAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaa                  <  2:4439782/76‑1 (MQ=255)

                                                        aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCTTAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                 <  1:2938590/76‑1 (MQ=255)

                                                        aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaca                 >  2:2238054/1‑74 (MQ=255)

                                                        aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAATTAAACGCTCACCAGAGATGCCCGCaaaa                 >  2:4235214/1‑76 (MQ=255)

                                                        aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCAGAGCTACTGAAATGCGCAACAGAGATGCCCGCaaaa                 >  2:4106290/1‑76 (MQ=255)

                                                        aaGGATCAGCAGCGGCATTTCCGGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                 <  2:1120179/76‑1 (MQ=255)

                                                        aaGGATCAGAAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                 >  2:3331257/1‑76 (MQ=255)

                                                        aaGGACCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                 <  2:676728/76‑1 (MQ=255)

                                                        aaGGACCAGCAGCGGCATCTCCAGCGGACCGAGAGCGCCAGAGCCACTTAAATGCTCAAAAGACATGCCACCaaaa                 >  2:2350412/1‑76 (MQ=255)

                                                         aGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                <  2:889225/76‑1 (MQ=255)

                                                          ggATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATc               >  1:1937414/1‑76 (MQ=255)

                                                          ggATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATc               <  1:831870/76‑1 (MQ=255)

                                                           gATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc              <  2:3594920/76‑1 (MQ=255)

                                                           gATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc              >  1:1996674/1‑76 (MQ=255)

                                                            aTCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCg             >  2:1623633/1‑76 (MQ=255)

                                                            aTCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCg             <  1:380379/76‑1 (MQ=255)

                                                              cagcagCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGa           >  2:1832708/1‑76 (MQ=255)

                                                               agcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGAt          >  2:3244069/1‑76 (MQ=255)

                                                                gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATg         >  1:1848792/1‑76 (MQ=255)

                                                                gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATg         >  2:3741974/1‑76 (MQ=255)

                                                                 cagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg        >  2:1671592/1‑76 (MQ=255)

                                                                 cagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg        >  1:2751719/1‑76 (MQ=255)

                                                                  agcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa       >  2:3373634/1‑76 (MQ=255)

                                                                    cGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAg     <  2:252097/76‑1 (MQ=255)

                                                                     ggCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGCAgc    >  1:1079776/1‑75 (MQ=255)

                                                                     ggCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAgc    >  2:2331528/1‑75 (MQ=255)

                                                                     ggCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt    <  2:2518921/76‑1 (MQ=255)

                                                                      gCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc   <  1:3592814/76‑1 (MQ=255)

                                                                      gCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc   <  1:3606034/76‑1 (MQ=255)

                                                                      gCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc   <  2:1829152/76‑1 (MQ=255)

                                                                       cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGACGCCCGCAAAATCCGGATGGCAGCCg  >  2:579969/1‑76 (MQ=255)

CGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCG  >  NC_000913/3115040‑3115185

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