BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,047	T→A	11.9%	intergenic (‑189/+10)	argQ ← / ← argZ	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,818,047	0	T	A	11.9%	105.9 / 7.3	42	intergenic (‑189/+10)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (20/16); new base A (3/2); total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

CTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAG  >  NC_000913/2817899‑2818192

ctaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGAttt                                                                                                                                                   <  1:314134/149‑1 (MQ=255)

        atatCAATTAAGGGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAat                                                                                                                                           <  1:555080/149‑1 (MQ=255)

           tCAATTAAGAGATGGTGCATCCGGGGGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                        <  1:123925/149‑1 (MQ=255)

           tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGGTCGCCCTTAGGGCCGTTGCGGGAGAAACGGTATCATCCCGAGTGGCTGCGGTTAGTTCCGCTACCTTTTCGATACCTTCATTGCTGAAGATTTAGGCGGATATg                                                                                                                                        >  2:627647/1‑149 (MQ=17)

            cAATTACGAGGTGTCGAATCCGGGAGTATTACCCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATgg                                                                                                                                       <  1:84313/149‑1 (MQ=18)

              aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCTTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAAATTTCGGAGAATATGGTg                                                                                                                                     >  2:395212/1‑149 (MQ=255)

              aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTAGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTg                                                                                                                                     >  1:476434/1‑149 (MQ=255)

                tAAGAGATGGTGCATCCGGGAGGATTACTCGCCTCGCGCCCCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTGGTGCCGGCGGTTAATGCCGCTACCGCCCCGATACCTTCATTGCAGAAGAGTTCGGAGAACATgggcca                                                                                                                                   >  2:462535/1‑145 (MQ=14)

                        ggTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGggag                                                                                                                           >  2:664494/1‑149 (MQ=255)

                             aTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTc                                                                                                                      <  1:607493/149‑1 (MQ=255)

                                cGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTGCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGAATTCGaa                                                                                                                   >  2:311441/1‑149 (MQ=255)

                                       aTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                            >  2:361728/1‑149 (MQ=255)

                                            tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGGGAGAATGGTGCATCCGGGAGGATTCGAACCCCCGACCGCt                                                                                                       >  1:548590/1‑149 (MQ=255)

                                               ccTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTGGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCgg                                                                                                    <  1:166142/149‑1 (MQ=255)

                                                cTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGt                                                                                                   <  2:580777/149‑1 (MQ=255)

                                                   gCGGCTCCCCTTTCGGGCCGTTGCTAAAGCAACGTTCTCACCCATTCTGCCCGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCg                                                                                                <  1:240076/149‑1 (MQ=17)

                                                      gCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                                             >  2:559531/1‑149 (MQ=255)

                                                                     cGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTCCTCTATcc                                                                              >  1:260297/1‑149 (MQ=255)

                                                                       tttcTAAAGCAACGTTATCATCCCTAGTGCTTGCGACATATTCCGTTGCCGCTGCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCCGGAGGATTCGAACCTCCGACTGCTCGGTTCGTAGCCGAGTACTCTATcccg                                                                            <  2:515084/146‑3 (MQ=17)

                                                                       ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCCCGGTTCGTAGCCGAGTACTCTATCCAg                                                                            >  2:321331/1‑149 (MQ=255)

                                                                           tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                                        >  2:263701/1‑149 (MQ=255)

                                                                           tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGGTCGTAGCCGAGTACTCTATCCACCtac                                                                        >  2:623500/1‑147 (MQ=255)

                                                                                  aCGTTATGTTACCGAGTGCGTCCGGTTAATTCCGCTACCTCTTCCATATCTTCATTGCTGAAGATTTCGGAGAATAGGGTGCACCCGGGAGGCTTCGAACCTCCGTCCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACgg                                                                 <  2:184366/149‑1 (MQ=14)

                                                                                  aCGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGGACCCCCTACCGCCCTTTTCGAACCCGAATACTCTATCCTGCTGAGcgccgg                                                                 >  1:698454/1‑144 (MQ=17)

                                                                                   cGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTGTATCCAGCTGGGCTACgcc                                                                >  2:348109/1‑147 (MQ=255)

                                                                                           tcCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                                                        >  1:590747/1‑149 (MQ=255)

                                                                                                    ccTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttat                                               >  1:728135/1‑149 (MQ=255)

                                                                                                     cTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGGACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttaat                                              >  2:108279/1‑147 (MQ=255)

                                                                                                                     tACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAGCGGGAAATTTATTTTACTGCTGATATTGa                              <  1:446432/149‑1 (MQ=255)

                                                                                                                        ctctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGAATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAc                           <  1:321331/149‑1 (MQ=255)

                                                                                                                         tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTCATTTTACTGCTGATATTAACAcc                          >  1:310562/1‑149 (MQ=255)

                                                                                                                             cGATACCTTCATTGCTGACGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                      <  1:32243/149‑1 (MQ=255)

                                                                                                                                 aCCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTa                                           <  1:265447/124‑1 (MQ=2)

                                                                                                                                 aCCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTa                                           >  2:265447/1‑124 (MQ=2)

                                                                                                                                  ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                 <  1:263701/149‑1 (MQ=255)

                                                                                                                                   cTTCATTGCTGAATGTTTCGTAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg                <  1:593282/149‑1 (MQ=255)

                                                                                                                                   cTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAGCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg                >  1:115315/1‑149 (MQ=255)

                                                                                                                                     tCATTGCTGAAGGTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgctg                                     <  1:164211/126‑1 (MQ=255)

                                                                                                                                     tCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgctg                                     <  1:164299/126‑1 (MQ=255)

                                                                                                                                     tCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgctg                                     >  2:164299/1‑126 (MQ=255)

                                                                                                                                     tCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgctg                                     >  2:164211/1‑126 (MQ=16)

                                                                                                                                          gCTGAAGATGTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGCCCGCTCGGTTCGTAGCCGAGTACTCGATCCAGCTGAGCTACGGTTGCATCTGGAAATTTGTTTTTCTGCTGATTTTGATACCGCTACTAAAGCCATATc         <  2:380894/149‑1 (MQ=14)

                                                                                                                                                 aTTGCGGAGTATATGGTGCATCCGGGAGGCTTCGGACCTCCGTCCGCTCGGTTTGTAGCCGACTACCCAAACCTGCTGGCCTACGGACGCAGTGGGAAATTTATTTGACTGCCGCAATTGATACCGCTACTAAAGCCATATCAAGTAag  <  2:229080/149‑1 (MQ=9)

CTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAG  >  NC_000913/2817899‑2818192

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