BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,814,967	A→G	5.0%	R70R (CGT→CGC)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,814,967	0	A	G	5.0%	358.5 / 8.1	119	R70R (CGT→CGC)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (57/56); new base G (3/3); total (60/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 139 total aligned reads displayed.

GCATTTCGAAGAAGGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCAAGGAACACATCCTTGA  >  Exported/1814897‑1815037

gCATTTCGAAGAAGGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgca                                                                        >  1:1049320/1‑71 (MQ=255)

     tCGAAGAAGGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgcacgcac                                                                   >  1:4279085/1‑71 (MQ=255)

     tCGAAGAAGGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgcacgcac                                                                   >  1:3626822/1‑71 (MQ=255)

       gaagaaGGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgcacgcacgc                                                                 <  1:6704483/71‑1 (MQ=255)

        aagaagGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgcacgcacgca                                                                <  1:5914252/71‑1 (MQ=255)

         agaagGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgcacgcacgca                                                                <  1:2968682/70‑1 (MQ=255)

            agGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            >  1:4963637/1‑71 (MQ=255)

                 tGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTggg                                                        <  1:2460648/70‑1 (MQ=255)

                 tGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTggg                                                        <  1:1708652/70‑1 (MQ=255)

                   gTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgcacgc                                                                     <  1:4536589/55‑1 (MQ=255)

                   gTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGaa                                                      <  1:4904720/70‑1 (MQ=255)

                    tGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGt                                                    >  1:6927302/1‑71 (MQ=255)

                    tGACGCGCAGTGTAGCCGACGTTTTCCAGGTCGTTGTGCTTGCCGCCGGCGCGCACGCAGCGCTggacc                                                      >  1:6776399/1‑66 (MQ=255)

                       cgcgcgGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGt                                                 <  1:2410880/71‑1 (MQ=255)

                       cgcgcgGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGt                                                 <  1:6796990/71‑1 (MQ=255)

                       cgcgcgGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGt                                                 <  1:2611455/71‑1 (MQ=255)

                        gcgcgGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGAAGCGTTGGGAAGTgg                                                  >  1:2047235/1‑69 (MQ=255)

                          gcgGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAgc                                              <  1:762287/71‑1 (MQ=255)

                            agtgtAGCCGACGTTTTCCAGGTCGTTGTGCTTGCCGCCGGCGCGCACGCAGc                                                              <  1:4039606/52‑1 (MQ=255)

                            agtgtAGCCGACGTTTTCCAGGTCGTTGTGCTTGCCGCCGGCGCGCACGCAGc                                                              <  1:693214/52‑1 (MQ=255)

                             gtgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGcgcg                                           >  1:2035440/1‑71 (MQ=255)

                             gtgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGcgcg                                           <  1:3153994/71‑1 (MQ=255)

                             gtgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACACGCACGCACGCAGCGTTGGGAAGTGGTAGcgcg                                           >  1:3290043/1‑71 (MQ=255)

                              tgtAGCCGACGTTTTCCAGGTCGTTGTGCTTGCCGCCggcgcgc                                                                     >  1:4612607/1‑44 (MQ=16)

                              tgtAGCCGACGTTTTCCAGGTCGTTGTGCTTGCCGCCGGCGCGCACGCAGCgc                                                            >  1:5933786/1‑52 (MQ=255)

                              tgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            <  1:2569325/53‑1 (MQ=255)

                              tgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            <  1:1440566/53‑1 (MQ=255)

                              tgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            <  1:5084803/53‑1 (MQ=255)

                              tgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            <  1:5104585/53‑1 (MQ=255)

                              tgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            <  1:3609329/53‑1 (MQ=255)

                              tgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            <  1:772218/53‑1 (MQ=255)

                              tgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            <  1:6706853/53‑1 (MQ=255)

                              tgtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGcgcg                                           <  1:1092097/70‑1 (MQ=255)

                               gtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgca                                                                        >  1:5368208/1‑40 (MQ=255)

                               gtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCcgca                                                                        >  1:3535069/1‑40 (MQ=255)

                               gtAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCgg                                          >  1:1956590/1‑70 (MQ=255)

                                tAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTggg                                                        <  1:870724/55‑1 (MQ=255)

                                tAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTggg                                                        <  1:2193344/55‑1 (MQ=255)

                                tAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCggg                                         >  1:6191533/1‑70 (MQ=255)

                                tAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCggg                                         >  1:3939106/1‑70 (MQ=255)

                                tAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCggg                                         >  1:4744540/1‑70 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            >  1:5966846/1‑49 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            >  1:2810246/1‑49 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            >  1:6525830/1‑49 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            >  1:6901468/1‑49 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            >  1:3306027/1‑49 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            >  1:5267482/1‑49 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGt                                                            >  1:2769467/1‑49 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGTGAAGTGGTAGCGCGGGaa                                       >  1:4181749/1‑70 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaat                                      <  1:5255385/71‑1 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaa                                       >  1:2317881/1‑70 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaa                                       <  1:5396176/70‑1 (MQ=255)

                                  aCCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaa                                       <  1:3151176/70‑1 (MQ=255)

                                  aCCGACGTTTTCCAGCTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaa                                       >  1:5959100/1‑70 (MQ=255)

                                   cctaCGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGcgcg                                           <  1:377531/62‑1 (MQ=255)

                                   ccGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaata                                     <  1:5786401/71‑1 (MQ=255)

                                   ccGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaata                                     <  1:2751974/71‑1 (MQ=255)

                                   ccGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaata                                     <  1:1030993/71‑1 (MQ=255)

                                   ccGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaata                                     <  1:6662543/71‑1 (MQ=255)

                                     gACGTTTTCCATGTAGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaataat                                   <  1:2549287/71‑1 (MQ=255)

                                      aCGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaataat                                   >  1:5447127/1‑70 (MQ=255)

                                         ttttCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaat                                      <  1:4441685/64‑1 (MQ=255)

                                         ttttCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACg                               <  1:6309461/71‑1 (MQ=255)

                                         ttttCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACg                               <  1:2849518/71‑1 (MQ=255)

                                          tttCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGc                              >  1:5503595/1‑71 (MQ=255)

                                           ttCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCt                             >  1:5008827/1‑71 (MQ=255)

                                           ttCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCt                             >  1:768313/1‑71 (MQ=255)

                                             ccAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCtt                            <  1:1205094/70‑1 (MQ=255)

                                              cAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTg                           <  1:3861902/70‑1 (MQ=255)

                                              cAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGt                          <  1:4984460/71‑1 (MQ=255)

                                                ggTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGTGAAGTGGTAGCGCGGGAATAATTACGCTTGTCg                        <  1:6218943/71‑1 (MQ=255)

                                                ggTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTc                         >  1:4971532/1‑70 (MQ=255)

                                                ggTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTc                         >  1:2443835/1‑70 (MQ=255)

                                                ggTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTc                         >  1:5320192/1‑70 (MQ=255)

                                                ggTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCg                        <  1:6474001/71‑1 (MQ=255)

                                                ggTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCg                        <  1:948317/71‑1 (MQ=255)

                                                 gTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTc                         >  1:6207465/1‑69 (MQ=255)

                                                 gTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTc                         >  1:1663899/1‑69 (MQ=255)

                                                 gTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTc                         >  1:6329873/1‑69 (MQ=255)

                                                 gTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGa                       <  1:5880571/71‑1 (MQ=255)

                                                  tCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAg                      <  1:5113619/71‑1 (MQ=255)

                                                  tCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCG‑GGGAATAATTACGCTTGTCGAg                      >  1:3609793/1‑70 (MQ=255)

                                                   cGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGc                     >  1:1313101/1‑71 (MQ=255)

                                                     ttGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaataat                                   <  1:5711060/55‑1 (MQ=255)

                                                     ttGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaataat                                   <  1:2561030/55‑1 (MQ=255)

                                                      tgtgTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGca                             >  1:464825/1‑59 (MQ=255)

                                                      tgtgTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGccc                   >  1:3382939/1‑70 (MQ=255)

                                                      tgtgTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGccc                   >  1:2654202/1‑70 (MQ=255)

                                                      tgtgTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGccc                   >  1:2617262/1‑70 (MQ=255)

                                                      tgtgTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGccc                   >  1:5219262/1‑70 (MQ=255)

                                                       gtgtTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGATTAATTACGCTTGTCGAGccc                   >  1:1595110/1‑69 (MQ=255)

                                                        tgtTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCaa                 >  1:6351743/1‑70 (MQ=255)

                                                         gtTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCAAg                >  1:3514400/1‑70 (MQ=255)

                                                            tACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCa                  >  1:4265336/1‑65 (MQ=255)

                                                            tACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCAAGGaa             >  1:3667977/1‑70 (MQ=255)

                                                             accaccCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGaataat                                   <  1:5158313/47‑1 (MQ=255)

                                                                     cacgcacgcaGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCAAGGAACACATCCtt    >  1:4618058/1‑70 (MQ=255)

                                                                     cacgcacgcaGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCAAGGAACACATCCtt    >  1:6869218/1‑70 (MQ=255)

                                                                     cacgcacgcaGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCAAGGAACACATCCtt    >  1:6879887/1‑70 (MQ=255)

                                                                      acgcacgcaGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCAAGGAACACATCCTTGa  <  1:1018359/71‑1 (MQ=255)

GCATTTCGAAGAAGGTATGGTGACGCGCGGTGTAACCGACGTTTTCCAGGTCGTTGTGTTTACCACCCGCACGCACGCAGCGTTGGGAAGTGGTAGCGCGGGAATAATTACGCTTGTCGAGCCCAAGGAACACATCCTTGA  >  Exported/1814897‑1815037

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