BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	795,728	(A)_7→6	coding (602/993 nt)	ECB_00737 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	795,722	0	A	.	98.0%	746.5 / 7.5	151	coding (608/993 nt)	ECB_00737	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/3); new base . (74/74); total (74/77)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 177 total aligned reads displayed.

TTTTGCTGGGTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTAAAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAATATACGG  >  REL606/795647‑795794

 tttGCTGGGTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaa                                                                         <  1:3959782/76‑1 (MQ=255)

  ttGCTGGGTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaaa                                                                        >  1:5342165/1‑76 (MQ=255)

    gCTGGGTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaaaaa                                                                      >  2:1514546/1‑76 (MQ=255)

     cTGGGTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaaaaaa                                                                     >  1:5524860/1‑76 (MQ=255)

     cTGGGTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaaaaaa                                                                     >  1:2245040/1‑76 (MQ=255)

      tGGGTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaaaaaag                                                                    >  1:2875556/1‑75 (MQ=255)

       gggTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaaaaaagt                                                                   >  2:1207042/1‑74 (MQ=255)

       gggTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaaaaaagt                                                                   >  1:6884759/1‑74 (MQ=255)

       gggTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTaaaaaaga                                                                   >  2:1845553/1‑74 (MQ=255)

          tggtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑ACAAAAGTAtt                                                               <  2:5680684/76‑1 (MQ=255)

          tggtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTAtt                                                               <  2:6519918/76‑1 (MQ=255)

          tggtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTAtt                                                               <  1:4750540/76‑1 (MQ=255)

          tggtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTAtt                                                               >  2:1093394/1‑76 (MQ=255)

           ggtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTg                                                              >  1:6945574/1‑76 (MQ=255)

           ggtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTg                                                              <  2:7855779/76‑1 (MQ=255)

           ggtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTg                                                              <  1:6406134/76‑1 (MQ=255)

            gtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGc                                                             <  2:1700182/76‑1 (MQ=255)

            gtAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGc                                                             >  2:3666284/1‑76 (MQ=255)

             tAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTAAAAAAAGTATTGc                                                             <  1:328431/76‑1 (MQ=255)

             tAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCg                                                            <  2:196942/76‑1 (MQ=255)

              aGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGa                                                           <  2:4845960/76‑1 (MQ=255)

              aGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGa                                                           >  1:7891077/1‑76 (MQ=255)

                 aTACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATAc                                                        <  1:4159377/76‑1 (MQ=255)

                 aTACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATAc                                                        <  1:5814441/76‑1 (MQ=255)

                 aTACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATAc                                                        >  2:6094653/1‑76 (MQ=255)

                 aTACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATAc                                                        >  2:7663813/1‑76 (MQ=255)

                 aTACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATAc                                                        >  1:1364316/1‑76 (MQ=255)

                  tACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACa                                                       <  1:5943697/76‑1 (MQ=255)

                  tACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACa                                                       <  1:4832975/76‑1 (MQ=255)

                   aCAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAg                                                      <  1:4407765/76‑1 (MQ=255)

                   aCAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAg                                                      >  1:867125/1‑76 (MQ=255)

                   aCAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAg                                                      >  2:6867796/1‑76 (MQ=255)

                    cAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGa                                                     <  2:2364923/76‑1 (MQ=255)

                     aGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGAc                                                    <  1:1101173/76‑1 (MQ=255)

                     aGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGAc                                                    <  1:7088107/76‑1 (MQ=255)

                      gAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACg                                                   >  1:6018727/1‑76 (MQ=255)

                      gAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACg                                                   <  1:4818408/76‑1 (MQ=255)

                      aaaaTAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACg                                                   <  2:1355672/75‑1 (MQ=255)

                       aaaTAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGc                                                  >  1:1685354/1‑76 (MQ=255)

                       aaaTAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGc                                                  >  1:8001077/1‑76 (MQ=255)

                         ataataAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGccc                                                >  1:424119/1‑76 (MQ=255)

                          taataaAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCcca                                               <  1:4116969/76‑1 (MQ=255)

                          taataaAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCcca                                               <  1:1465600/76‑1 (MQ=255)

                           aataaAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCccac                                              <  1:626872/76‑1 (MQ=255)

                           aataaAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCccac                                              >  1:7679996/1‑76 (MQ=255)

                            ataaAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCccacc                                             >  2:3294524/1‑76 (MQ=255)

                             taaAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCccacca                                            >  2:675578/1‑76 (MQ=255)

                              aaaaTGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCcaccac                                           <  1:80716/76‑1 (MQ=255)

                              aaaaTGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCAcaccac                                           >  2:7310351/1‑76 (MQ=255)

                               aaaTGTTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACcaca                                          <  2:5343485/76‑1 (MQ=255)

                                aaTGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACaa                                         >  2:1268070/1‑76 (MQ=255)

                                  tGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAAt                                       >  1:5942746/1‑76 (MQ=255)

                                    cTTAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACACATTc                                     >  2:4158903/1‑76 (MQ=255)

                                     ttAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCt                                    >  2:796657/1‑76 (MQ=255)

                                     ttAATAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCt                                    >  2:1575971/1‑76 (MQ=255)

                                       aaTAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTgg                                  >  2:6492950/1‑76 (MQ=255)

                                        aTAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGt                                 >  2:1281416/1‑76 (MQ=255)

                                        aTAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGt                                 <  2:5845158/76‑1 (MQ=255)

                                         tAGCCATTTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTa                                <  2:1765173/76‑1 (MQ=255)

                                         tAGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTa                                <  2:4215637/76‑1 (MQ=255)

                                          aGCCATGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTaa                               >  1:2467728/1‑76 (MQ=255)

                                              aTGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTg                           >  1:2298574/1‑76 (MQ=255)

                                              aTGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTg                           >  2:2992303/1‑76 (MQ=255)

                                               tGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGt                          >  1:3046466/1‑76 (MQ=255)

                                               tGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGt                          >  1:7270197/1‑76 (MQ=255)

                                               tGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGt                          >  1:8202367/1‑76 (MQ=255)

                                               tGTATAATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGt                          >  2:5389447/1‑76 (MQ=255)

                                                 tataATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGtat                        >  2:4040686/1‑76 (MQ=255)

                                                 tataATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGtat                        >  1:3574643/1‑76 (MQ=255)

                                                  ataATGTTTACCACGATAGCCATTT‑AAAAAAGTATTGCAAAACAGACGCCCACCACAAATTCTGGTAACCTGtata                       <  2:564878/76‑1 (MQ=255)

                                                  ataATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGtata                       >  1:1057502/1‑76 (MQ=255)

                                                  ataATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGtata                       >  1:2886832/1‑76 (MQ=255)

                                                  ataATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGtata                       >  2:1021460/1‑76 (MQ=255)

                                                  ataATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGtata                       <  1:6699181/76‑1 (MQ=255)

                                                   taATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTatat                      >  1:7462458/1‑76 (MQ=255)

                                                   taATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTatat                      >  2:7515900/1‑76 (MQ=255)

                                                   taATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTatat                      >  2:2624240/1‑76 (MQ=255)

                                                   taATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTatat                      <  1:5634036/76‑1 (MQ=255)

                                                   taATGGTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACCCCCACCACAAATTCTGGTAACCTGTatat                      <  1:181084/76‑1 (MQ=255)

                                                       ggTTAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGa                  >  1:1232578/1‑76 (MQ=255)

                                                          tAACACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGtt               <  2:2253426/76‑1 (MQ=255)

                                                           aaCACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTg              <  2:5158561/76‑1 (MQ=255)

                                                           aaCACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTg              >  1:4394324/1‑76 (MQ=255)

                                                           aaCACGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTg              >  2:7500236/1‑76 (MQ=255)

                                                              acGATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGt           <  2:5276534/76‑1 (MQ=255)

                                                              acGAAAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGt           <  1:2499179/76‑1 (MQ=255)

                                                               caaTAGCTATTT‑AAAAAATTTTGGCGATGCAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTa          <  2:6365998/74‑1 (MQ=255)

                                                                gATAGCCATTT‑AAAAAAGTATTGCGATACAGCCGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTaa         <  1:257809/76‑1 (MQ=255)

                                                                gATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTaa         <  1:7819608/76‑1 (MQ=255)

                                                                gATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTaa         <  1:7188983/76‑1 (MQ=255)

                                                                gATAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTaa         <  1:5837292/76‑1 (MQ=255)

                                                                 aTAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAat        >  2:4586169/1‑76 (MQ=255)

                                                                 aTAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAat        >  2:2163890/1‑76 (MQ=255)

                                                                 aTAGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAat        >  2:3422778/1‑76 (MQ=255)

                                                                   aGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAatat      <  1:1948102/76‑1 (MQ=255)

                                                                   aGCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAatat      <  2:6686562/76‑1 (MQ=255)

                                                                    gCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAAtata     <  1:5496975/76‑1 (MQ=255)

                                                                    gCCATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAAtata     >  2:8149637/1‑76 (MQ=255)

                                                                      cATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAATATACg   <  1:3793660/76‑1 (MQ=255)

                                                                      cATTT‑AAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAATATACg   >  2:3136389/1‑76 (MQ=255)

TTTTGCTGGGTGGTAGCATACAGAAATAATAAAATGCTTAATAGCCATGTATAATGGTTAACACGATAGCCATTTAAAAAAAGTATTGCGATACAGACGCCCACCACAAATTCTGGTAACCTGTATATTGGAGTTGCGTAATATACGG  >  REL606/795647‑795794

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