BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	870,805	C→T	F2F (TTC→TTT)	yliH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	870,805	0	C	T	99.5%	800.2 / ‑5.9	194	F2F (TTC→TTT)	yliH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (124/69); total (124/70)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

CAGCATTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTCGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACGTGCAA  >  REL606/870732‑870877

cAGCATTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGt                                                                        <  1:4467314/76‑1 (MQ=255)

cAGCATTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGt                                                                        <  2:2539811/76‑1 (MQ=255)

  gCATTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTttgttg                                                                      <  2:264058/76‑1 (MQ=255)

   cATTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTga                                                                     <  1:4152712/76‑1 (MQ=255)

   cATTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTga                                                                     <  1:5967163/76‑1 (MQ=255)

   cATTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTga                                                                     >  1:4533518/1‑76 (MQ=255)

    aTTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgac                                                                    >  2:246045/1‑76 (MQ=255)

     ttttGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgaca                                                                   >  1:4655116/1‑76 (MQ=255)

     ttttGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgaca                                                                   >  2:5340406/1‑76 (MQ=255)

      tttGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgacag                                                                  >  2:1260558/1‑76 (MQ=255)

      tttGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgacag                                                                  <  1:291448/76‑1 (MQ=255)

       ttGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgacaga                                                                 <  2:726413/76‑1 (MQ=255)

       ttGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgacaga                                                                 <  2:4343523/76‑1 (MQ=255)

       ttGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgacaga                                                                 >  1:5426063/1‑76 (MQ=255)

       ttGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgacaga                                                                 >  1:6053584/1‑76 (MQ=255)

        tGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTgacagac                                                                >  1:1621873/1‑76 (MQ=255)

        tGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTCGTTgacagac                                                                <  2:269019/76‑1 (MQ=255)

          cTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGacagacag                                                              >  1:644820/1‑76 (MQ=255)

           tCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGc                                                             <  2:2026475/76‑1 (MQ=255)

           tCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGc                                                             >  1:3823784/1‑76 (MQ=255)

           tCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGc                                                             <  1:3991544/76‑1 (MQ=255)

              aaaaaTAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGaa                                                          >  2:1403142/1‑76 (MQ=255)

              aaaaaTAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGaa                                                          >  1:1340891/1‑76 (MQ=255)

               aaaaTAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAAt                                                         <  2:5134943/76‑1 (MQ=255)

               aaaaTAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAAt                                                         <  1:232907/76‑1 (MQ=255)

               aaaaTAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAAt                                                         <  1:3686599/76‑1 (MQ=255)

                aaaTATCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATc                                                        >  1:4058265/1‑76 (MQ=255)

                aaaTAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATc                                                        <  1:2125128/76‑1 (MQ=255)

                 aaTAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCg                                                       <  2:925277/76‑1 (MQ=255)

                   tAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGAt                                                     <  1:2772153/76‑1 (MQ=255)

                    aGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATc                                                    <  1:2825893/76‑1 (MQ=255)

                      ccATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctg                                                  >  1:954763/1‑76 (MQ=255)

                       cATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctgc                                                 >  2:4035342/1‑76 (MQ=255)

                       cATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctgc                                                 >  2:3063338/1‑76 (MQ=255)

                       cATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctgc                                                 >  1:6118672/1‑76 (MQ=255)

                        aTACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctgct                                                >  2:2765965/1‑76 (MQ=255)

                        aTACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctgct                                                >  1:5473682/1‑76 (MQ=255)

                        aTACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctgct                                                >  2:2104325/1‑76 (MQ=255)

                        aTACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctgct                                                >  2:3715413/1‑76 (MQ=255)

                        aTACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATctgct                                                >  2:5818213/1‑76 (MQ=255)

                           cTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGaa                                             >  2:5805057/1‑76 (MQ=255)

                           cTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGaa                                             >  1:2625780/1‑76 (MQ=255)

                           cTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGaa                                             >  2:3808822/1‑76 (MQ=255)

                            tATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAAc                                            <  2:4267310/76‑1 (MQ=255)

                             aTTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAAcc                                           >  1:3933667/1‑76 (MQ=255)

                               ttAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCgg                                         >  1:13357/1‑76 (MQ=255)

                                tAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGt                                        <  1:3317246/76‑1 (MQ=255)

                                tAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGt                                        <  2:2272127/76‑1 (MQ=255)

                                tAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGt                                        >  2:311698/1‑76 (MQ=255)

                                tAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGt                                        >  1:1005920/1‑76 (MQ=255)

                                tAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGt                                        <  1:3420475/76‑1 (MQ=255)

                                 aaTTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGtt                                       >  2:2642632/1‑76 (MQ=255)

                                 aaTTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGtt                                       <  1:175407/76‑1 (MQ=255)

                                 aaTTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGtt                                       >  2:3310321/1‑76 (MQ=255)

                                   ttGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGa                                     >  2:493128/1‑76 (MQ=255)

                                   ttGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGa                                     <  1:5835181/76‑1 (MQ=255)

                                   ttGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGa                                     >  2:335196/1‑76 (MQ=255)

                                    tGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGAt                                    >  1:4567479/1‑76 (MQ=255)

                                     gCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATc                                   <  1:1809163/76‑1 (MQ=255)

                                     gCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATc                                   >  2:2234580/1‑76 (MQ=255)

                                      caacaaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCg                                  <  1:5799664/76‑1 (MQ=255)

                                      caacaaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCg                                  <  2:4768706/76‑1 (MQ=255)

                                      caacaaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCg                                  >  2:531501/1‑76 (MQ=255)

                                      caacaaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCg                                  <  2:2871980/76‑1 (MQ=255)

                                        acaaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGac                                <  1:4707020/76‑1 (MQ=255)

                                         caaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGacg                               >  2:5855855/1‑76 (MQ=255)

                                         caaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGacg                               <  1:5648796/76‑1 (MQ=255)

                                         caaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGacg                               >  2:6261575/1‑76 (MQ=255)

                                          aaGGCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGacgc                              >  1:5529335/1‑76 (MQ=255)

                                           aGGCTGGAAAGAGGAGGACCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGacgca                             >  2:5626273/1‑76 (MQ=255)

                                            ggCTGGAAAGAGGAGGATCTAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGacgcac                            <  2:4399747/76‑1 (MQ=255)

                                            ggCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGacgcac                            >  1:236276/1‑76 (MQ=255)

                                            ggCTGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGacgcac                            >  1:5738202/1‑76 (MQ=255)

                                               tGGAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCAcgcg                         >  1:1442959/1‑76 (MQ=255)

                                                 gAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGtt                       >  1:1903430/1‑76 (MQ=255)

                                                 gAAAGAGGAGGATCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGtt                       >  2:1467305/1‑76 (MQ=255)

                                                       ggaggaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTc                 >  2:3615532/1‑76 (MQ=255)

                                                       ggaggaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTc                 <  1:480693/76‑1 (MQ=255)

                                                       ggaggaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTc                 <  1:5168064/76‑1 (MQ=255)

                                                        gaggaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTcg                >  1:6252241/1‑76 (MQ=255)

                                                        gaggaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTcg                >  1:1101621/1‑76 (MQ=255)

                                                        gaggaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTcg                <  1:3565826/76‑1 (MQ=255)

                                                          ggaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTcgcc              <  2:5048592/76‑1 (MQ=255)

                                                           gaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTcgccg             <  1:799637/76‑1 (MQ=255)

                                                           gaTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTcgccg             <  2:5192665/76‑1 (MQ=255)

                                                            aTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTcgccgc            <  2:1753369/76‑1 (MQ=255)

                                                            aTCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTcgccgc            >  1:3095257/1‑76 (MQ=255)

                                                             tCGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCa           <  1:3341343/76‑1 (MQ=255)

                                                              cGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCAt          <  2:3410028/76‑1 (MQ=255)

                                                              cGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCAt          >  1:3462298/1‑76 (MQ=255)

                                                              cGAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCAt          <  1:5706584/76‑1 (MQ=255)

                                                               gAAGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATa         <  2:5922660/76‑1 (MQ=255)

                                                                 aGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACg       <  1:5606890/76‑1 (MQ=255)

                                                                 aGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACg       >  1:4760597/1‑76 (MQ=255)

                                                                 aGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACg       <  1:894625/76‑1 (MQ=255)

                                                                 aGTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACg       >  2:4671873/1‑76 (MQ=255)

                                                                  gTATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACGt      >  2:4097672/1‑76 (MQ=255)

                                                                   tATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACgag     >  1:5911141/1‑74 (MQ=255)

                                                                   tATGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACGTg     >  2:100525/1‑76 (MQ=255)

                                                                    aTGTTTGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACGTGc    >  1:1030935/1‑76 (MQ=255)

CAGCATTTTGCTCAAAAAATAGCCATACTATTTAATTGCAACAAGGCTGGAAAGAGGAGGATCGAAGTATGTTCGTTGACAGACAGCGAATCGATCTGCTGAACCGGTTGATCGACGCACGCGTTGACCTCGCCGCATACGTGCAA  >  REL606/870732‑870877

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