BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,636	A→G	17.3%	noncoding (208/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,636	0	A	G	17.3%	227.0 / 41.7	98	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (36/45); new base G (12/5); total (48/50)
Fisher's exact test for biased strand distribution p-value = 6.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 104 total aligned reads displayed.

GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAT  >  Exported/1763567‑1763705

gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                       <  1:5597178/70‑1 (MQ=255)

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

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

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

  gttggtggggtaacggctctccaaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                                    >  1:4753940‑M2/25‑71 (MQ=255)

   ttggtggggtaacggctcaccaaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                   >  1:4061089‑M2/24‑71 (MQ=255)

   ttggtggggtaacggctcaccaaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                   >  1:6139877‑M2/24‑71 (MQ=255)

   TTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                      <  1:4056081‑M1/68‑1 (MQ=255)

   TTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                      <  1:6420567‑M1/68‑1 (MQ=255)

   TTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                      <  1:5211833‑M1/68‑1 (MQ=255)

   TTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                      <  1:1269480‑M1/68‑1 (MQ=255)

   TTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                      <  1:230057‑M1/68‑1 (MQ=255)

     GGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                                 <  1:5186271‑M1/71‑1 (MQ=255)

      GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                                 >  1:6786038‑M1/1‑70 (MQ=255)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                       >  1:2208239/1‑61 (MQ=40)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc                                                             >  1:4959065/1‑71 (MQ=255)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc                                                             >  1:5119212/1‑71 (MQ=255)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc                                                             >  1:6550657/1‑71 (MQ=255)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                              <  1:400937/70‑1 (MQ=255)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                                 >  1:2451989/1‑67 (MQ=40)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCGCct                                                              <  1:4594601/70‑1 (MQ=255)

         aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                              <  1:3468332/70‑1 (MQ=40)

           ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt                                                           <  1:5204161/71‑1 (MQ=255)

             gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg                                                         >  1:5002111/1‑71 (MQ=255)

             gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg                                                         >  1:4572443/1‑71 (MQ=255)

             gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg                                                         >  1:2862609/1‑71 (MQ=255)

             gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                         <  1:3646992/71‑1 (MQ=39)

             gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                         <  1:864607/71‑1 (MQ=39)

              ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                        >  1:5841430/1‑71 (MQ=39)

                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGt                                                       <  1:642248/70‑1 (MQ=255)

                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                      >  1:3004064/1‑71 (MQ=255)

                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                      >  1:1976476/1‑71 (MQ=255)

                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                       <  1:5514979/70‑1 (MQ=40)

                aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                       <  1:3522141/70‑1 (MQ=40)

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

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

                 ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                     <  1:5144198/71‑1 (MQ=39)

                  cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                     <  1:4949350/70‑1 (MQ=40)

                    caccaaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc                                                      >  1:3474528‑M2/7‑67 (MQ=255)

                     accaaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                 <  1:4844489‑M2/66‑1 (MQ=255)

                      ccaaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                       >  1:2257199‑M2/5‑48 (MQ=255)

                      ccaaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                       >  1:6426969‑M2/5‑48 (MQ=255)

                       caaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAc                                               >  1:4582595‑M2/4‑71 (MQ=255)

                       caaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCACCTCTTTGGTCTTGc                                                  >  1:1956287‑M2/4‑68 (MQ=255)

                        aaCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                     <  1:364118‑M2/62‑1 (MQ=255)

                         aCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGa                                                <  1:6848400‑M2/67‑1 (MQ=255)

                         aCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGa                                                <  1:2058613‑M2/67‑1 (MQ=255)

                         aCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:3600293‑M2/70‑1 (MQ=255)

                         aCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:5827897‑M2/70‑1 (MQ=255)

                         aCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             >  1:2967874‑M2/2‑71 (MQ=255)

                         aCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:3381385‑M2/70‑1 (MQ=255)

                         aCAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:3761682‑M2/70‑1 (MQ=255)

                          caaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAtt                                            <  1:2248209/71‑1 (MQ=38)

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

                             gatAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                         <  1:3337342/69‑1 (MQ=21)

                             gCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                         <  1:797245/71‑1 (MQ=34)

                               tAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCg                                       >  1:151066/1‑71 (MQ=34)

                                 aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              >  1:2003736/1‑62 (MQ=31)

                                 aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              >  1:900833/1‑62 (MQ=31)

                                 aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              >  1:5285351/1‑62 (MQ=31)

                                 aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              >  1:5279141/1‑62 (MQ=31)

                                 aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:6588579/70‑1 (MQ=35)

                                    ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAtt                                  <  1:5293048/71‑1 (MQ=34)

                                    ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                   <  1:5288614/70‑1 (MQ=35)

                                       tCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                                     <  1:1830012/65‑1 (MQ=35)

                                        cTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                               >  1:3339360/1‑70 (MQ=35)

                                         tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                        >  1:2117387/1‑60 (MQ=35)

                                         tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                        >  1:3654497/1‑60 (MQ=35)

                                         tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                        >  1:2591598/1‑60 (MQ=35)

                                         tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                             >  1:4359207/1‑71 (MQ=34)

                                              acaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt                        <  1:1992679/71‑1 (MQ=34)

                                                aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:4144498/48‑1 (MQ=35)

                                                aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:6674689/48‑1 (MQ=35)

                                                 tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:6610230/47‑1 (MQ=35)

                                                 tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:4353704/47‑1 (MQ=35)

                                                 tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:6291518/47‑1 (MQ=35)

                                                 tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:544000/47‑1 (MQ=35)

                                                 tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                           >  1:803421/1‑49 (MQ=38)

                                                  ccGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                             <  1:3305356/46‑1 (MQ=35)

                                                   cGATGGCAAGAGGCCCGAAGGTCCACCTCTTTGGTCTTGCGACATTATGCGGTATTa                                 <  1:181476/57‑1 (MQ=255)

                                                    gATGGCAAGAGGCCCGAATGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                                 >  1:2237041/1‑56 (MQ=255)

                                                    gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                                 >  1:1751539/1‑56 (MQ=35)

                                                       ggCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              >  1:5028103/1‑40 (MQ=33)

                                                       ggCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                              >  1:3157910/1‑40 (MQ=33)

                                                        gCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:1744428/47‑1 (MQ=38)

                                                        gCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:2252863/47‑1 (MQ=38)

                                                         cAAGACGCCCGAAGGGCCCCCTCTTTGGTCTTGCGACAt                                             <  1:2183756/39‑1 (MQ=11)

                                                          aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                        >  1:4058435/1‑43 (MQ=38)

                                                           agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                       >  1:491836/1‑59 (MQ=35)

                                                               gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccccc       >  1:1318083/1‑71 (MQ=34)

                                                               gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccccc       >  1:239395/1‑71 (MQ=34)

                                                               gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc        >  1:5612629/1‑70 (MQ=35)

                                                                 ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:6373311/38‑1 (MQ=40)

                                                                 ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:1896628/38‑1 (MQ=40)

                                                                 ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                      <  1:5643427/38‑1 (MQ=40)

                                                                 ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc     <  1:3232421/71‑1 (MQ=34)

                                                                 ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc     >  1:4405309/1‑71 (MQ=34)

                                                                 ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc     >  1:606269/1‑71 (MQ=34)

                                                                    aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt  <  1:1911191/71‑1 (MQ=34)

                                                                    aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt  >  1:5834363/1‑71 (MQ=34)

GCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAT  >  Exported/1763567‑1763705

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