BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,594	A→T	6.5%	noncoding (250/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,594	0	A	T	6.5%	458.0 / 14.8	150	noncoding (250/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (51/89); new base T (5/5); total (56/94)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 158 total aligned reads displayed.

TGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAT  >  Exported/1763529‑1763662

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

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

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

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

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

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

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

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

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

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

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

    gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                                                             >  1:2479877/1‑71 (MQ=255)

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

      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)

         tCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                                                           <  1:386002/68‑1 (MQ=255)

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

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

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

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

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

              ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc                                                    <  1:4074086/70‑1 (MQ=255)

              ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc                                                    <  1:4464196/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                             gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAag                                     <  1:3984355/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)

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

                                 tcttcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGc                                 <  1:4286242/67‑1 (MQ=255)

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

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

                                  cgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                                                            <  1:4713566/42‑1 (MQ=255)

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

                                   gtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCg                              <  1:2155986/71‑1 (MQ=255)

                                   gtcgCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGATAGCGTGAGGTCCg                              <  1:3956989/71‑1 (MQ=255)

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

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

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

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

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

                                         ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                         <  1:596083/70‑1 (MQ=255)

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

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

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

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

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

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

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

                                             tGAGCCGTTTCCCCACCAACAAGCTAATCCCTTCTGGGCACATCCGATGGCAAGAGGCCCGa                             <  1:907629/62‑1 (MQ=255)

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

                                                    ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg             <  1:3529859/71‑1 (MQ=255)

                                                    ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg             >  1:3055949/1‑71 (MQ=255)

                                                      aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGcc                                >  1:166114/1‑50 (MQ=255)

                                                      aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc           >  1:469849/1‑71 (MQ=39)

                                                       ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCg                              <  1:3584796/51‑1 (MQ=255)

                                                       ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCg                              <  1:4460798/51‑1 (MQ=255)

                                                       ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt          <  1:3615093/71‑1 (MQ=255)

                                                        cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAaga                                    <  1:2336716/44‑1 (MQ=255)

                                                        cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAaga                                    <  1:3027657/44‑1 (MQ=255)

                                                         ccaccaACAAGCTAATCCCATCTGGGCACATCCGATGTCAAGAGGCCCg                              <  1:2519365/49‑1 (MQ=255)

                                                         ccaccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCg                              <  1:4026176/49‑1 (MQ=255)

                                                         ccaccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCg                              <  1:2763391/49‑1 (MQ=255)

                                                            ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGc                                 >  1:4520048/1‑43 (MQ=255)

                                                            ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGc                                 >  1:801493/1‑43 (MQ=255)

                                                            ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGc                                 >  1:1339575/1‑43 (MQ=255)

                                                            ccaACAAGCTAATCCCATCTGGGCACATCAGATGGCAAGAGGCCCGAAGGTCCCCctct                 <  1:3667579/59‑1 (MQ=255)

                                                             caacaaGCTAATCCCATCTGGGCACATCCGATGGCa                                       <  1:763185/36‑1 (MQ=255)

                                                               acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGccc                               >  1:824054/1‑42 (MQ=38)

                                                               acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCtt         <  1:2742813/64‑1 (MQ=35)

                                                                caaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt  >  1:4151963/1‑70 (MQ=37)

                                                                 aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCg                              >  1:1935655/1‑41 (MQ=32)

TGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAT  >  Exported/1763529‑1763662

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