BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,556	T→C	48.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	48.8%	1.9 / 262.1	168	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (46/40); new base C (41/41); total (87/81)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 169 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGG  >  Exported/1763502‑1763623

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

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

ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                      <  1:3898091/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:3639358/71‑1 (MQ=34)

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

    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   >  1:1961792/1‑69 (MQ=35)

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

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

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

       ggACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                              >  1:1883677/1‑71 (MQ=21)

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

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                               <  1:2686093/69‑1 (MQ=25)

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

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                            >  1:3155690/1‑70 (MQ=35)

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

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

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

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

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

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

               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                     <  1:3593602/56‑1 (MQ=33)

               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcccc                                      <  1:2069517/71‑1 (MQ=34)

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

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

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

                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                    >  1:43676/1‑71 (MQ=34)

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

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

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

                   aGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacc                                   <  1:2180373/70‑1 (MQ=11)

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

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

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                  <  1:4278501/71‑1 (MQ=25)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                             tgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                        <  1:74407/71‑1 (MQ=11)

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

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

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                        >  1:1558553/1‑71 (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:1849771/71‑2 (MQ=11)

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

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

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

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

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

                                    tCATCCTCTCAGACCAGCCAGGGATCGGCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                          <  1:896226/62‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                            tcAGACCAGCCAGGGATCGTCGCCTTGTTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt         >  1:2453252/1‑71 (MQ=255)

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

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

                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc        >  1:414979/1‑71 (MQ=255)

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

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

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

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc       >  1:3478404/1‑71 (MQ=255)

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc        >  1:620793/1‑70 (MQ=255)

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

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

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

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

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                        >  1:4411094/1‑38 (MQ=255)

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

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

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

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

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

                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTAccc                                       <  1:3280584/36‑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)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGG  >  Exported/1763502‑1763623

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