BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,145	T→C	12.0%	intergenic (‑10/+53)	argZ ← / ← argY	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,818,145	0	T	C	12.0%	88.3 / 2.0	35	intergenic (‑10/+53)	argZ/argY	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (7/23); new base C (2/2); total (10/25)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTT  >  NC_000913/2818003‑2818286

tGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTAtttt                                                                                                                                         <  1:1245393/150‑1 (MQ=255)

tGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAATTAgttt                                                                                                                                         >  1:446065/1‑150 (MQ=255)

  cGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAc                                                                                                                                       <  1:1281963/150‑1 (MQ=255)

  cGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAc                                                                                                                                       <  2:708747/150‑1 (MQ=255)

    tttATTTCCGCTTGATCTTCGTTACCTTCATTTCTGAATATTTCGGATTATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTGGTAGCCGAGTACTCTATCCAGCTGAGCTACGGTTGGATCGGGAAACTTACTTTActg                                                                                                                                     <  2:446065/149‑1 (MQ=14)

      tAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgct                                                                                                                                   <  2:716490/150‑1 (MQ=255)

                  ctctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGTTGATATTG‑ATAcc                                                                                                                      <  2:850246/150‑1 (MQ=255)

                  ctctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTAATTTACTGCAGAtttttttaaa                                                                                                                       >  2:525760/1‑144 (MQ=11)

                   tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATAcc                                                                                                                      <  1:1445889/149‑1 (MQ=255)

                   tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATAcc                                                                                                                      >  2:1445889/1‑149 (MQ=255)

                     ttCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGct                                                                                                                   <  1:1505097/150‑1 (MQ=255)

                          tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTaa                                                                                                              <  1:849217/150‑1 (MQ=255)

                            ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTGCGACCGCTCGGTTAGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATAGGGATTTTTAGTTTAATGATGATATTG‑TTAGCGCTAATATAg                                                                                                            >  2:586220/1‑150 (MQ=14)

                            ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGctac                                                                                                                 <  1:1478008/145‑1 (MQ=255)

                            ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGctac                                                                                                                 >  2:1478008/1‑145 (MQ=255)

                            ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAg                                                                                                            <  1:886903/150‑1 (MQ=255)

                               tCATTGCTGAAGTTTTCGGATAATATGGTGCATCCGGTAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCa                                                                                                         <  1:1167479/150‑1 (MQ=255)

                               tCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCa                                                                                                         <  2:403567/150‑1 (MQ=255)

                               tCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCa                                                                                                         <  2:28266/150‑1 (MQ=255)

                                 aTTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCata                                                                                                       <  2:1431896/150‑1 (MQ=255)

                                           aTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAaga                                                                                             <  2:1490350/150‑1 (MQ=255)

                                              acGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCTCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa                                                                                                                                 <  1:849823/111‑1 (MQ=255)

                                              acGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa                                                                                                                                 >  2:849823/2‑112 (MQ=40)

                                               cGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACttactttact                                                                                                                                      >  2:137788/1‑106 (MQ=9)

                                               cGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACttactttact                                                                                                                                      <  1:137788/106‑1 (MQ=9)

                                                  agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTg                                                                                      <  2:1180106/150‑1 (MQ=40)

                                                  agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTg                                                                                      <  2:840877/150‑1 (MQ=40)

                                                            gCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGggag                                                                            <  2:111789/150‑1 (MQ=40)

                                                                    gaggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGaa                                                                    <  2:1387512/150‑1 (MQ=255)

                                                                     ctcaatCGAACCTCAGACCTCTGGTCTCGTATCCGATTACTCAATCTATCTGTCCTTCTTATGAAGCGGGATATTTAATTTCCTGTTGATATTGTATA‑CGTTAATTTCTCAATATTAAGTATTATATGGTGCATCCGGGAGGATTCGAAc                                                                   <  2:977427/145‑1 (MQ=255)

                                                                         atCGAACCTCCGACCGCTCGGGGCGCTTCCGCGTACTCCATACAGCAGCGCTACGTATGCAGCGGGTCATTTAATTTACTGCTGATCTTG‑CTTCCGCTACTATAGCCATATCTCGTACTATATGGTGCATCCGGGAGGATTCGAACCTcc                                                               <  2:1184939/149‑1 (MQ=255)

                                                                                                                  ccAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTT‑ATAACGCTAATAAAGCCATATCAAGTAAGAGATGGTGAATCAGGGAGGATTAGAAACGAAGACCGCAGAGATCGTAACCGAGTACTCTATCAAgaggatct                      >  2:1197485/1‑143 (MQ=17)

                                                                                                                                     aTCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                       >  2:542991/1‑130 (MQ=255)

                                                                                                                                     aTCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                       <  1:542991/130‑1 (MQ=255)

                                                                                                                                      tCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTGAGCGCGAGCGCATACTAGCCGAGTACACTGACCCGATGAGCTCCGGATGCCGAGGGAAaatt  >  1:729028/1‑150 (MQ=12)

TGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTG‑ATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTT  >  NC_000913/2818003‑2818286

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