BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,813,392	G→A	15.1%	R595R (CGC→CGT)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,813,392	0	G	A	15.1%	316.4 / 41.0	119	R595R (CGC→CGT)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base G (47/54); new base A (3/15); total (50/69)
Fisher's exact test for biased strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 125 total aligned reads displayed.

TGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTGGCGCAGCGCAG  >  Exported/1813325‑1813462

tgtTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAgcgc                                                                     <  1:3210150/71‑1 (MQ=255)

tgtTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAgcg                                                                      <  1:1937903/70‑1 (MQ=255)

  ttCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAgccc                                                                     >  1:4459594/1‑67 (MQ=255)

  ttCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCgcag                                                                   >  1:3026768/1‑71 (MQ=255)

  ttCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCgca                                                                    <  1:2435897/70‑1 (MQ=255)

    cACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCgcagca                                                                 >  1:4338319/1‑71 (MQ=255)

    cACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCgcagca                                                                 >  1:2303315/1‑71 (MQ=255)

    cACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCgcagca                                                                 >  1:1559329/1‑71 (MQ=255)

      ccAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCAcc                                                               >  1:3081440/1‑71 (MQ=255)

      ccAGGTATTCGCCCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCAcc                                                               >  1:4474561/1‑71 (MQ=255)

       cAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCt                                                              >  1:3598104/1‑71 (MQ=255)

        aGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCtt                                                             >  1:824002/1‑71 (MQ=255)

        aGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCtt                                                             >  1:4756175/1‑71 (MQ=255)

          gTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGt                                                           <  1:1475550/71‑1 (MQ=255)

            cTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCg                                                         >  1:2757230/1‑71 (MQ=255)

               cGACCGCACTAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt                                                       <  1:812238/70‑1 (MQ=255)

               cGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt                                                       <  1:2175638/70‑1 (MQ=255)

               cGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt                                                       <  1:2961565/70‑1 (MQ=255)

               cGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt                                                       <  1:1990196/70‑1 (MQ=255)

               cGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt                                                       >  1:4132456/1‑70 (MQ=255)

               cGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAACGCAGCACCTTGTCGtt                                                       >  1:532210/1‑70 (MQ=255)

               cGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAACGCAGCACCTTGTCGtt                                                       >  1:600447/1‑70 (MQ=255)

                gACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTaa                                                     >  1:3015679/1‑71 (MQ=255)

                gACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTa                                                      <  1:1794054/70‑1 (MQ=255)

                   cGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAAc                                                    >  1:4488538/1‑69 (MQ=255)

                   cGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCa                                                  <  1:2166767/71‑1 (MQ=255)

                     cACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAACGCAGCACCTTGTCGTTAACCAGt                                                <  1:1477983/71‑1 (MQ=255)

                      aCGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGt                                                >  1:2092612/1‑70 (MQ=255)

                      aCGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTg                                               <  1:3469626/71‑1 (MQ=255)

                       cGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt                                                       <  1:1800769/62‑1 (MQ=255)

                       cGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt                                                       <  1:2672546/62‑1 (MQ=255)

                       cGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGtt                                                       <  1:2576698/62‑1 (MQ=255)

                       cGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAACGCAGCACCTTGTCGTTAACCAGTGa                                              <  1:1644884/71‑1 (MQ=255)

                          aTCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCgcagc                                                                  <  1:3883721/48‑1 (MQ=255)

                           tctcTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCtt                                                             <  1:845037/52‑1 (MQ=255)

                           tctcTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAAc                                                    <  1:4735949/61‑1 (MQ=255)

                             tcttctGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCttt                                        >  1:1520715/1‑71 (MQ=255)

                             tcttctGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCtt                                         >  1:1390709/1‑70 (MQ=255)

                               ttctGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCtt                                         <  1:2000100/68‑1 (MQ=255)

                               ttctGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCt                                          <  1:2785610/67‑1 (MQ=255)

                               ttctGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCt                                          <  1:4068867/67‑1 (MQ=255)

                                tctGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGt                                                           <  1:3410541/49‑1 (MQ=255)

                                tctGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGt                                                           <  1:926613/49‑1 (MQ=255)

                                 ctGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTc                                       <  1:3658253/68‑1 (MQ=255)

                                  tGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCg                                                         >  1:1564717/1‑49 (MQ=255)

                                    gTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGc                                    <  1:4238235/68‑1 (MQ=255)

                                    gTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGc                                    <  1:3512760/68‑1 (MQ=255)

                                    gTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGc                                    <  1:2063305/68‑1 (MQ=255)

                                    gTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGa                                  <  1:2452525/70‑1 (MQ=255)

                                     tttCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCtt                                         <  1:4722236/62‑1 (MQ=255)

                                      ttCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCt                                                              <  1:2114539/40‑1 (MQ=255)

                                      ttCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCt                                                              <  1:1934049/40‑1 (MQ=255)

                                      ttCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATa                                >  1:3674155/1‑70 (MQ=255)

                                      ttCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATa                                >  1:1166470/1‑70 (MQ=255)

                                         aTCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACAt                             <  1:3698422/70‑1 (MQ=255)

                                         aTCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACAt                             <  1:2837704/70‑1 (MQ=255)

                                          tCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCa                                                  >  1:3805191/1‑48 (MQ=255)

                                          tCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCa                                                  >  1:3549655/1‑48 (MQ=255)

                                          tCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGa                           <  1:892333/71‑1 (MQ=255)

                                            gCTTCGTTGTGTGAGAAGTCGAAGCGCAGCTCCTTGTCGTTaa                                                     >  1:1129752/1‑43 (MQ=255)

                                            gCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAAc                                                    >  1:3428535/1‑44 (MQ=255)

                                            gCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAAc                                                    >  1:1184079/1‑44 (MQ=255)

                                            gCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAAc                                                    >  1:2321691/1‑44 (MQ=255)

                                            gCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAAc                                                    >  1:2435606/1‑44 (MQ=255)

                                            gCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAg                          >  1:3140089/1‑70 (MQ=255)

                                              ttCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTAc                       >  1:880791/1‑71 (MQ=255)

                                                cGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTAcc                      >  1:554723/1‑70 (MQ=255)

                                                    gtgtgAGAAGTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACGCAGa                  <  1:672959/70‑1 (MQ=255)

                                                    gtgtgAGAAGTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGa                  <  1:3791972/70‑1 (MQ=255)

                                                    gtgtgAGAAGTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGa                  <  1:2840640/70‑1 (MQ=255)

                                                    gtgtgAGAAGTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGa                  <  1:2744659/70‑1 (MQ=255)

                                                         agaATTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTgg            <  1:2186236/71‑1 (MQ=255)

                                                         agaAGTCGAAGCGCAGCACCTTGTCGTTAATCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTgg            <  1:1772808/71‑1 (MQ=255)

                                                         agaAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAAGCTTTCTGCGATACATGAGTACCCAGAACCTgg            <  1:649527/71‑1 (MQ=255)

                                                         agaAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTc                                       <  1:4687636/44‑1 (MQ=255)

                                                         agaAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTgg            <  1:3864708/71‑1 (MQ=255)

                                                         agaAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTgg            <  1:283548/71‑1 (MQ=255)

                                                         agaAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTgg            <  1:2498579/71‑1 (MQ=255)

                                                         agaAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTgg            <  1:3767203/71‑1 (MQ=255)

                                                           aaGTCGAAGCGCAGCACCTTGTCGTTAACAAGTGAACCTTTCTg                                     >  1:3402192/1‑44 (MQ=255)

                                                             gTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTc                                       >  1:686048/1‑40 (MQ=255)

                                                             gTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATa                                <  1:565822/47‑1 (MQ=255)

                                                             gTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATa                                <  1:1116878/47‑1 (MQ=255)

                                                             gTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATa                                <  1:4011387/47‑1 (MQ=255)

                                                             gTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATa                                <  1:2506274/47‑1 (MQ=255)

                                                             gTCGAAACGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATa                                <  1:4292671/47‑1 (MQ=255)

                                                              tCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAAcc               <  1:578525/63‑1 (MQ=255)

                                                              tCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAAcc               <  1:3099868/63‑1 (MQ=255)

                                                              tCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAAcc               <  1:1488319/63‑1 (MQ=255)

                                                              tCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAAcc               <  1:3572457/63‑1 (MQ=255)

                                                               cGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACAt                             >  1:1906490/1‑48 (MQ=255)

                                                               cGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACAt                             >  1:4398002/1‑48 (MQ=255)

                                                               cGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTg             >  1:4283098/1‑64 (MQ=255)

                                                               cGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTg             >  1:697239/1‑64 (MQ=255)

                                                               cGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTg             >  1:2830033/1‑64 (MQ=255)

                                                                gAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTGgcgcagcg     >  1:475682/1‑71 (MQ=255)

                                                                gAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTGgcgcagcg     >  1:2988175/1‑71 (MQ=255)

                                                                  aGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTGgcgcagcgc    <  1:1880616/70‑1 (MQ=255)

                                                                   gcgcAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTGGcgcagcgcag  <  1:778867/71‑1 (MQ=255)

                                                                   gcgcAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTGGcgcagcgcag  <  1:1751621/71‑1 (MQ=255)

TGTTCACCAGGTCTTCGACCGCACGAATCTCTTCTGGTTTCATCGCTTCGTTGTGTGAGAAGTCGAAGCGCAGCACCTTGTCGTTAACCAGTGAACCTTTCTGCGATACATGAGTACCCAGAACCTGGCGCAGCGCAG  >  Exported/1813325‑1813462

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