BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	225,003	0	G	A	76.0%	2.0 / 12.2	25	noncoding (1233/1542 nt)	rrsH	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (1/5); new base A (3/16); total (4/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01
Rejected as consensus: E-value score below prediction cutoff.

CCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224867‑225018

ccGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTa              >  1:1150884/1‑140 (MQ=2)

 cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtactac             >  2:372530/1‑140 (MQ=2)

 cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtactac             >  2:106839/1‑140 (MQ=2)

  gCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACa            <  1:1540853/140‑1 (MQ=32)

  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:1211099/140‑1 (MQ=2)

  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:526550/140‑1 (MQ=2)

  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:1078955/140‑1 (MQ=2)

  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:372530/140‑1 (MQ=2)

  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  2:342908/140‑1 (MQ=2)

  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:71439/140‑1 (MQ=2)

   cAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa            <  1:1067416/140‑1 (MQ=2)

   cAACGAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACaa           >  2:394804/1‑140 (MQ=2)

      cGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATgg        <  1:1721801/140‑1 (MQ=32)

      cGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATgg        <  1:1497198/140‑1 (MQ=32)

       gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:870557/140‑1 (MQ=2)

       gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:948868/140‑1 (MQ=2)

       gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:869305/140‑1 (MQ=2)

       gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:977036/140‑1 (MQ=2)

       gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:359469/140‑1 (MQ=2)

       gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:466811/140‑1 (MQ=2)

       gAGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  1:1131151/140‑1 (MQ=2)

        ggcgaAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCGAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGgcg      <  2:1243256/139‑1 (MQ=12)

        aGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:294938/140‑1 (MQ=2)

        aGCGCAACCCCTATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc      >  1:1193388/1‑139 (MQ=2)

            cAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATa  <  2:1410256/140‑1 (MQ=32)

            cAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATa  >  1:1739724/1‑140 (MQ=33)

                ccctATTTTTAGTTGCCAGC‑ACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc       <  2:427183/128‑1 (MQ=2)

CCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGC‑GGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATA  >  NC_000913/224867‑225018

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