BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,580	G→T	9.1%	noncoding (264/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,580	0	G	T	9.1%	610.8 / 26.6	185	noncoding (264/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (76/92); new base T (10/7); total (86/99)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 188 total aligned reads displayed.

ACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTG  >  minE/1763511‑1763650

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

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

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

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

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

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

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

        tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                                                    >  1:2042295/1‑66 (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)

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

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

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

              ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                                                       >  1:3342261/1‑57 (MQ=255)

              ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                                                       >  1:3198613/1‑57 (MQ=255)

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

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

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

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

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

                    tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                                                   <  1:2810103/71‑1 (MQ=255)

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

                    tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                                                   >  1:1558553/1‑71 (MQ=255)

                    tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                                                   >  1:3116563/1‑71 (MQ=255)

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

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

                     ggCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                                               <  1:222373/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:2479877/1‑71 (MQ=255)

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

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

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

                          gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                                               <  1:473589/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                     aGACCTGCCAGGGATCGTCGCCTTGGTGAGCCGTTTCCCCACCAACAAGCTAATCCCATCTGGGCACATcc                                  <  1:3782029/71‑1 (MQ=14)

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

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

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

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

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

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

                                     aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                                               >  1:253790/1‑58 (MQ=255)

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

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

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

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

                                        ccagccagGGATCGTCGCCTTGGTGAGCCGTTAccc                                                                  <  1:3280584/36‑1 (MQ=255)

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

                                           gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg                             >  1:97095/1‑70 (MQ=255)

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

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

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

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

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

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

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

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

                                                 gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAaga                        >  1:1205961/1‑69 (MQ=255)

                                                 gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAaga                        >  1:2855736/1‑69 (MQ=255)

                                                  aTCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGc                     <  1:2377411/71‑1 (MQ=255)

                                                   tcgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGc                     <  1:1807675/70‑1 (MQ=255)

                                                    cgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGccc                   >  1:3673151/1‑71 (MQ=255)

                                                    cgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGccc                   >  1:4581910/1‑71 (MQ=255)

                                                      tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                                  <  1:3542142/54‑1 (MQ=255)

                                                      tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                                  <  1:3456009/54‑1 (MQ=255)

                                                      tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa                          <  1:514597/62‑1 (MQ=255)

                                                      tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa                          <  1:168675/62‑1 (MQ=255)

                                                      tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGa                 <  1:2389661/71‑1 (MQ=255)

                                                      tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGa                 <  1:3078324/71‑1 (MQ=255)

                                                      tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGa                 <  1:1788043/71‑1 (MQ=255)

                                                        gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGc                     <  1:297968/65‑1 (MQ=255)

                                                        gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg               <  1:1053618/71‑1 (MQ=255)

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

                                                            tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAag                         <  1:4700308/57‑1 (MQ=255)

                                                            tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAag                         <  1:3952314/57‑1 (MQ=255)

                                                            tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAgg                      >  1:4138414/1‑60 (MQ=255)

                                                            tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAgg                      >  1:711066/1‑60 (MQ=255)

                                                            tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAgg                      >  1:3172866/1‑60 (MQ=255)

                                                            tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc           >  1:4708828/1‑71 (MQ=255)

                                                             ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc          >  1:630173/1‑71 (MQ=39)

                                                             ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc          >  1:2290360/1‑71 (MQ=39)

                                                             ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc          >  1:245894/1‑71 (MQ=39)

                                                                 aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct       >  1:493510/1‑70 (MQ=40)

                                                                    cGTTACCCCACCAACAAGCTAATCTCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt    <  1:264689/70‑1 (MQ=255)

                                                                     gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg  >  1:1945801/1‑71 (MQ=39)

ACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTG  >  minE/1763511‑1763650

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