BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	218,387	A→C	5.5%	L130R (CTG→CGG)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	218,387	0	A	C	5.5%	287.4 / 11.1	109	L130R (CTG→CGG)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (58/45); new base C (3/3); total (61/48)
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.17e-01

Only 100 of 112 total aligned reads displayed.

CGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCC  >  AP009048/218324‑218457

cGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAg                                                                 >  1:3024779/1‑71 (MQ=255)

cGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAg                                                                 >  1:4019046/1‑71 (MQ=255)

cGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAg                                                                 >  1:7845199/1‑71 (MQ=255)

cGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAg                                                                 >  1:5210570/1‑71 (MQ=255)

   aTGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTg                                                              >  1:1455681/1‑71 (MQ=255)

    tGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGt                                                             <  1:4025211/71‑1 (MQ=255)

     gACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGtt                                                            >  1:5690573/1‑71 (MQ=255)

     gACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGtt                                                            >  1:692820/1‑71 (MQ=255)

     gACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGt                                                             >  1:6040030/1‑70 (MQ=255)

       cGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGTCAGCTGTTTg                                                          <  1:4850872/71‑1 (MQ=255)

       cGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTg                                                          >  1:6008168/1‑71 (MQ=255)

        gCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGt                                                         <  1:5154705/71‑1 (MQ=255)

        gCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGt                                                         <  1:7289623/71‑1 (MQ=255)

         ccGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGt                                                         >  1:2217620/1‑70 (MQ=255)

         ccGAAACGCGGACGCACTTCGTCGAGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGt                                                         >  1:5081512/1‑70 (MQ=255)

          cGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCATCTGTTTGTaa                                                       <  1:7920691/71‑1 (MQ=255)

          cGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTaa                                                       <  1:5016081/71‑1 (MQ=255)

            aaaCGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTaa                                                       >  1:582428/1‑69 (MQ=255)

            aaaCGAGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTaa                                                       >  1:107697/1‑69 (MQ=255)

               cgcgGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAg                                                                 >  1:3390775/1‑56 (MQ=255)

                gcgGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGt                                                             <  1:928186/59‑1 (MQ=255)

                gcgGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAag                                                      <  1:6857852/66‑1 (MQ=255)

                gcgGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCt                                                  <  1:1863428/70‑1 (MQ=255)

                  ggACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCa                                                                        >  1:2952442/1‑46 (MQ=255)

                   gACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCt                                                  <  1:6419768/67‑1 (MQ=255)

                     cGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTc                                            >  1:2779051/1‑71 (MQ=255)

                     cGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTc                                            <  1:6434358/71‑1 (MQ=255)

                      gCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCg                                           >  1:5661595/1‑71 (MQ=255)

                      gCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCg                                           >  1:1414550/1‑71 (MQ=255)

                      gCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCg                                           >  1:5851033/1‑71 (MQ=255)

                      gCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCg                                           <  1:4225714/71‑1 (MQ=255)

                        aCTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGt                                          >  1:6964618/1‑70 (MQ=255)

                         cTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGt                                                         <  1:6011183/54‑1 (MQ=255)

                         cTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTa                                        >  1:2592613/1‑71 (MQ=255)

                           tcgtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCa                                                                        >  1:2147863/1‑37 (MQ=255)

                           tcgtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAg                                                                       >  1:7757865/1‑38 (MQ=255)

                           tcgtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAg                                                                       >  1:6117969/1‑38 (MQ=255)

                           tcgtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAagag                                                    >  1:4461647/1‑57 (MQ=255)

                           tcgtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTAc                                       >  1:4125926/1‑70 (MQ=255)

                           tcgtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTAc                                       >  1:4116563/1‑70 (MQ=255)

                            cgtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGt                                          >  1:7905706/1‑66 (MQ=255)

                            cgtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACg                                      >  1:3223692/1‑70 (MQ=255)

                             gtcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGaa                                    <  1:3680193/71‑1 (MQ=255)

                                gcgGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGa                                     <  1:3086800/67‑1 (MQ=255)

                                gcgGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGa                                     <  1:2072540/67‑1 (MQ=255)

                                 cgGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAag                                                      <  1:677432/49‑1 (MQ=255)

                                  ggAACTTGGTCTGTATCTGATAGAAGTTGATCGGCAGTTGCTTGTAGCTGTTCAGCTCGTTGCgcgc                                   <  1:2288271/67‑4 (MQ=255)

                                  ggAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCt                                                               <  1:5537175/39‑1 (MQ=255)

                                  ggAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAAt                                   >  1:3314243/1‑67 (MQ=255)

                                    aaCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGt                                                         >  1:6294556/1‑43 (MQ=255)

                                    aaCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGt                                                         >  1:6680076/1‑43 (MQ=255)

                                    aaCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAagag                                                    >  1:1418989/1‑48 (MQ=255)

                                    aaCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAagag                                                    >  1:7208817/1‑48 (MQ=255)

                                    aaCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAagag                                                    >  1:4133012/1‑48 (MQ=255)

                                    aaCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAagag                                                    >  1:361109/1‑48 (MQ=255)

                                    aaCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAAt                                   >  1:3293076/1‑65 (MQ=255)

                                    aaCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAAt                                   >  1:6696357/1‑65 (MQ=255)

                                     aCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGt                              >  1:1999299/1‑69 (MQ=255)

                                     aCTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCa                            <  1:6410613/71‑1 (MQ=255)

                                        tGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTaa                                                       >  1:3849387/1‑41 (MQ=255)

                                        tGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTa                                                        >  1:4823084/1‑40 (MQ=255)

                                        tGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGt                          >  1:2058648/1‑70 (MQ=255)

                                          gTCTGGATCTGATAGAAGTTCAGCGTCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTGAGTGAt                       <  1:10405/71‑1 (MQ=255)

                                               gATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAAt                                   <  1:3296964/54‑1 (MQ=255)

                                               gATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAAct                   <  1:1517594/70‑1 (MQ=255)

                                                aTCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAActtc                 >  1:6110765/1‑71 (MQ=255)

                                                 tCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTa                                                 <  1:5261824/38‑1 (MQ=255)

                                                    gATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATaa                     >  1:3385821/1‑63 (MQ=255)

                                                    gATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATaa                     >  1:831206/1‑63 (MQ=255)

                                                    gATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATaa                     >  1:2785976/1‑63 (MQ=255)

                                                     aTAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAActtct                <  1:4666057/67‑1 (MQ=255)

                                                        gAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGa          <  1:6085443/70‑1 (MQ=255)

                                                        gAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGa          <  1:6288792/70‑1 (MQ=255)

                                                        gAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGa          <  1:1898471/70‑1 (MQ=255)

                                                        caaGTTCCGCGACAGCTGTTTGTAAGAGCTAAGCTCGTTACGa                                     <  1:1596678/42‑1 (MQ=255)

                                                        caaGTTCCGCGACAGCTGTTTGTAAGAGCTAAGCTCGTTACGa                                     <  1:710519/42‑1 (MQ=255)

                                                        caaGTTCCGCGACAGCTGTTTGTAAGAGCTAAGCTCGTTACGa                                     <  1:2662424/42‑1 (MQ=255)

                                                         aaGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCTTTACGAATCAGGTCAGTGATAACTTCTTCATGAGt        <  1:4451670/71‑1 (MQ=255)

                                                         aaGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGt        >  1:321359/1‑71 (MQ=255)

                                                          aGTTCCGCGACAGCTGTTTGTAAGAGCTAAGCTCGt                                          >  1:2747630/1‑36 (MQ=255)

                                                          aGTTCCGCGACAGCTGTTTGTAAGAGCTAAGCTCGTTACGAAt                                   >  1:27411/1‑43 (MQ=255)

                                                          aGTTCCGCGACAGCTGTTTGTAAGAGCTAAGCTCGTTACGAAt                                   >  1:3379736/1‑43 (MQ=255)

                                                          aGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCATTGATAACTTCTTCATGa          <  1:351332/68‑1 (MQ=255)

                                                          aGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGtt       <  1:7677700/71‑1 (MQ=255)

                                                          aGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGtt       <  1:1182500/71‑1 (MQ=255)

                                                          aGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGtt       <  1:5371099/71‑1 (MQ=255)

                                                           gTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTa                                        >  1:4080378/1‑37 (MQ=255)

                                                           gTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTa                                        >  1:1905171/1‑37 (MQ=255)

                                                           gTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTAc                                       <  1:6512577/38‑1 (MQ=255)

                                                           gTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTAc                                       <  1:7515506/38‑1 (MQ=255)

                                                           gTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTAc                                       <  1:5058852/38‑1 (MQ=255)

                                                           gTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTg      <  1:7725902/71‑1 (MQ=255)

                                                           gTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTg      >  1:3021830/1‑71 (MQ=255)

                                                           gTTCAGAGGCAGCTGTTTGTAAGAGCTAAGCTCGTTa                                        >  1:6526426/1‑37 (MQ=255)

                                                            ttCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACg                                      <  1:3349966/38‑1 (MQ=255)

                                                            ttCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACg                                      <  1:5368019/38‑1 (MQ=255)

                                                            ttCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACg                                      <  1:1608344/38‑1 (MQ=255)

                                                            ttCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACg                                      <  1:1786252/38‑1 (MQ=255)

                                                             tCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTggg    >  1:5322607/1‑71 (MQ=255)

                                                               aGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGcc  <  1:3705983/71‑1 (MQ=255)

CGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCC  >  AP009048/218324‑218457

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