BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,558	G→C	8.4%	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	8.4%	460.6 / 25.6	154	noncoding (286/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (63/78); new base C (11/2); total (74/80)
Fisher's exact test for biased strand distribution p-value = 7.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

Only 100 of 164 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGA  >  Exported/1763502‑1763628

ggAGTCTGGCCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                          >  1:4306071/1‑71 (MQ=255)

ggAGTCTGGACCGTGTCTCAGTTCCGGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                           <  1:3374459/70‑1 (MQ=14)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                           <  1:4213758/70‑1 (MQ=31)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                           <  1:2765117/70‑1 (MQ=31)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                           <  1:3661689/70‑1 (MQ=31)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                          <  1:25147/71‑1 (MQ=32)

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                          <  1:3872712/71‑1 (MQ=17)

 gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                         <  1:993772/71‑1 (MQ=21)

 gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                         <  1:164891/71‑1 (MQ=21)

 gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                         <  1:3796229/71‑1 (MQ=21)

   gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                        >  1:4418622/1‑70 (MQ=35)

     cTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                                     >  1:22422/1‑71 (MQ=34)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGc                                                    >  1:1453281/1‑71 (MQ=11)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGc                                                    >  1:3978263/1‑71 (MQ=11)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                           >  1:2472343/1‑64 (MQ=31)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                                                  >  1:2192201/1‑71 (MQ=34)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                     <  1:2937736/68‑1 (MQ=25)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                     <  1:2394198/68‑1 (MQ=25)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                                   <  1:3486400/70‑1 (MQ=25)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCg                                                  <  1:765621/71‑1 (MQ=21)

         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                   <  1:443892/69‑1 (MQ=35)

         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                                 <  1:1205157/71‑1 (MQ=34)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                <  1:2656238/68‑1 (MQ=35)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                <  1:2569415/68‑1 (MQ=35)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                <  1:1599223/68‑1 (MQ=35)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                <  1:1580753/68‑1 (MQ=35)

             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                <  1:1325082/68‑1 (MQ=35)

                ctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                                 >  1:3422176/1‑64 (MQ=35)

                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGcctt                                                 >  1:3972107/1‑61 (MQ=11)

                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGcctt                                                 >  1:1658493/1‑61 (MQ=11)

                 tcAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTgg                                                         <  1:4404606/55‑1 (MQ=11)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCTTCGCCTTGGTGAGCCGTTACCccacc                                        <  1:2844254/70‑1 (MQ=25)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCTTCGCCTTGGTGAGCCGTTACCccacc                                        <  1:2839255/70‑1 (MQ=25)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                        <  1:4277247/70‑1 (MQ=35)

                    gTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGa                                                      >  1:2046916/1‑55 (MQ=11)

                    gTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGa                                                      >  1:3508419/1‑55 (MQ=11)

                    gTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGa                                                      >  1:4323930/1‑55 (MQ=11)

                    gTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                      <  1:3046462/71‑1 (MQ=38)

                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                     >  1:1126538/1‑70 (MQ=40)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                    <  1:3736093/70‑1 (MQ=255)

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

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                    <  1:1814882/70‑1 (MQ=255)

                         aGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                           <  1:2713607/45‑1 (MQ=32)

                           tgtgGCTGGTCATCCTCTCAGCCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                >  1:1410691/1‑70 (MQ=255)

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

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCg                                                  >  1:786689/1‑52 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                >  1:326147/1‑70 (MQ=255)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                >  1:2387843/1‑70 (MQ=255)

                                cTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAATcag                          <  1:195338/71‑3 (MQ=11)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                          <  1:3468487/71‑1 (MQ=255)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                          <  1:3410354/71‑1 (MQ=255)

                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcca                          <  1:3615707/71‑2 (MQ=255)

                                  ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                        >  1:553642/1‑71 (MQ=255)

                                  ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                        >  1:1688533/1‑71 (MQ=255)

                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                          <  1:4347705/36‑1 (MQ=37)

                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                                     <  1:2128821/41‑1 (MQ=38)

                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                                     <  1:4053808/41‑1 (MQ=38)

                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                     >  1:1521361/1‑71 (MQ=255)

                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                     >  1:248158/1‑71 (MQ=255)

                                     cATCCTCTCAGACCAGCTAGGGATAGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                     >  1:3401955/1‑71 (MQ=255)

                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                                  >  1:4360695/1‑58 (MQ=255)

                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                                  >  1:2008808/1‑58 (MQ=255)

                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                                  >  1:3214754/1‑58 (MQ=255)

                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                                  >  1:644805/1‑58 (MQ=255)

                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                <  1:1035113/60‑1 (MQ=255)

                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                    >  1:1403740/1‑71 (MQ=255)

                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                    >  1:3924749/1‑71 (MQ=255)

                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                    >  1:1899084/1‑71 (MQ=255)

                                      aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                <  1:1886495/59‑1 (MQ=255)

                                       tCCTCTCAGACCAGCCAGGGATCGTCTCCTTGGTGAGCCGTTACCCCACcaa                                      <  1:3884132/52‑1 (MQ=14)

                                       tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                      <  1:3668713/52‑1 (MQ=255)

                                       tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                      <  1:278968/52‑1 (MQ=255)

                                       tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                   <  1:1791628/71‑1 (MQ=255)

                                       tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                   <  1:2977841/71‑1 (MQ=255)

                                        ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                   >  1:3541776/1‑70 (MQ=255)

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

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

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                          >  1:1736484/1‑46 (MQ=38)

                                          tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                                 >  1:1562084/1‑54 (MQ=255)

                                          tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                              >  1:924709/1‑41 (MQ=255)

                                          tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                              >  1:3056509/1‑41 (MQ=255)

                                          tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                              >  1:4334207/1‑41 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt              >  1:2669363/1‑71 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt              >  1:4121739/1‑71 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc            <  1:2076955/71‑1 (MQ=255)

                                                aCCAGCTAGGGATCGTCGCCTTGGTTAGCCGTTACCCCACcaacaa                                   <  1:1232134/46‑1 (MQ=255)

                                                    gccaGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg       >  1:3751758/4‑70 (MQ=255)

                                                     ccaGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                   <  1:3605262/55‑1 (MQ=255)

                                                     ccaGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCa     <  1:4171305/69‑1 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                    <  1:1943773/38‑1 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                    <  1:4426729/38‑1 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc      >  1:3237754/1‑68 (MQ=255)

                                                       aGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc      >  1:1751260/1‑68 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc      >  1:3496072/1‑67 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc      >  1:2045981/1‑67 (MQ=255)

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

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

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCa     <  1:2893360/68‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCa     <  1:2818554/68‑1 (MQ=255)

                                                        gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAaga  <  1:3569173/71‑1 (MQ=255)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGA  >  Exported/1763502‑1763628

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