BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,636	A→G	17.6%	noncoding (208/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,763,636	0	A	G	17.6%	113.9 / 19.7	51	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (29/13); new base G (4/5); total (33/18)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCC  >  Exported/1763567‑1763703

gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                    <  1:1053618/71‑1 (MQ=255)

gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                    >  1:4255611/1‑71 (MQ=255)

   ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                                  <  1:596083/70‑1 (MQ=255)

    tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc                                                                >  1:516459/1‑71 (MQ=255)

    tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc                                                                >  1:4708828/1‑71 (MQ=255)

     ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                               >  1:630173/1‑71 (MQ=39)

     ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                               >  1:245894/1‑71 (MQ=39)

     ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                               >  1:4703849/1‑71 (MQ=39)

     ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                               >  1:2346465/1‑71 (MQ=39)

     ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                               >  1:2290360/1‑71 (MQ=39)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                            >  1:493510/1‑70 (MQ=40)

            cGTTACCCCACCAACAAGCTAATCTCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt                                                         <  1:264689/70‑1 (MQ=255)

             gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg                                                       <  1:4394891/71‑1 (MQ=255)

             gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt                                                        >  1:2466649/1‑70 (MQ=40)

             gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                       >  1:1945801/1‑71 (MQ=39)

              ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                      <  1:3529859/71‑1 (MQ=255)

              ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                      >  1:3055949/1‑71 (MQ=255)

                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc                                                    >  1:469849/1‑71 (MQ=39)

                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc                                                    >  1:4468803/1‑71 (MQ=39)

                 ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt                                                   <  1:3615093/71‑1 (MQ=255)

                 ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt                                                   <  1:1335174/71‑1 (MQ=255)

                 ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc                                                    <  1:1540704/70‑1 (MQ=40)

                  cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt                                                        >  1:3792273/1‑65 (MQ=40)

                    caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                                >  1:2515560/1‑71 (MQ=255)

                    caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                <  1:3596735/71‑1 (MQ=39)

                      ccaACAAGCTAATCCCATCTGGGCACATCAGATGGCAAGAGGCCCGAAGGTCCCCctct                                                          <  1:3667579/59‑1 (MQ=255)

                         acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCtt                                                  <  1:2742813/64‑1 (MQ=35)

                          caaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                           >  1:4151963/1‑70 (MQ=37)

                              cTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                         >  1:4485221/1‑68 (MQ=35)

                              cTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                         >  1:2018780/1‑68 (MQ=35)

                               tAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                      <  1:1928594/70‑1 (MQ=35)

                                           ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                         >  1:340388/1‑71 (MQ=34)

                                                aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                    <  1:3812205/71‑1 (MQ=34)

                                                 tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                         >  1:3357235/1‑65 (MQ=35)

                                                  ccGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                    >  1:1185001/1‑69 (MQ=35)

                                                   cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                 <  1:378356/71‑1 (MQ=34)

                                                   cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                 >  1:2535436/1‑71 (MQ=34)

                                                      tGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                         >  1:2959417/1‑60 (MQ=35)

                                                      tGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                         >  1:1655281/1‑60 (MQ=35)

                                                      tGGCAAGACGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag              >  1:1863891/1‑71 (MQ=21)

                                                        gCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt            <  1:1848356/71‑1 (MQ=34)

                                                          aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                            >  1:4246139/1‑37 (MQ=33)

                                                          aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                            >  1:2231027/1‑37 (MQ=33)

                                                           agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAg                              >  1:3972735/1‑50 (MQ=38)

                                                           agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                           >  1:523583/1‑53 (MQ=35)

                                                           agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                           >  1:1965757/1‑53 (MQ=35)

                                                            gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcc        >  1:3107810/1‑71 (MQ=34)

                                                             agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc       >  1:2492939/1‑71 (MQ=34)

                                                              ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc      <  1:3757237/71‑1 (MQ=34)

                                                                  cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTcc  <  1:1207742/71‑1 (MQ=34)

                                                                   gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                    <  1:2818165/36‑1 (MQ=40)

                                                                    aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                                  >  1:1324300/1‑37 (MQ=40)

GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCC  >  Exported/1763567‑1763703

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