BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,556	T→C	38.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
*	minE	1,763,556	0	T	C	38.8%	112.2 / 148.4	134	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (31/51); new base C (21/31); total (52/82)
Fisher's exact test for biased strand distribution p-value = 8.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 136 total aligned reads displayed.

AGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  minE/1763501‑1763624

aGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       >  1:231395/1‑71 (MQ=31)

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

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

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

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

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

  gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                     <  1:3875200/71‑1 (MQ=34)

  gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      <  1:3475129/70‑1 (MQ=18)

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

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

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

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

       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                    >  1:3116312/1‑67 (MQ=25)

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

         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                                  >  1:3156770/1‑67 (MQ=35)

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

         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                               <  1:3711867/70‑1 (MQ=35)

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

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

          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:2256116/71‑1 (MQ=21)

           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      >  1:1375748/1‑61 (MQ=33)

           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:968054/71‑1 (MQ=34)

           ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:1879222/71‑1 (MQ=34)

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

            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                       >  1:275890/1‑59 (MQ=9)

            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:37970/70‑1 (MQ=25)

            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTAAGACCAGCCAGGGATCGTCGCCtt                                                       >  1:4458129/1‑59 (MQ=11)

            cGTGTCTCAGTTCCAGTGTGGCGGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:2337357/70‑1 (MQ=11)

             gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:1179135/69‑1 (MQ=35)

              tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       >  1:104115/1‑57 (MQ=31)

               gtCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcct                                        >  1:2340067/1‑70 (MQ=11)

               gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:2270110/66‑1 (MQ=35)

                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGTGATCGTCGCCTTGGTGAGCCGTTAcccc                                       <  1:1796926/71‑1 (MQ=255)

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

                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         <  1:1340872/69‑1 (MQ=25)

                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAcccc                                       >  1:1834110/1‑69 (MQ=25)

                   cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      <  1:121742/53‑1 (MQ=16)

                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:4026514/71‑1 (MQ=35)

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

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

                       tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:1248684/1‑70 (MQ=255)

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

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

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

                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:442141/1‑69 (MQ=255)

                        ccAGTGTGACTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                <  1:3710619/70‑1 (MQ=255)

                           gtgtGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:2522513/71‑1 (MQ=255)

                           gtgtGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:1423163/71‑1 (MQ=255)

                           gtgtGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCCCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:1415899/71‑1 (MQ=255)

                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:949729/46‑1 (MQ=38)

                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGAGCGTCGCCtt                                                       <  1:1645852/43‑1 (MQ=12)

                            tgtgGCTGGTCATCCTCTCAGACAAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:3506360/70‑1 (MQ=255)

                             gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           >  1:3625238/1‑70 (MQ=255)

                              tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:2992166/41‑1 (MQ=32)

                              tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:3022897/68‑1 (MQ=255)

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

                                gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGTCGTTACCCCACCAACAAGCTAATcc                       >  1:392958/1‑71 (MQ=14)

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

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

                                   ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:1227319/63‑1 (MQ=255)

                                    gtcgtcCTCTCAGACCAGCTACTGATCGTCGCCTTGGt                                                    >  1:3235503/1‑38 (MQ=14)

                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                      <  1:88788/68‑1 (MQ=255)

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

                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                   <  1:2151579/71‑1 (MQ=255)

                                     tCATCCTCTCAGACGAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                    <  1:31207/69‑1 (MQ=255)

                                     tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:4252038/61‑1 (MQ=255)

                                      cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                 >  1:15039/1‑71 (MQ=255)

                                       aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:3392074/1‑44 (MQ=38)

                                       aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                 <  1:4431201/70‑1 (MQ=255)

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

                                        tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                           <  1:1957442/43‑1 (MQ=38)

                                        tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:2463882/1‑71 (MQ=255)

                                          ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         >  1:1499679/1‑59 (MQ=255)

                                          ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca             >  1:466717/1‑71 (MQ=255)

                                           tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                <  1:822446/35‑1 (MQ=40)

                                           tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                <  1:2725331/35‑1 (MQ=40)

                                           tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca             >  1:4356638/1‑70 (MQ=255)

                                             tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                      <  1:2395770/59‑1 (MQ=255)

                                             tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcac            <  1:1003515/69‑1 (MQ=255)

                                             tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGaaca           >  1:1637830/1‑70 (MQ=255)

                                             tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:555795/36‑1 (MQ=255)

                                             tcAGACCAGCCAGGGAGCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:4034789/53‑1 (MQ=17)

                                             tcAGAACAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACACCACCAACAAGCTAATCCCATCTGGGCACAt          >  1:1528597/1‑71 (MQ=255)

                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          >  1:2954028/1‑53 (MQ=255)

                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:2463025/71‑1 (MQ=255)

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

                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGACACATc         <  1:5788/70‑1 (MQ=255)

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

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

                                                gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg       <  1:3573253/71‑1 (MQ=255)

                                                gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg       <  1:444013/71‑1 (MQ=255)

                                                gACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                <  1:3406510/46‑1 (MQ=255)

                                                 aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa      <  1:3886377/71‑1 (MQ=255)

                                                 aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                             <  1:747355/48‑1 (MQ=255)

                                                 aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCAATCTGGGCACATCCg       >  1:2750356/1‑70 (MQ=255)

                                                  ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa      <  1:3905719/70‑1 (MQ=255)

                                                   cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt     <  1:1697908/70‑1 (MQ=255)

                                                   cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           >  1:3759173/1‑48 (MQ=255)

                                                    aGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          <  1:4245106/48‑1 (MQ=255)

                                                       tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  >  1:1783441/1‑69 (MQ=255)

AGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  minE/1763501‑1763624

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