BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,636	A→G	18.8%	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	18.8%	144.8 / 30.7	69	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (23/33); new base G (5/8); total (28/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTC  >  minE/1763568‑1763702

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

 cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                                 >  1:1762923/1‑71 (MQ=255)

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

   tGGTGAGCCGTTACCCCACCAACAAGCTAATACCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGt                                                                 <  1:133192/69‑1 (MQ=255)

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

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

      tgtgCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc                                                              <  1:2255985/66‑1 (MQ=255)

        aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc                                                          >  1:3909055/1‑71 (MQ=255)

        aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                           <  1:3690607/70‑1 (MQ=255)

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

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

             ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcccc                                                             <  1:3082532/63‑1 (MQ=255)

               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGt                                                    >  1:1076993/1‑70 (MQ=255)

               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                   >  1:88779/1‑71 (MQ=255)

               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                    <  1:2293750/70‑1 (MQ=40)

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

                 cccACCAACAAGCTAATCCCATCTGGGCAGATCGGATGGCAAGTGGGCCGAAGGTCCCCCTCTTTGGTCt                                                  <  1:3748826/70‑1 (MQ=14)

                 cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCtt                                                 <  1:332698/71‑1 (MQ=39)

                 cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                  <  1:3937828/70‑1 (MQ=40)

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                               <  1:1414810/71‑1 (MQ=255)

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                           >  1:951158/1‑59 (MQ=255)

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                      >  1:2119290/1‑64 (MQ=40)

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                               >  1:2238374/1‑71 (MQ=39)

                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcccc                                                             <  1:172782/56‑1 (MQ=255)

                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc                                                              >  1:3764097/1‑54 (MQ=255)

                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                             <  1:260412/71‑1 (MQ=39)

                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCGTTGGTCTTGCGa                                             <  1:1041696/71‑1 (MQ=255)

                     ccaACAAGCTAATCCAATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                             <  1:2399992/71‑1 (MQ=255)

                       aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCtt                                                 <  1:3121980/65‑1 (MQ=255)

                       aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcg                                           >  1:1057873/1‑70 (MQ=19)

                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                    >  1:1334533/1‑45 (MQ=38)

                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcgt                                          <  1:2609696/71‑3 (MQ=9)

                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                          <  1:2496339/71‑1 (MQ=38)

                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                          <  1:2855112/71‑1 (MQ=38)

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

                          aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                      >  1:1389220/1‑57 (MQ=31)

                          aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                           <  1:2240768/68‑1 (MQ=33)

                             cTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                     <  1:3582871/71‑1 (MQ=34)

                                aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                   <  1:106073/70‑1 (MQ=35)

                                aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                                  <  1:2082555/71‑1 (MQ=34)

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

                                   ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                                 >  1:1117293/1‑69 (MQ=35)

                                   ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                >  1:337115/1‑70 (MQ=35)

                                     aTCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAg                             <  1:3254652/71‑1 (MQ=34)

                                      tCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                                  <  1:2737423/65‑1 (MQ=35)

                                       cTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                          <  1:699049/56‑1 (MQ=33)

                                        tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                                 >  1:2146554/1‑64 (MQ=35)

                                          ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                        <  1:351105/55‑1 (MQ=35)

                                          ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCt                           <  1:933922/68‑1 (MQ=35)

                                          ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                        <  1:333473/71‑1 (MQ=34)

                                          cgCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                            <  1:333959/66‑1 (MQ=25)

                                             acaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt                     <  1:3551618/71‑1 (MQ=34)

                                              caTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                          <  1:3133329/49‑1 (MQ=35)

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

                                                tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                                 >  1:756254/1‑56 (MQ=35)

                                                tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                   >  1:3111954/1‑70 (MQ=35)

                                                tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                   >  1:913459/1‑70 (MQ=35)

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

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

                                                    aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagta              <  1:2189513/71‑1 (MQ=34)

                                                     tGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAACGTTTc                   <  1:1744531/65‑1 (MQ=25)

                                                       gCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt            <  1:3407046/70‑1 (MQ=35)

                                                        cAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                   <  1:3138472/46‑1 (MQ=38)

                                                           gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                   <  1:2473538/59‑1 (MQ=35)

                                                           gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt            >  1:3544608/1‑66 (MQ=35)

                                                            agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAg                             >  1:86028/1‑48 (MQ=38)

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

                                                               cccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                   <  1:1388426/55‑1 (MQ=35)

                                                                 cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt               <  1:3858806/57‑1 (MQ=35)

                                                                 cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc  <  1:1313207/70‑1 (MQ=35)

                                                                  gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt           <  1:1005334/60‑1 (MQ=35)

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTC  >  minE/1763568‑1763702

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