BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,556	T→C	43.8%	noncoding (288/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,556	0	T	C	43.8%	38.4 / 204.4	146	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (43/39); new base C (28/36); total (71/75)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 146 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATG  >  Exported/1763502‑1763622

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

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:4210717/71‑1 (MQ=21)

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

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

  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                  >  1:310359/1‑71 (MQ=34)

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

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

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

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

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

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                              >  1:4163493/1‑71 (MQ=34)

      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                              >  1:2039585/1‑71 (MQ=34)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                             >  1:1669673/1‑70 (MQ=35)

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

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

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

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

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

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

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

            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                        <  1:4408300/71‑1 (MQ=21)

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

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

               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                     <  1:3672187/55‑1 (MQ=9)

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

                ctcAGTTCCAGTGTGGCTGGCCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                    >  1:156804/1‑71 (MQ=21)

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

                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGcctt                                           >  1:135568/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)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                  >  1:989805/1‑70 (MQ=25)

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

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

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

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                               <  1:3764557/71‑1 (MQ=255)

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

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

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

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

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

                          gtgtGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGc                           >  1:1778176/1‑70 (MQ=11)

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

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

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

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                         <  1:3921856/71‑1 (MQ=255)

                               gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                      <  1:4303449/70‑1 (MQ=255)

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

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

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

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

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

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

                                cTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                    <  1:510530/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                      aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg              >  1:2637910/1‑71 (MQ=255)

                                      aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg               <  1:1422979/70‑1 (MQ=255)

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

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

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

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

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

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

                                        ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg             >  1:2805940/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)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca           >  1:1999065/1‑71 (MQ=255)

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

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

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

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

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

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc      >  1:1612976/1‑71 (MQ=255)

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

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

                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc       <  1:346886/68‑1 (MQ=255)

                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc       <  1:434683/68‑1 (MQ=255)

                                                  cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg  <  1:3976402/71‑1 (MQ=255)

                                                  cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg     >  1:4020848/1‑68 (MQ=255)

                                                  cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg     >  1:2242058/1‑68 (MQ=255)

                                                  cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg  >  1:2997063/1‑71 (MQ=255)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATG  >  Exported/1763502‑1763622

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