BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,636	A→G	10.0%	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	10.0%	132.6 / 6.8	50	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (23/22); new base G (3/2); total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCC  >  Exported/1763568‑1763699

ccTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAggg                                                               >  1:365259/1‑71 (MQ=255)

 cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                              <  1:70406/71‑1 (MQ=255)

 cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                              <  1:1499909/71‑1 (MQ=255)

   tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc                                                            >  1:619897/1‑71 (MQ=40)

      tGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccgcc                                                         <  1:3428773/71‑1 (MQ=255)

      tGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccccc                                                         <  1:1808259/71‑1 (MQ=39)

       gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                        >  1:562958/1‑71 (MQ=255)

       gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                        >  1:3553934/1‑71 (MQ=255)

       gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                        >  1:3059507/1‑71 (MQ=39)

       gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                        >  1:1043056/1‑71 (MQ=39)

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

            gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg                                                   <  1:1981182/71‑1 (MQ=255)

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

              tACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                  >  1:3519608/1‑70 (MQ=40)

              tACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                 >  1:1472931/1‑71 (MQ=39)

              tACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTAGGt                                                 >  1:2398251/1‑71 (MQ=255)

              tACCCCACCAACAAGCTAATCACATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc                                                            >  1:291792/1‑60 (MQ=255)

                ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTc                                                <  1:2903739/70‑1 (MQ=40)

                 cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCtt                                                     <  1:2280313/64‑1 (MQ=40)

                 cccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCtt                                              <  1:1445576/71‑1 (MQ=39)

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

                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg                                           <  1:236455/71‑1 (MQ=255)

                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                             >  1:231286/1‑53 (MQ=255)

                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                             >  1:162511/1‑53 (MQ=255)

                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTg                                             >  1:3682720/1‑69 (MQ=40)

                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGa                                          >  1:264813/1‑71 (MQ=255)

                          aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                     <  1:1017426/71‑1 (MQ=35)

                           aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:1708801/71‑1 (MQ=34)

                                   ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                             >  1:1744232/1‑70 (MQ=35)

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

                                     aTCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                               <  1:2351618/66‑1 (MQ=35)

                                        tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                        >  1:2342178/1‑54 (MQ=31)

                                          ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                           >  1:1416810/1‑65 (MQ=35)

                                           gCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAcc                     <  1:2017838/70‑1 (MQ=35)

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

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

                                             acaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCCTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt                  <  1:1038094/71‑1 (MQ=21)

                                              caTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                 <  1:464715/71‑1 (MQ=34)

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

                                                   gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                             <  1:683515/54‑1 (MQ=35)

                                                    aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagta           <  1:1134827/71‑1 (MQ=34)

                                                       gCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt         <  1:2489207/70‑1 (MQ=35)

                                                         aaGCGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt      >  1:1332589/1‑71 (MQ=21)

                                                           gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTa       <  1:1372305/68‑1 (MQ=35)

                                                             ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                 <  1:3045166/56‑1 (MQ=35)

                                                             ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                 <  1:534421/56‑1 (MQ=35)

                                                               cccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTcc               >  1:2401362/1‑56 (MQ=35)

                                                               cccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTcc               >  1:406220/1‑56 (MQ=35)

                                                                 cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcc    >  1:2548283/1‑65 (MQ=35)

                                                                  gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc  >  1:3176647/1‑66 (MQ=35)

                                                                  gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc  >  1:2815220/1‑66 (MQ=35)

                                                                   aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              <  1:2032476/37‑1 (MQ=40)

                                                                   aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              <  1:199729/37‑1 (MQ=40)

CCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCC  >  Exported/1763568‑1763699

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