BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,942,094	C→T	H35Y (CAT→TAT)	ygfA →	predicted ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,942,094	0	C	T	98.8%	659.3 / ‑1.5	162	H35Y (CAT→TAT)	ygfA	predicted ligase
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (101/59); total (103/59)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 178 total aligned reads displayed.

GCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCG  >  REL606/2942021‑2942166

gCCGAAAAAACAGCAGGAAAGGGATAAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCatat                                                                        <  2:15331/76‑1 (MQ=255)

 ccGGAACAACAGCAGGAAAGGGGTCAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCAtata                                                                       <  1:7040396/76‑1 (MQ=255)

  cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATAc                                                                      >  1:455591/1‑76 (MQ=255)

  cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATAc                                                                      >  2:1485384/1‑76 (MQ=255)

  cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATAc                                                                      >  2:5003581/1‑76 (MQ=255)

   ggAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACg                                                                     >  2:364584/1‑76 (MQ=255)

   ggAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACg                                                                     >  2:5782507/1‑76 (MQ=255)

     aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGt                                                                   >  2:6614150/1‑76 (MQ=255)

     aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGt                                                                   >  2:3197458/1‑76 (MQ=255)

      acaacaGAAGGAAATGGTAAAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTc                                                                  <  1:4766366/76‑1 (MQ=255)

      acaaaaGAAGGAAAGGGGTAAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTc                                                                  <  1:4902163/76‑1 (MQ=255)

      aaaacaGCAGAAAGTGGGAAAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTc                                                                  <  1:1571549/74‑1 (MQ=255)

       caacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCg                                                                 >  1:6680348/1‑76 (MQ=255)

         acaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCt                                                               >  2:1616295/1‑76 (MQ=255)

         acaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCt                                                               >  2:1988780/1‑76 (MQ=255)

            gcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTa                                                            >  2:3002122/1‑76 (MQ=255)

             cagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTAt                                                           >  1:6443604/1‑76 (MQ=255)

             cagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTAt                                                           >  2:140223/1‑76 (MQ=255)

             aagGAACTGCGTAAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTAt                                                           <  1:804060/75‑1 (MQ=255)

              gggAAAGGGGAAAAAAAGCCGCTACCCGGATGATGACTTACCCCCCGGTGGTGATGGCATATACGGTCGCTGTAtt                                                          <  1:3810100/75‑1 (MQ=255)

              agGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTAtt                                                          >  2:6535906/1‑76 (MQ=255)

              agGAAATGGAAAAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTAtt                                                          <  1:5048165/76‑1 (MQ=255)

               agAAATGGGGCAAAAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTc                                                         <  2:7015063/75‑1 (MQ=255)

                gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTcc                                                        >  2:2977375/1‑76 (MQ=255)

                gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTcc                                                        >  1:6117810/1‑76 (MQ=255)

                 aaaTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCct                                                       >  1:832671/1‑76 (MQ=255)

                 aaaTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCct                                                       >  2:2758232/1‑76 (MQ=255)

                    tGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCctctc                                                    >  1:2503382/1‑76 (MQ=255)

                       gTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCttt                                                 >  1:5459854/1‑76 (MQ=255)

                          aacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGa                                              >  1:2963839/1‑76 (MQ=255)

                          aacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGa                                              >  2:5925861/1‑76 (MQ=255)

                          aacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGa                                              >  2:939056/1‑76 (MQ=255)

                           aaaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGAt                                             <  2:3019445/74‑1 (MQ=255)

                           aaaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGAt                                             <  2:4705552/74‑1 (MQ=255)

                            aaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATg                                            <  1:3349827/75‑1 (MQ=255)

                             aaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATgg                                           >  1:1012273/1‑76 (MQ=255)

                             aaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATgg                                           >  1:885915/1‑76 (MQ=255)

                             aaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATgg                                           >  2:5538328/1‑76 (MQ=255)

                              aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGc                                          <  1:3829360/76‑1 (MQ=255)

                              aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGc                                          >  1:6998623/1‑76 (MQ=255)

                                 cGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGaa                                       <  1:1101303/76‑1 (MQ=255)

                                  gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAAc                                      >  1:3626360/1‑76 (MQ=255)

                                  gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAAc                                      <  1:6535906/76‑1 (MQ=255)

                                   cTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACt                                     <  2:4560391/76‑1 (MQ=255)

                                   cTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACt                                     >  1:838015/1‑76 (MQ=255)

                                    tACCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTc                                    <  1:6753313/76‑1 (MQ=255)

                                     aCCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCg                                   >  1:2319213/1‑76 (MQ=255)

                                     aCCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCg                                   <  1:2670012/76‑1 (MQ=255)

                                     aCCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCg                                   <  2:4406981/76‑1 (MQ=255)

                                     aCCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCg                                   >  1:3341785/1‑76 (MQ=255)

                                     aCCCGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCg                                   >  1:5723501/1‑76 (MQ=255)

                                        cGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGaca                                >  2:5656930/1‑76 (MQ=255)

                                        cGGATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGaca                                >  1:2066689/1‑76 (MQ=255)

                                        cGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGaca                                >  1:4641101/1‑76 (MQ=255)

                                         ggATGATGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGacac                               >  1:2227383/1‑76 (MQ=255)

                                          gatgatGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACAcc                              <  1:4788758/76‑1 (MQ=255)

                                          gatgatGACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACAcc                              >  2:2261334/1‑76 (MQ=255)

                                           atgatgACTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACAccc                             <  2:6694686/76‑1 (MQ=255)

                                            tgatgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                            >  1:4013252/1‑76 (MQ=255)

                                            tgatgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                            <  2:565954/76‑1 (MQ=255)

                                            tgatgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                            <  2:6307238/76‑1 (MQ=255)

                                            tgatgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                            >  1:3739838/1‑76 (MQ=255)

                                            tgatgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                            >  1:628162/1‑76 (MQ=255)

                                            tgatgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                            >  1:3536991/1‑76 (MQ=255)

                                            tgatgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                            <  2:6874686/76‑1 (MQ=255)

                                             gatgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTTTTCCTCTCTTTTGATGGCGAACTCGACACCCAg                           >  1:1389749/1‑76 (MQ=255)

                                              atgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc                          >  1:12796/1‑76 (MQ=255)

                                              atgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc                          >  1:3182646/1‑76 (MQ=255)

                                              atgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc                          >  2:6689988/1‑76 (MQ=255)

                                               tgaCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGcc                         <  1:6362035/76‑1 (MQ=255)

                                                 aCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCAc                       <  1:1451952/76‑1 (MQ=255)

                                                 aCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCAc                       >  2:4902339/1‑76 (MQ=255)

                                                 aCTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCAc                       >  1:1381527/1‑76 (MQ=255)

                                                  cTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACt                      >  1:5921317/1‑76 (MQ=255)

                                                  cTTATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACt                      >  1:3612069/1‑76 (MQ=255)

                                                   ttATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTc                     >  1:1025171/1‑76 (MQ=255)

                                                    tATCCCCCGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCa                    >  1:3373388/1‑76 (MQ=255)

                                                       cccccGGTGGTGATGGCATATACGGTCGCTTTATTCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCACTCATAg                 <  1:5891773/76‑1 (MQ=255)

                                                       cccccGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAg                 >  1:6886331/1‑76 (MQ=255)

                                                       cccccGGGGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGAGGGCGAACTCGACCCCCAGCCACTCATAg                 <  1:4027824/76‑1 (MQ=255)

                                                        ccccGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCACTCATAGa                <  1:396387/76‑1 (MQ=255)

                                                        ccccGGTGGTGATGCCATAGCCGGTCGCTGTTTTCCCCTCTTTGGAGGGCGAACTCGACCCCCAGCCACTCATAGa                <  2:2537721/76‑1 (MQ=255)

                                                         ccgggtggtGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCACTCATAGaa               <  2:5732707/73‑1 (MQ=255)

                                                         cccGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCACTCATAGaa               <  2:5263581/76‑1 (MQ=255)

                                                         cccGGTGGTGATGGCATATACGGTCGCTGGATTCCTCTCTTTTGATGGCGAACTCGACACCCAAGCACTCATAGaa               >  2:216329/1‑76 (MQ=255)

                                                          ccGGTGGTGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGAGGGCGAACTCGACCCCCAGCCACTCATAGaac              <  2:5433071/76‑1 (MQ=255)

                                                            ggtggtGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCACTCATAGaacaa            <  1:2194331/76‑1 (MQ=255)

                                                            ggtggtGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaacaa            >  2:3519216/1‑76 (MQ=255)

                                                            ggtggtGATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGAACACGAGGCACTCATAGaacac            >  2:1217143/1‑75 (MQ=255)

                                                                gtgATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCACTCATAGAACAActct        <  1:4271813/76‑1 (MQ=255)

                                                                gtgATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCACTCATAGAACAActct        <  1:3076465/76‑1 (MQ=255)

                                                                gtgATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAActct        >  2:1795805/1‑76 (MQ=255)

                                                                gtgATGCCCTATCCGGTCGTTGTTTTCCTCTCTTTTGACGGCGACCTCGACCCCCAGCCACTCATAGAACAActct        <  2:1013166/76‑1 (MQ=255)

                                                                 tgATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTg       >  2:5346176/1‑76 (MQ=255)

                                                                 tgATGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTg       >  2:2420704/1‑76 (MQ=255)

                                                                 tgATGGCATATACGGTCGCTGTATCCCTCTCTTTTGATGGCGAACTCGACCCCCAGCCACTCATAGAACAACTCTg       <  2:5872480/76‑1 (MQ=255)

                                                                   aTGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgc     >  1:6952994/1‑76 (MQ=255)

                                                                   aTGGCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgc     >  1:4519242/1‑76 (MQ=255)

                                                                    tGGCATATCCGGTCCTTGTTTCCCTCTCTTTTGATGGCCAACCCGACACCCAGCCACTCATAGAACAACTCTGgcg    <  1:4182311/76‑1 (MQ=255)

                                                                     ggCATATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgcgc   >  2:2583314/1‑76 (MQ=255)

GCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCG  >  REL606/2942021‑2942166

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