BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,636	A→G	28.9%	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	28.9%	58.0 / 34.7	45	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (16/16); new base G (3/10); total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTC  >  Exported/1763568‑1763702

ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAggg                                                                  <  1:2751618/71‑1 (MQ=255)

    ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                <  1:2588313/69‑1 (MQ=255)

         gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct                                                         <  1:1010439/71‑1 (MQ=255)

         gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct                                                         <  1:1006842/71‑1 (MQ=255)

         gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct                                                         <  1:1499787/71‑1 (MQ=255)

         gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctct                                                         <  1:752786/71‑1 (MQ=39)

          ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAggg                                                                  <  1:248207/61‑1 (MQ=255)

          ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCtt                                                        >  1:2297303/1‑71 (MQ=39)

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

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

             ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                      <  1:3009831/70‑1 (MQ=40)

              tACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                     <  1:2104634/70‑1 (MQ=255)

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

               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                      >  1:1926428/1‑68 (MQ=40)

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

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

                 cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                                <  1:2766984/56‑1 (MQ=255)

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                    >  1:1561265/1‑50 (MQ=255)

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt                                                  <  1:2827123/68‑1 (MQ=255)

                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                               >  1:316080/1‑71 (MQ=39)

                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct                                                         >  1:1149884/1‑59 (MQ=255)

                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt                                                        <  1:2824125/60‑1 (MQ=255)

                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTg                                                >  1:3035806/1‑68 (MQ=255)

                      caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                     <  1:236821/62‑1 (MQ=255)

                      caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                     <  1:1486736/62‑1 (MQ=255)

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

                                aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                   >  1:2187050/1‑70 (MQ=35)

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

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

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

                                      tCTGGGCACATCCGATGGCAAGAGGCCATAATGTCCGCCTCTTTGGTCTTGCGACATTAttc                                     <  1:1912344/62‑3 (MQ=11)

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

                                          ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                            <  1:1213337/67‑1 (MQ=35)

                                            cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGt                      <  1:176263/71‑1 (MQ=34)

                                            cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGt                      >  1:435996/1‑71 (MQ=34)

                                             acaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt                     >  1:770106/1‑71 (MQ=34)

                                              caTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                       <  1:2324991/52‑1 (MQ=35)

                                               aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                   <  1:1198541/71‑1 (MQ=34)

                                                tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                   >  1:1715817/1‑70 (MQ=35)

                                                  cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                >  1:1095097/1‑71 (MQ=34)

                                                   gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt               <  1:1334790/71‑1 (MQ=34)

                                                     tGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag             >  1:288036/1‑71 (MQ=34)

                                                         aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt         <  1:1231930/71‑1 (MQ=34)

                                                          agagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                        <  1:253570/55‑1 (MQ=35)

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

                                                                 cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc  >  1:1609059/1‑70 (MQ=35)

                                                                  gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                          >  1:2810188/1‑45 (MQ=38)

                                                                  gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                          >  1:1016604/1‑45 (MQ=38)

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTC  >  Exported/1763568‑1763702

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