BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	980,269	A→G	10.4%	F12F (TTT→TTC)	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	980,269	0	A	G	10.4%	303.6 / 20.2	115	F12F (TTT→TTC)	cspI	cold shock protein
Reads supporting (aligned to +/- strand): ref base A (51/52); new base G (6/6); total (57/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCT  >  minE/980199‑980333

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

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

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

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

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

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

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

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

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

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

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

               acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                                                     >  1:1052620/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)

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

                      tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                                               <  1:5616151/52‑1 (MQ=255)

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

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

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

                         gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                                          <  1:4479250/70‑1 (MQ=255)

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

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

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

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

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

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

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

                                  tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                                    >  1:5261761/1‑67 (MQ=255)

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

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

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

                                    tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                               >  1:2333624/1‑70 (MQ=255)

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

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

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

                                       gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                          >  1:4568935/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)

                                           gatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                          <  1:3978299/68‑1 (MQ=255)

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

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

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

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

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

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

                                             tgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                            >  1:2128480/1‑64 (MQ=255)

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

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

                                              gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt                           <  1:4098546/64‑1 (MQ=255)

                                              gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt                           <  1:1845520/64‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                    aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACAAGTCATTTTGt                                      >  1:2076230/1‑47 (MQ=255)

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

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

                                                     aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:3166722/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                        ccTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa          >  1:691397/1‑71 (MQ=255)

                                                         cTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAt         <  1:2764800/71‑1 (MQ=255)

                                                          tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:5566213/38‑1 (MQ=12)

                                                          tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:2054396/38‑1 (MQ=12)

                                                          tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:4874612/38‑1 (MQ=12)

                                                          tttatCAGCGTTAAACCATTTTACTAAACGAGTCAttt                                         <  1:447067/38‑1 (MQ=21)

                                                          tttatCAGCGTTAAACCATTTTACTAAACCAGTCAttt                                         <  1:5506022/38‑1 (MQ=25)

                                                          tttatCAGCGTTAAACCATTTTACTAAACCAGTCAttt                                         <  1:2365632/38‑1 (MQ=25)

                                                          tttatCAGCGTTAAACCATTTTACTAAACCAGTCATTTTATTTGa                                  <  1:1796779/45‑1 (MQ=255)

                                                            tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAtt        <  1:2627508/69‑1 (MQ=255)

                                                            tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTa       <  1:165492/70‑1 (MQ=255)

                                                            tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTa       <  1:5659966/70‑1 (MQ=255)

                                                            tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAg      >  1:697425/1‑71 (MQ=255)

                                                              tCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc     <  1:1376651/70‑1 (MQ=255)

                                                              tCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc     <  1:2869713/70‑1 (MQ=255)

                                                               cAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc     <  1:778949/69‑1 (MQ=255)

                                                                aGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCt  <  1:1923509/71‑1 (MQ=255)

                                                                    ttAAACCATTTTACTAAACCAGTCATTTTATTTGa                                  <  1:5405127/35‑1 (MQ=25)

                                                                    ttAAACCATTTTACTAAACCAGTCATTTTATTAGACAt                               <  1:3734753/38‑1 (MQ=38)

                                                                    ttAAACCATTTTACTAAACCAGTCATTTTATTAGACAt                               <  1:26667/38‑1 (MQ=38)

                                                                    ttAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                            <  1:2705755/41‑1 (MQ=255)

                                                                      aaaCCATTTTACTAAACCAGTCATTTTATTAGACATat                             >  1:208707/1‑37 (MQ=37)

                                                                      aaaCCATTTTACTAAACCAGTCATTTTATTAGACATAtttta                         >  1:4194764/1‑40 (MQ=37)

CTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCT  >  minE/980199‑980333

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: `‑`