BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,558	G→C	10.2%	noncoding (286/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,558	0	G	C	10.2%	505.5 / 39.5	176	noncoding (286/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (76/82); new base C (12/6); total (88/88)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

Only 100 of 183 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAG  >  Exported/1763502‑1763627

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                         >  1:3628801/1‑71 (MQ=32)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                          <  1:1004722/70‑1 (MQ=14)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                         >  1:2319412/1‑71 (MQ=17)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                         >  1:1921084/1‑71 (MQ=17)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                         >  1:1858234/1‑71 (MQ=17)

  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                       <  1:4079601/71‑1 (MQ=34)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                   >  1:3367420/1‑71 (MQ=21)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                   >  1:4731660/1‑71 (MQ=21)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                                    <  1:3369704/68‑1 (MQ=35)

         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                                >  1:2474083/1‑71 (MQ=34)

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                               <  1:428441/71‑1 (MQ=21)

          ccGTGTCTCAGGTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                               <  1:3223644/71‑1 (MQ=21)

           cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                              >  1:3741835/1‑71 (MQ=34)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                            >  1:2500466/1‑71 (MQ=34)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                             <  1:3859730/70‑1 (MQ=25)

                ctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGCCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                         >  1:4619897/1‑71 (MQ=255)

                ctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                         >  1:798499/1‑71 (MQ=21)

                ctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                         >  1:4199549/1‑71 (MQ=21)

                ctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                         >  1:2101466/1‑71 (MQ=21)

                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                          <  1:1898411/53‑1 (MQ=31)

                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                             >  1:2042295/1‑66 (MQ=35)

                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                       >  1:2831627/1‑56 (MQ=255)

                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                             >  1:3139172/1‑66 (MQ=25)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                        >  1:3943943/1‑53 (MQ=35)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                       >  1:2318077/1‑70 (MQ=35)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                         <  1:3678858/68‑1 (MQ=35)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                         <  1:3415938/68‑1 (MQ=35)

                   aGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCcca                                         <  1:4298232/68‑1 (MQ=11)

                    gTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                     >  1:961651/1‑71 (MQ=38)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                    <  1:891497/71‑1 (MQ=39)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                                  <  1:1138459/57‑1 (MQ=255)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                                  <  1:4511650/57‑1 (MQ=255)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                                  <  1:4560797/57‑1 (MQ=255)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                                  <  1:4663587/57‑1 (MQ=255)

                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                    >  1:3579975/1‑70 (MQ=255)

                      tCCAGTGTGGCTAGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                    >  1:2180233/1‑70 (MQ=255)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                    >  1:3785609/1‑69 (MQ=40)

                           tgtgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcctc                                          >  1:974829/1‑57 (MQ=11)

                           tgtgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcctc                                          >  1:3933247/1‑57 (MQ=11)

                           tgtgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcctc                                          >  1:930577/1‑57 (MQ=11)

                           tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                          <  1:2212591/43‑1 (MQ=32)

                           tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                          <  1:3611651/43‑1 (MQ=32)

                           tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                          <  1:716173/43‑1 (MQ=32)

                           tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                          <  1:1535058/43‑1 (MQ=32)

                           tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                              <  1:3040396/71‑1 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                              <  1:561457/71‑1 (MQ=255)

                           tgtgGCTGGGCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                          <  1:1203898/43‑1 (MQ=12)

                            gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                             >  1:3454117/1‑71 (MQ=255)

                            gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                             >  1:1212042/1‑71 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                             <  1:1324711/70‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                             <  1:3389987/70‑1 (MQ=255)

                              ggCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                            <  1:3099215/70‑1 (MQ=255)

                              ggCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                        <  1:3971345/58‑1 (MQ=255)

                              ggCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                        <  1:1067824/58‑1 (MQ=255)

                              ggCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                        <  1:1644756/58‑1 (MQ=255)

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAATca                          <  1:327476/71‑2 (MQ=11)

                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAATca                          <  1:4369799/71‑2 (MQ=11)

                               gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                            >  1:1876064/1‑53 (MQ=255)

                               gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                          <  1:1848359/71‑1 (MQ=255)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                         >  1:4217589/1‑39 (MQ=37)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                         >  1:4061959/1‑39 (MQ=37)

                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                         >  1:3643056/1‑54 (MQ=35)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                         <  1:662157/70‑1 (MQ=255)

                                  ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                       >  1:2265316/1‑71 (MQ=255)

                                   gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                             >  1:2485762/1‑64 (MQ=255)

                                   gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                              <  1:21607/63‑1 (MQ=255)

                                   gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                        <  1:3785622/69‑1 (MQ=255)

                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCAATCTg                    >  1:1240109/1‑71 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                  <  1:4446918/71‑1 (MQ=255)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                            >  1:2706952/1‑43 (MQ=38)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                              >  1:4642159/1‑57 (MQ=255)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                       >  1:1965886/1‑64 (MQ=255)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                       >  1:476005/1‑64 (MQ=255)

                                         ctctCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACcaa                                     >  1:133083/1‑50 (MQ=14)

                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                     >  1:387119/1‑50 (MQ=255)

                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca              >  1:1543872/1‑71 (MQ=255)

                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                   <  1:1699503/49‑1 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                              <  1:3178927/38‑1 (MQ=255)

                                              aGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATccga                        >  1:79743/1‑56 (MQ=11)

                                              aGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATccga                        >  1:252049/1‑56 (MQ=11)

                                              aGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATccga                        >  1:3170688/1‑56 (MQ=11)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                               >  1:2788129/1‑35 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                               >  1:3767065/1‑35 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca                >  1:2292478/1‑66 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt             <  1:3975296/69‑1 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt             <  1:3004847/69‑1 (MQ=255)

                                                  cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                         <  1:2532926/53‑1 (MQ=255)

                                                  cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTA‑CCCACCAACAAGCTAATCCCATCTGGGCACATGCGATgg      <  1:4418/71‑1 (MQ=14)

                                                   aGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg      <  1:1185447/71‑1 (MQ=255)

                                                   aGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg      <  1:177761/71‑1 (MQ=255)

                                                   aGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg      <  1:2387595/71‑1 (MQ=255)

                                                   aGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg      <  1:1303607/71‑1 (MQ=255)

                                                      caGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCACCAAGCTAATCCCATCTGGGCACATCCGATGGCaa   <  1:2434602/70‑1 (MQ=255)

                                                      caGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc     >  1:3312778/2‑69 (MQ=255)

                                                      caGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa   <  1:1202027/70‑1 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa   <  1:4175146/70‑1 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAag  <  1:4435866/71‑1 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAag  <  1:2593443/71‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa   <  1:2837172/69‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAag  <  1:3984355/70‑1 (MQ=255)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAG  >  Exported/1763502‑1763627

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