BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	980,208	A→G	80.5%	Y33H (TAT→CAT)	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	980,208	0	A	G	80.5%	248.9 / 48.3	118	Y33H (TAT→CAT)	cspI	cold shock protein
Reads supporting (aligned to +/- strand): ref base A (8/15); new base G (41/54); total (49/69)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 119 total aligned reads displayed.

AATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCA  >  Exported/980140‑980274

aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATa                                                                    <  1:4530235/69‑1 (MQ=255)

  tACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGa                                                                  >  1:1178598/1‑69 (MQ=255)

  tACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACa                                                                <  1:2791768/71‑1 (MQ=255)

    ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACaaa                                                              <  1:1904003/71‑1 (MQ=255)

    ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACaaa                                                              <  1:4883414/71‑1 (MQ=255)

     cAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAac                                                             <  1:4409717/71‑1 (MQ=255)

       aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACaa                                                               <  1:1403887/67‑1 (MQ=255)

       aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACaa                                                               <  1:4254926/67‑1 (MQ=255)

       aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAaca                                                            >  1:121020/1‑70 (MQ=255)

       aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAaca                                                            >  1:5433044/1‑70 (MQ=255)

       aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAaca                                                            >  1:4185147/1‑70 (MQ=255)

       aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaaa                                                              >  1:4621400/1‑68 (MQ=255)

        aTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGATGCTCTGAATTGCTGAGAAATAGACAAAcaca                                                          >  1:1047789/1‑71 (MQ=255)

         ttCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACAt                                                         <  1:5535246/71‑1 (MQ=255)

             aCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCttt                                                     <  1:4068873/71‑1 (MQ=255)

             aCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCttt                                                     <  1:943845/71‑1 (MQ=255)

              cTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtg                                                    <  1:1760973/71‑1 (MQ=255)

                tCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGAc                                                                 <  1:520180/56‑1 (MQ=255)

                tCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgct                                                  >  1:1177808/1‑71 (MQ=255)

                tCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgct                                                  >  1:1130899/1‑71 (MQ=255)

                 ccTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTTCTGAGAAATGGACAAACACAt                                                         <  1:1430526/63‑1 (MQ=255)

                     aTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGa                                                                  >  1:2254090/1‑50 (MQ=255)

                     aTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGa                                                                  >  1:3612358/1‑50 (MQ=255)

                           aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCttt                                         <  1:3185473/69‑1 (MQ=255)

                           aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                         <  1:785838/69‑1 (MQ=255)

                           aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCtt                                          <  1:4769465/68‑1 (MQ=255)

                           aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCtt                                          <  1:1350852/68‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCtt                                                      <  1:1667448/54‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCtt                                                      <  1:1969854/54‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCtt                                                      <  1:1762082/54‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCAt                                             <  1:3933525/63‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCAt                                             <  1:3939995/63‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCAt                                             <  1:3485250/63‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCAt                                             <  1:1768278/63‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCAt                                             <  1:2025014/63‑1 (MQ=255)

                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                     <  1:1635304/71‑1 (MQ=255)

                              tAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCg                                    >  1:3901107/1‑71 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgc                                                   >  1:2417013/1‑55 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgc                                                   >  1:1380491/1‑55 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgc                                                   >  1:4136755/1‑55 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgc                                                   >  1:1883774/1‑55 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgc                                                   >  1:3680699/1‑55 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtg                                                    >  1:940751/1‑54 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtg                                                    >  1:5287761/1‑54 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtg                                                    >  1:1532988/1‑54 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCt                                           >  1:1827069/1‑63 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCt                                           >  1:1162587/1‑63 (MQ=255)

                               aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                   >  1:2682005/1‑71 (MQ=255)

                                aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                   <  1:4822840/70‑1 (MQ=255)

                                  gTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCtt                                                      <  1:4983481/49‑1 (MQ=255)

                                  gTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCtt                                                      <  1:1201370/49‑1 (MQ=255)

                                   ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAaca                                                            >  1:792594/1‑42 (MQ=255)

                                   ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                <  1:5266302/70‑1 (MQ=255)

                                   ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                <  1:2937666/70‑1 (MQ=255)

                                   ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                <  1:2931050/70‑1 (MQ=255)

                                   ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                <  1:805262/70‑1 (MQ=255)

                                   ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                     <  1:1655403/65‑1 (MQ=255)

                                    tttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAAcaca                                                          <  1:2691464/43‑1 (MQ=255)

                                    tttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAAcaca                                                          <  1:4183817/43‑1 (MQ=255)

                                     ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                              <  1:2261872/70‑1 (MQ=255)

                                     ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                              <  1:3551699/70‑1 (MQ=255)

                                         aaTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                         <  1:5691882/71‑1 (MQ=255)

                                          aTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCttt                                                     >  1:2178736/1‑42 (MQ=255)

                                          aTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                         >  1:4960457/1‑70 (MQ=255)

                                          aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                        <  1:4837128/71‑1 (MQ=255)

                                             gTTGCTCTGAATTGCTGAGAAATGGACAAACACATCt                                                       >  1:3358172/1‑37 (MQ=255)

                                               tGCTCTGAATTGCTGAGAAATGGACAAACACATCttt                                                     <  1:3266698/37‑1 (MQ=255)

                                               tGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                         >  1:5331891/1‑49 (MQ=255)

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

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

                                                   ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                 >  1:2937930/1‑53 (MQ=255)

                                                   ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                 >  1:4610993/1‑53 (MQ=255)

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

                                                    tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTc                                        <  1:3741904/45‑1 (MQ=255)

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

                                                    tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCATGATGAAACCaaa                       <  1:5561120/62‑1 (MQ=255)

                                                    tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                   <  1:5022575/50‑1 (MQ=255)

                                                     gAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                       >  1:4012266/1‑61 (MQ=255)

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

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

                                                      aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                         >  1:3194578/1‑42 (MQ=255)

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

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

                                                        ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatg                               >  1:2411473/1‑50 (MQ=255)

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

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

                                                        ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa             <  1:2919994/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)

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

                                                           cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt               <  1:380875/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)

                                                               gaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                   <  1:4193400/39‑1 (MQ=255)

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

                                                                aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCg                                    >  1:4403533/1‑37 (MQ=255)

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

AATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCA  >  Exported/980140‑980274

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