BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,173,321	A→T	33.2%	intergenic (‑12/+3)	trnB‑Leu2 ← / ← trnB‑Gly1	tRNA‑Leu/tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,173,321	0	A	T	33.2%	17.8 / 2.6	29	intergenic (‑12/+3)	trnB‑Leu2/trnB‑Gly1	tRNA‑Leu/tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base A (6/12); new base T (3/7); total (9/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTCTATCCAATTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGTATAGTTGGAGCGGAAGACGGGATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCACGTCTGTAAAGACACTAGAGCCTG  >  BS168_reference_genome/3173185‑3173456

aaaCGCTCTATCCAATTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATTGTATAGt                                                                                                                                       <  6:147564/139‑1 (MQ=17)

aaaCGCTCTATCCAATTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGTATAGt                                                                                                                                       >  5:312645/1‑139 (MQ=32)

  aCGCTCTATCCAATTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGTATAGTTg                                                                                                                                     <  6:312645/139‑1 (MQ=32)

      cagaTCCAATTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGTCTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGGGGATGGTATAGTTGGAg                                                                                                                                  <  5:232080/135‑1 (MQ=12)

        tATCCAATTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGGTGGTTTTGTTGGAGCgg                                                                                                                               >  7:394548/1‑139 (MQ=12)

            cAATTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGTATAGTTGGAGCGGAAGa                                                                                                                           <  5:216483/139‑1 (MQ=32)

              aTTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGTATAGTTGGAGCGGAAGACg                                                                                                                         <  3:100412/139‑1 (MQ=32)

                  aGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGGTGGTATTGTTGGAGCGGAAGACgggct                                                                                                                     >  6:157557/1‑137 (MQ=12)

                  aGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATTGTATTGTTGGGCGGGAAGACGGGAt                                                                                                                     >  3:118131/1‑139 (MQ=11)

                  aGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGTATAGTTGGAGCGGAAGACGGGAt                                                                                                                     >  3:326301/1‑139 (MQ=32)

                       cGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGGATATTTGGAGCGGAAGACGGGATTctca                                                                                                                >  2:367962/1‑136 (MQ=12)

                                aaCTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGTATAGTTGGAGCGGAAGACGGGATTCGAACCCGCGAcc                                                                                                       <  4:326301/139‑1 (MQ=32)

                                                        gacccGGCACGGTAGGCACCTACCGCAGGATTTTACGTCCTGTGTGGCTGCCAATACCACCACCCCGGCGGGGTTGGTATAGTGGGAGCGGAAGACGGGATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTc                                                                                <  7:148847/135‑1 (MQ=1)

                                                                                   agATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGGGGATGGTATAGTTGGAGCGGAAGACGGGATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATgg                                                    <  1:333937/138‑1 (MQ=18)

                                                                                    taTTTTAATTCCTGTGTGTCTCCCACTTCCACCCCCCCGGCGTGGTTGTTATTGTTGGAGCGCAACACGGAATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGTACTACTTCCGCAGAAATGGt                                                   <  1:173847/138‑1 (MQ=11)

                                                                                      ttttATGTCCGGTGTGTGTGACAATCCCCCCCCCCCGGGGTGGATGGTTTTGTGGGCGCGGAAGACGGGATGCGACCCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTg                                                  <  6:110494/138‑1 (MQ=11)

                                                                                      ttttAGGTCCTGTGTGTCTGCCAATGCAACCACCCCGGCGTGGTTGGTATTGTTGTCGCGGAACCGGGTATTCTAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTg                                                  <  4:347904/138‑1 (MQ=11)

                                                                                                  tgtgtCTGCCAATTCCACCACCCCGGCGTGGCTGGTTTAGTTGGGGCGGAAGACCGGATTCGCACCCCCCCCCCCCCCCCTGGCAAGGTGGTGTTGGAGCACAGGACTACTTCAGCACAAAAGGTGGGGCTGAAggggc                                     >  1:296656/1‑137 (MQ=2)

                                                                                                    tgtCGGCAAATTCCACCCCCCCGGGGTGTTTGGTATTGTTGGAGCGGAAGACGGTATTCCAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACtt                                   <  6:133612/139‑1 (MQ=17)

                                                                                                           cAATACCACCACCCCGGGGGGTATGGTATAGTGGGCGCGGACGCCGGGATTCGACCCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAAcccc                            <  5:204565/139‑1 (MQ=17)

                                                                                                           aaaTTCCACCCCCCCGGCGGGTATGTTTTTGTGGCAGCGCAAGAGGGGATTCGACCCCGCGACCCCCACCTTGGCAAGGGGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAAcccc                            <  5:388375/138‑1 (MQ=17)

                                                                                                            aaTTCCACCACCCCGGCGTGGATGGTATATTTGGAGCGGAAGACGGGATTAGAACCCCCGACCCCCCCCATGGCAAGGTGGTGTTGTACCACAGAACTACTTACGCAGAAATTGTGCGGATGAAGGGACTTGAAccccc                           >  4:386610/1‑139 (MQ=17)

                                                                                                             aTTCCCCCCCCCCGGGGTGTATGGTATAGTTGGAGCGGACGACGGGATGCGACCCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCa                          <  5:96358/139‑1 (MQ=17)

                                                                                                                    caccCCGGCGTGTATTGAATAGTTGGAGCGGAAGACGGGATACGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTGCTACCACTGAACTACTACCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCACGTCTGt                   <  7:76700/139‑1 (MQ=18)

                                                                                                                    caccCCGGCGGGGATGTTATAGTTGGAGCGGACGACGGGATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCACGTCTGt                   <  1:3769/139‑1 (MQ=255)

                                                                                                                     accCGGGCGTGTATGGTATCGTTGCAGCGTACCCGGGGATTCCACCCCGCGACCCCCACCTTGGCAAGGTGGTGGTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCACGTCTGTa                  <  1:284612/139‑1 (MQ=18)

                                                                                                                      ccccGGGGGGGATGGTATTGTGGTCGCGGAAGACGGGATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCACGTCTGTaa                 <  2:267053/139‑1 (MQ=18)

                                                                                                                                ggtGGTATAGTTGGAGCGGGAGACGGGATTCGAACCCGCCACCCCCACCGGGGCAAGGTGGTGTTGTACCACGGAACTTCTTCCGCAGAAATGGTGCGGATGATGGGACGTGATCCCCCCCCTCAGTAAAGacacga         >  6:179947/3‑135 (MQ=17)

                                                                                                                                 aTGGTATTGTTGGAGCGGAAGACGGGATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCACGTCTGTAAAGACACTagag      <  3:386610/139‑1 (MQ=255)

                                                                                                                                     tataGTTGGAGCGGAAGACGGGATTCGCACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGTGACTTGAACCCCCACGTCTGTAAAGACACTAGAGCCTg  >  4:94824/1‑139 (MQ=255)

                                                                                                                                     tataGTTGGAGCGGAAGACGGGATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCACGTCTGTAAAGACACTAGAGCCt   >  8:75672/1‑138 (MQ=255)

AAACGCTCTATCCAATTGAGCTACGGGCGCAAAACTTAATGCCGAGGGCCGGACTTGAACCGGCACGGTAGTCACCTACCGCAGGATTTTAAGTCCTGTGTGTCTGCCAATTCCACCACCCCGGCGTGGATGGTATAGTTGGAGCGGAAGACGGGATTCGAACCCGCGACCCCCACCTTGGCAAGGTGGTGTTCTACCACTGAACTACTTCCGCAGAAATGGTGCGGATGAAGGGACTTGAACCCCCACGTCTGTAAAGACACTAGAGCCTG  >  BS168_reference_genome/3173185‑3173456

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