BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,636	A→G	27.3%	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	27.3%	85.2 / 40.4	55	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (15/25); new base G (6/9); total (21/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAT  >  minE/1763567‑1763705

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

 ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAgg                                                                     <  1:735488/71‑1 (MQ=255)

 ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAgg                                                                     <  1:1562994/71‑1 (MQ=255)

 ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAgg                                                                     <  1:3594988/71‑1 (MQ=255)

 ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAgg                                                                     <  1:2363364/71‑1 (MQ=255)

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

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

     ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                   <  1:2487938/69‑1 (MQ=255)

      gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATCGCAAGAGGCCCGAAGGTccc                                                                 <  1:1125099/70‑1 (MQ=255)

        gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                              >  1:3396555/1‑71 (MQ=255)

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

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

          gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctct                                                            <  1:3961615/71‑1 (MQ=39)

           ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct                                                            <  1:4008810/70‑1 (MQ=255)

                 ccccACCAACAAGCTAATCCCATCTGGGGACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                      <  1:382335/70‑1 (MQ=14)

                 ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgggc                                                      <  1:387385/70‑3 (MQ=255)

                 ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                      <  1:3967093/70‑1 (MQ=255)

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

                     accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                  >  1:2737418/1‑70 (MQ=40)

                     accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                  >  1:3358908/1‑70 (MQ=40)

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

                       caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAc                                               >  1:2898841/1‑71 (MQ=255)

                       caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAc                                               >  1:4113424/1‑71 (MQ=255)

                       caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                >  1:505316/1‑70 (MQ=39)

                       caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                >  1:2796496/1‑70 (MQ=39)

                       caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                >  1:4049284/1‑70 (MQ=39)

                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                        >  1:1734690/1‑61 (MQ=37)

                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                 >  1:835715/1‑68 (MQ=37)

                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              >  1:2804635/1‑71 (MQ=38)

                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              >  1:3624577/1‑71 (MQ=38)

                        aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              <  1:364991/71‑1 (MQ=38)

                         acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg                                                 <  1:1236592/67‑1 (MQ=25)

                         acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcgt                                             >  1:873090/1‑69 (MQ=9)

                           aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                           >  1:3741120/1‑71 (MQ=35)

                                 aTCCCATCTGGGTACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:3845165/70‑1 (MQ=25)

                                    ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              <  1:3841098/59‑1 (MQ=31)

                                     cATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                                 <  1:1869310/71‑1 (MQ=34)

                                     cATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                                 <  1:3199562/71‑1 (MQ=34)

                                     cATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                                 <  1:2565597/71‑1 (MQ=34)

                                     cATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                                 <  1:752409/71‑1 (MQ=34)

                                           ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                               <  1:2111610/67‑1 (MQ=35)

                                           ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                               <  1:3208288/67‑1 (MQ=35)

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

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

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

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

                                                   cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                      >  1:2569805/1‑68 (MQ=35)

                                                    gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt                  <  1:2505170/71‑1 (MQ=34)

                                                     aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt                  <  1:455525/70‑1 (MQ=35)

                                                     aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                   >  1:1150216/1‑69 (MQ=35)

                                                     aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                   >  1:992166/1‑69 (MQ=35)

                                                     aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCa                    >  1:45847/1‑68 (MQ=35)

                                                            gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:997725/43‑1 (MQ=38)

                                                             agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:3945700/42‑1 (MQ=38)

                                                             agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc         <  1:4096762/71‑1 (MQ=34)

                                                                    aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt  <  1:346915/71‑1 (MQ=34)

                                                                    aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt  <  1:2616592/71‑1 (MQ=34)

GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAT  >  minE/1763567‑1763705

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