BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,580	G→T	11.2%	noncoding (264/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,580	0	G	T	11.2%	562.5 / 31.2	177	noncoding (264/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (91/66); new base T (8/12); total (99/78)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

ACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTC  >  Exported/1763511‑1763646

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

    tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                                               <  1:241780/71‑1 (MQ=34)

        tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                                           >  1:1452450/1‑71 (MQ=21)

         cAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcac                                                           >  1:1022699/1‑70 (MQ=11)

         cAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcac                                                           >  1:4572446/1‑70 (MQ=11)

         cAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcac                                                           >  1:4985712/1‑70 (MQ=11)

                    tgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                                               >  1:4628818/1‑71 (MQ=11)

                    tgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                                               >  1:1421740/1‑71 (MQ=11)

                    tgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                                               >  1:4048186/1‑71 (MQ=11)

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

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

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:4992091/51‑1 (MQ=11)

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:31845/51‑1 (MQ=11)

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:1834253/51‑1 (MQ=11)

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:212854/51‑1 (MQ=11)

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:1762217/51‑1 (MQ=11)

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:3681734/51‑1 (MQ=11)

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:2505864/51‑1 (MQ=11)

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:4646036/51‑1 (MQ=11)

                      gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                                                 <  1:1336918/51‑1 (MQ=11)

                       cTGGTCATCCTCTCAGCCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                                             >  1:4443353/1‑70 (MQ=255)

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

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

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

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

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

                              tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                                         <  1:4379955/67‑1 (MQ=255)

                               ccTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                                         <  1:3524311/66‑1 (MQ=255)

                                ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca                                   >  1:1537970/1‑71 (MQ=255)

                                 tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                                  <  1:1320844/55‑1 (MQ=255)

                                 tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                                  <  1:1946977/55‑1 (MQ=255)

                                  ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                          >  1:3195243/1‑62 (MQ=255)

                                  ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                          <  1:2132786/62‑1 (MQ=255)

                                  ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                          <  1:3149788/62‑1 (MQ=255)

                                  ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                          <  1:5102724/62‑1 (MQ=255)

                                  ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                          <  1:2857380/62‑1 (MQ=255)

                                  ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                          <  1:3928392/62‑1 (MQ=255)

                                  ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                          <  1:1723317/62‑1 (MQ=255)

                                  ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcac                                  >  1:4233531/1‑70 (MQ=255)

                                   tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                                        >  1:2044302/1‑47 (MQ=255)

                                   tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                                        >  1:133872/1‑47 (MQ=255)

                                   tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACaaa                                                        >  1:4117983/1‑47 (MQ=255)

                                    cAGACCATCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc                               >  1:4784450/1‑71 (MQ=255)

                                    cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc                               >  1:4398980/1‑71 (MQ=255)

                                    cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc                               >  1:1234723/1‑71 (MQ=255)

                                    cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc                               >  1:4464581/1‑71 (MQ=255)

                                    cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                                        >  1:1866600/1‑62 (MQ=255)

                                    cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                                        >  1:612096/1‑62 (MQ=255)

                                    cAGACCAGCCAGGGATCGTCGCCTTGGTGACCCGTTACCCCACCAACAAGCAAATCCCATCt                                        >  1:4051324/1‑62 (MQ=17)

                                     aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca                                 <  1:2642919/68‑1 (MQ=255)

                                     aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca                                 <  1:4346734/68‑1 (MQ=255)

                                     aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca                                 <  1:4866292/68‑1 (MQ=255)

                                     aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca                                 <  1:2888818/68‑1 (MQ=255)

                                     aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:4733271/1‑71 (MQ=255)

                                     aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:3721379/1‑71 (MQ=255)

                                           gccaGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg                         >  1:3858277/4‑70 (MQ=255)

                                           gCTGGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              <  1:3432643/65‑1 (MQ=255)

                                           gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              <  1:4416314/65‑1 (MQ=255)

                                           gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              <  1:1960879/65‑1 (MQ=255)

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

                                                ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAagag                   >  1:3070544/1‑71 (MQ=255)

                                                ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAagag                   >  1:5191429/1‑71 (MQ=255)

                                                 gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa                            >  1:2264995/1‑61 (MQ=255)

                                                 gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCa                       >  1:210461/1‑66 (MQ=255)

                                                  aTCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg                          >  1:2907781/1‑62 (MQ=255)

                                                  aTCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg                          >  1:4323347/1‑62 (MQ=255)

                                                  aTCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg                          >  1:3287829/1‑62 (MQ=255)

                                                   tcgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc                        >  1:2727332/1‑63 (MQ=255)

                                                   tcgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc                        >  1:2027733/1‑63 (MQ=255)

                                                   tcgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc                        >  1:4956541/1‑63 (MQ=255)

                                                   tcgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc                        >  1:339224/1‑63 (MQ=255)

                                                     gtcgCCTTGGTGCGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc                               >  1:2231841/1‑54 (MQ=255)

                                                     gtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                                            <  1:4999882/41‑1 (MQ=255)

                                                     gtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg                          >  1:3014534/1‑59 (MQ=255)

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

                                                       cgCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATCTGGGCACATc                               >  1:4173860/1‑52 (MQ=255)

                                                       cgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                                         <  1:2266510/42‑1 (MQ=255)

                                                       cgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa                            >  1:153873/1‑55 (MQ=255)

                                                       cgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa                            >  1:2203022/1‑55 (MQ=255)

                                                       cgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa                            >  1:720136/1‑55 (MQ=255)

                                                       cgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa                            >  1:3663828/1‑55 (MQ=255)

                                                        gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                          <  1:3083725/40‑1 (MQ=255)

                                                        gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGagg           >  1:4551279/1‑69 (MQ=255)

                                                        gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGagg           >  1:1370720/1‑69 (MQ=255)

                                                        gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGagg           >  1:2352392/1‑69 (MQ=255)

                                                          cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCTAGGGt         >  1:1028553/1‑71 (MQ=255)

                                                           ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc        >  1:1251940/1‑71 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:4491526/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:3748794/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:4827231/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:3537462/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:3041228/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:1797883/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:5184593/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:1731288/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:1139682/1‑46 (MQ=255)

                                                              gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                              >  1:847606/1‑46 (MQ=255)

                                                                 aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCa                       <  1:2764828/50‑1 (MQ=255)

                                                                 aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc  <  1:2904722/71‑1 (MQ=39)

                                                                 aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc  <  1:2615005/71‑1 (MQ=39)

ACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTC  >  Exported/1763511‑1763646

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