BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	980,254	A→G	8.5%	G17G (GGT→GGC)	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	980,254	0	A	G	8.5%	289.3 / 12.4	106	G17G (GGT→GGC)	cspI	cold shock protein
Reads supporting (aligned to +/- strand): ref base A (52/45); new base G (6/3); total (58/48)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 126 total aligned reads displayed.

GCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAA  >  Exported/980188‑980324

gCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCt                                                                     <  1:3456392/70‑1 (MQ=255)

 ctctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                                   >  1:4420352/1‑71 (MQ=255)

   ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                 >  1:104868/1‑71 (MQ=255)

    tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                 >  1:1620991/1‑70 (MQ=255)

      aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                 <  1:3675894/68‑1 (MQ=255)

      aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                 <  1:1508089/68‑1 (MQ=255)

      aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                                <  1:2949205/69‑1 (MQ=255)

      aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                              <  1:4446119/71‑1 (MQ=255)

        ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                               <  1:5582488/68‑1 (MQ=255)

        ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                               <  1:2919994/68‑1 (MQ=255)

        ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                               <  1:4550533/68‑1 (MQ=255)

        ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                               <  1:4578964/68‑1 (MQ=255)

          gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                                       >  1:5680592/1‑58 (MQ=255)

          gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                                           >  1:3022718/1‑70 (MQ=255)

           cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                 <  1:5457685/63‑1 (MQ=255)

           cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                                         <  1:3922659/71‑1 (MQ=255)

           cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                               >  1:2389092/1‑65 (MQ=255)

            tGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                 <  1:3410184/62‑1 (MQ=255)

             gagaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                                       >  1:1216871/1‑71 (MQ=255)

               gaAATGGACAAACTCATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                                       >  1:2788336/1‑53 (MQ=37)

                aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                                    <  1:2167906/71‑1 (MQ=255)

                aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                                      >  1:220450/1‑69 (MQ=255)

                  atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                                  >  1:1039078/4‑71 (MQ=255)

                   tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                                 <  1:1379447/69‑1 (MQ=255)

                     gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                              >  1:2381801/1‑56 (MQ=255)

                     caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttta                                               >  1:1453682/2‑69 (MQ=255)

                     caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACAAtttta                                               >  1:1298922/2‑69 (MQ=255)

                      aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                                <  1:4048922/53‑1 (MQ=255)

                      aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                               <  1:4158835/54‑1 (MQ=255)

                      aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                               <  1:395464/54‑1 (MQ=255)

                        aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                                             >  1:4306035/1‑54 (MQ=255)

                        aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                            <  1:5510559/71‑1 (MQ=255)

                        aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                            <  1:720516/71‑1 (MQ=255)

                        aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                            <  1:2315796/71‑1 (MQ=255)

                        aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                                             >  1:5153911/1‑70 (MQ=255)

                          acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                                            >  1:1052620/1‑67 (MQ=9)

                          acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                                            >  1:4621677/1‑67 (MQ=9)

                           cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                                 >  1:1567424/1‑63 (MQ=255)

                           cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                                >  1:1156252/1‑64 (MQ=255)

                           cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                                         <  1:1513145/71‑1 (MQ=255)

                             caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                              >  1:5185732/1‑48 (MQ=255)

                             caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                              >  1:1659681/1‑48 (MQ=255)

                             caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                              >  1:631642/1‑48 (MQ=255)

                                cTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCa                                       >  1:5511648/1‑68 (MQ=255)

                                 tttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                                      >  1:1263479/1‑36 (MQ=255)

                                 tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCa                                   >  1:3707685/1‑71 (MQ=255)

                                  ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGt                                     >  1:3644293/1‑68 (MQ=255)

                                  ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAt                                  >  1:5221626/1‑71 (MQ=255)

                                   tgctgcCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAc                                         >  1:4022803/1‑63 (MQ=9)

                                    gctgcCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCAtt                                 <  1:2100745/70‑1 (MQ=12)

                                     gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                                <  1:2003423/68‑1 (MQ=255)

                                     gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                                <  1:3892167/68‑1 (MQ=255)

                                     ctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                                <  1:1591862/70‑1 (MQ=255)

                                       gcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                                      <  1:3448645/46‑1 (MQ=255)

                                        ccATCTTTGGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                             >  1:2737445/1‑70 (MQ=255)

                                          aTCTTTCGGTGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                            <  1:3328034/53‑1 (MQ=255)

                                          aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                            <  1:3092487/53‑1 (MQ=255)

                                          aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                            <  1:160521/53‑1 (MQ=255)

                                             tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                       >  1:1425652/1‑71 (MQ=255)

                                              ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                       >  1:4028751/1‑70 (MQ=255)

                                              ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                       >  1:2488971/1‑70 (MQ=255)

                                                cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa                    >  1:151545/1‑71 (MQ=255)

                                                  gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAGCCATTTTACGATACCAGTCAtt                                 >  1:3514997/3‑56 (MQ=255)

                                                  gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                 >  1:4568935/3‑56 (MQ=255)

                                                  gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                 >  1:1857598/3‑56 (MQ=255)

                                                    gTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                 >  1:2103233/1‑70 (MQ=255)

                                                     tgatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCTCTAAACCAGTCATTTTGTTAGACATAatta                 <  1:4387484/69‑1 (MQ=255)

                                                     tgataaAGCCAAAACCTTTATCTGCGTTAAACCAt                                                   <  1:5295674/35‑1 (MQ=12)

                                                       atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                                <  1:3203337/52‑1 (MQ=255)

                                                       atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                             <  1:210327/55‑1 (MQ=255)

                                                       atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                             <  1:8492/55‑1 (MQ=255)

                                                       atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                             <  1:5633836/55‑1 (MQ=255)

                                                       atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                             <  1:4753505/55‑1 (MQ=255)

                                                       atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                             <  1:4986770/55‑1 (MQ=255)

                                                        tgaAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                 >  1:899380/1‑50 (MQ=255)

                                                         gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                              >  1:4645415/1‑52 (MQ=255)

                                                          aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAATCCAGTCATTTTGTTAGACATAATTATTAc             >  1:3841039/1‑68 (MQ=255)

                                                          aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtt          >  1:923052/1‑71 (MQ=255)

                                                           aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCttt         >  1:4390881/1‑71 (MQ=255)

                                                             ccAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg       <  1:2256864/71‑1 (MQ=255)

                                                             ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTc                                    >  1:1237644/1‑42 (MQ=25)

                                                             ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCa                                   >  1:1923790/1‑43 (MQ=25)

                                                             ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCa                                   >  1:3179028/1‑43 (MQ=25)

                                                              cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtttt                               >  1:5024372/1‑46 (MQ=255)

                                                              cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg       >  1:5017903/1‑70 (MQ=255)

                                                               aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                             >  1:5375720/1‑47 (MQ=255)

                                                               aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg       >  1:490479/1‑69 (MQ=255)

                                                               aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg       >  1:3169106/1‑69 (MQ=255)

                                                               aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg       >  1:2074788/1‑69 (MQ=255)

                                                               aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg       >  1:5448399/1‑69 (MQ=255)

                                                               aaaaCCTTTATCTGCGTTAAACCATTTTACTAAACCAGt                                     <  1:4281685/39‑1 (MQ=25)

                                                                aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                           >  1:2629414/1‑48 (MQ=255)

                                                                 aaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaaga   <  1:405913/71‑1 (MQ=255)

                                                                 aaCCTTTATCAGCGTTAAACCATTTTACTAAACCa                                       >  1:1012348/1‑35 (MQ=25)

                                                                  aCCTTTTTCAGGGTTAAACCATTTCACTATACCAGTCATTTTGt                             <  1:5152605/44‑1 (MQ=255)

                                                                  aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                >  1:2418489/1‑57 (MQ=255)

                                                                  aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                >  1:4752905/1‑57 (MQ=255)

                                                                  aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                >  1:4194263/1‑57 (MQ=255)

                                                                  aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa  <  1:5581304/71‑1 (MQ=255)

                                                                  aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa  <  1:2837519/71‑1 (MQ=255)

GCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAA  >  Exported/980188‑980324

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