BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	980,265	G→C	27.2%	P14A (CCT→GCT) ‡	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	980,265	0	G	C	27.2%	192.7 / 78.1	114	P14A (CCT→GCT) ‡	cspI	cold shock protein
Reads supporting (aligned to +/- strand): ref base G (44/39); new base C (16/15); total (60/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

Only 100 of 115 total aligned reads displayed.

GCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCT  >  Exported/980198‑980333

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

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

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

   gagaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                                                >  1:1216871/1‑71 (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)

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

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

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

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

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

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

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

                 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)

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

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

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

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

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

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

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

                           gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:82126/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)

                              ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      <  1:4351269/70‑1 (MQ=255)

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

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

                                aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                     <  1:160521/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: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)

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

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

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

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

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

                                             atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                      <  1:8492/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)

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

                                                aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtt                   >  1:923052/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)

                                                    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:2074788/1‑69 (MQ=255)

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

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

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

                                                     aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:490479/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)

                                                      aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                >  1:3407416/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:4752905/1‑57 (MQ=255)

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

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

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

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

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

                                                           tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                         <  1:5566213/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:5659966/70‑1 (MQ=255)

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

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

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

                                                               tCAGCGTTAAACCATTTTACTAAACCAGTCATTTTAttt                                    <  1:3098439/39‑2 (MQ=25)

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

GCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCT  >  Exported/980198‑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: `‑`