BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,232,364	C→A	10.2%	G131V (GGA→GTA)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,232,364	0	C	A	10.2%	220.4 / 7.9	68	G131V (GGA→GTA)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base C (31/30); new base A (3/4); total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCATCCCA  >  Exported/1232295‑1232434

ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                                       <  1:3652393/71‑1 (MQ=255)

ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                                       <  1:1742057/71‑1 (MQ=255)

  gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                     >  1:242790/1‑71 (MQ=255)

  gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAAAGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                     >  1:2515542/1‑71 (MQ=255)

             gggAAAGAAAATCTTCTGGCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaa                                                          <  1:2173542/71‑1 (MQ=255)

             gggAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaa                                                          <  1:1189974/71‑1 (MQ=255)

               gaaagaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaaaa                                                        <  1:1435691/71‑1 (MQ=255)

                  agaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGc                                                                     >  1:1766811/1‑55 (MQ=255)

                   gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    >  1:546977/1‑71 (MQ=255)

                   gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    <  1:1310724/71‑1 (MQ=255)

                   gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    <  1:2323041/71‑1 (MQ=255)

                   gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    >  1:1637556/1‑71 (MQ=255)

                    aaaaTCTTCTGCCAGCGACATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                                   >  1:2319817/1‑71 (MQ=255)

                     aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                                   >  1:2356153/1‑70 (MQ=255)

                     aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                                   <  1:1746360/70‑1 (MQ=255)

                     aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                                  <  1:1673387/71‑1 (MQ=255)

                     aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                                  <  1:14212/71‑1 (MQ=255)

                        tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                                   >  1:2025436/1‑51 (MQ=21)

                        tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                                   >  1:3135235/1‑51 (MQ=21)

                        tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                                   >  1:1842334/1‑51 (MQ=21)

                         cttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTAt                                                >  1:107123/1‑69 (MQ=255)

                           tctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGc                                            >  1:1840121/1‑71 (MQ=255)

                            ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                               <  1:2257400/67‑1 (MQ=255)

                            ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                               <  1:1772999/67‑1 (MQ=255)

                            ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                               <  1:43657/67‑1 (MQ=255)

                                cAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTgcgc                                       <  1:1553584/71‑1 (MQ=255)

                                  gcgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                                    <  1:3221466/56‑1 (MQ=255)

                                   cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                                    >  1:329111/1‑71 (MQ=255)

                                   cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                                    >  1:355919/1‑71 (MQ=255)

                                    gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                                   <  1:2452894/71‑1 (MQ=255)

                                    gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                                   <  1:3243234/71‑1 (MQ=255)

                                     ccATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTgcgcgc                                     <  1:2305181/68‑1 (MQ=255)

                                     ccATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTa                                  >  1:2349938/1‑71 (MQ=255)

                                      cATAATGGTGTTATCGGGCGTAATCGCCATTCCGCACATCGAATAAAACTCACTATGAGCTGCGCGCTTAt                                 <  1:3102388/71‑1 (MQ=255)

                                       aTAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTAtt                                >  1:664892/1‑71 (MQ=255)

                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTa                                                 <  1:957888/53‑1 (MQ=255)

                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:3282933/1‑71 (MQ=255)

                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:3287211/1‑71 (MQ=255)

                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:1004755/1‑71 (MQ=255)

                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:1857135/1‑71 (MQ=255)

                                        tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                               >  1:1891808/1‑71 (MQ=255)

                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGGGCGCTTATTgg                              >  1:2246104/1‑71 (MQ=255)

                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:472771/1‑71 (MQ=255)

                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:3656216/1‑71 (MQ=255)

                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:3619948/1‑71 (MQ=255)

                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:2117673/1‑71 (MQ=255)

                                         aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                              >  1:3205421/1‑71 (MQ=255)

                                                  aTCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGa                        >  1:1101784/1‑68 (MQ=255)

                                                   tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa                    <  1:2859813/71‑1 (MQ=255)

                                                   tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa                    <  1:1338074/71‑1 (MQ=255)

                                                    cGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa                    <  1:470352/70‑1 (MQ=255)

                                                     gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                            <  1:1042718/45‑1 (MQ=18)

                                                     gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                            <  1:113402/45‑1 (MQ=18)

                                                     gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                            <  1:3324596/45‑1 (MQ=18)

                                                     gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                            <  1:1601209/45‑1 (MQ=18)

                                                        cGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGa                                              <  1:1999930/40‑1 (MQ=255)

                                                         gTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTc              <  1:2189975/71‑1 (MQ=255)

                                                         gTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTc              >  1:1735470/1‑71 (MQ=255)

                                                          tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGc                                            >  1:1805924/1‑40 (MQ=255)

                                                          tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc             <  1:2043933/71‑1 (MQ=255)

                                                          tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc             >  1:238452/1‑71 (MQ=255)

                                                             tCGCCATTCCGCCCATCGAATAAAACTCACTATGa                                              >  1:2978945/1‑35 (MQ=255)

                                                               gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGc                      <  1:3155017/57‑1 (MQ=255)

                                                               gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGc                      <  1:1708764/57‑1 (MQ=255)

                                                               gCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTc        <  1:2489568/71‑1 (MQ=255)

                                                                ccATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTc        >  1:1117393/1‑70 (MQ=255)

                                                                ccATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCa       >  1:3390240/1‑71 (MQ=255)

                                                                   ttCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTa                                  <  1:1467919/41‑1 (MQ=255)

                                                                   ttCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTa                                  <  1:952355/41‑1 (MQ=255)

                                                                   ttCCGCCCATCGAATAAAACTCACTATGAGCTGCCCGCTTa                                  <  1:2792/41‑1 (MQ=255)

                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTAtt                                >  1:601354/1‑41 (MQ=255)

                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAAt                   >  1:2418162/1‑54 (MQ=255)

                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAAt                   >  1:2815644/1‑54 (MQ=255)

                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAAt                   >  1:3458154/1‑54 (MQ=255)

                                                                     ccGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCATCCCa  >  1:2107557/1‑71 (MQ=255)

CCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCCCACTCATCCCA  >  Exported/1232295‑1232434

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