BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,763,636	0	A	G	30.0%	60.0 / 36.7	50	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (19/16); new base G (8/7); total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCA  >  minE/1763568‑1763704

ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAggg                                                                    <  1:622856/71‑1 (MQ=255)

ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                     <  1:704825/70‑1 (MQ=255)

  ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                                  >  1:708752/1‑71 (MQ=255)

  ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                                  >  1:3863326/1‑71 (MQ=255)

   tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc                                                                 >  1:124560/1‑71 (MQ=40)

    ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc                                                                >  1:2766186/1‑71 (MQ=255)

    ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc                                                                >  1:1741501/1‑71 (MQ=255)

        aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                            <  1:2098878/71‑1 (MQ=39)

        aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                            <  1:531098/71‑1 (MQ=39)

         gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct                                                           <  1:3868565/71‑1 (MQ=255)

           cGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct                                                           >  1:578946/1‑69 (MQ=255)

           cGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCttt                                                         >  1:2429107/1‑71 (MQ=255)

           cGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCttt                                                         <  1:3548909/71‑1 (MQ=255)

            gTTACCCCACCAACAAGCTAATCCCATCTGGGTACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg                                                        <  1:3651301/71‑1 (MQ=255)

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

             ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                       <  1:2562394/71‑1 (MQ=39)

               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCttt                                                         <  1:414024/67‑1 (MQ=255)

                ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                    <  1:3431846/71‑1 (MQ=39)

                 cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCtt                                                   <  1:2923103/71‑1 (MQ=255)

                  ccaccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                             >  1:1484477/1‑60 (MQ=255)

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGt                                                      >  1:3650530/1‑66 (MQ=255)

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

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTg                                                  <  1:1239968/70‑1 (MQ=40)

                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                                 >  1:1881231/1‑70 (MQ=255)

                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                <  1:1940678/71‑1 (MQ=39)

                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGa                                               <  1:3457002/71‑1 (MQ=255)

                       aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCttt                                                         <  1:2219522/59‑1 (MQ=255)

                        acaaGCTAATCCCATCTGTGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                            >  1:3443857/1‑71 (MQ=255)

                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                               <  1:2778671/68‑1 (MQ=35)

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

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

                           aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                         <  1:1295563/71‑1 (MQ=34)

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

                               aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCGCCTCTTTGGTCTTGCGACATTATGCgg                                     >  1:3441107/1‑71 (MQ=21)

                               aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                     >  1:2199305/1‑71 (MQ=34)

                                  cccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                  >  1:625199/1‑71 (MQ=34)

                                        tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                            >  1:911784/1‑71 (MQ=34)

                                           gCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCg                         <  1:3000366/71‑1 (MQ=34)

                                              caTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                      >  1:124889/1‑71 (MQ=34)

                                              caTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                      >  1:120180/1‑71 (MQ=34)

                                                tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTcc                    <  1:1543063/71‑1 (MQ=34)

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

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

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

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

                                                   gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                  >  1:1740675/1‑70 (MQ=35)

                                                         aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt           <  1:726670/71‑1 (MQ=34)

                                                          agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                  <  1:2626874/47‑1 (MQ=38)

                                                          agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATc          >  1:2675385/1‑71 (MQ=34)

                                                          agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATc          >  1:2670904/1‑71 (MQ=34)

                                                                  gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa  >  1:1158890/1‑71 (MQ=34)

                                                                   aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAtt                                 >  1:2634037/1‑39 (MQ=40)

                                                                   aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCt     >  1:3516304/1‑67 (MQ=35)

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCA  >  minE/1763568‑1763704

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: `‑`