BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,763,556	T→C	44.3%	noncoding (288/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,556	0	T	C	44.3%	37.4 / 182.9	131	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (43/30); new base C (31/27); total (74/57)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 131 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  Exported/1763502‑1763624

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

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

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

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

  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:5579719/71‑1 (MQ=34)

   gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      >  1:1851735/1‑68 (MQ=18)

    tCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGtag                                                  <  1:5872820/71‑3 (MQ=11)

    tCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGtag                                                  <  1:3191511/71‑3 (MQ=11)

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

       ggACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                               >  1:3133810/1‑71 (MQ=34)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       >  1:360905/1‑62 (MQ=31)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       >  1:1824045/1‑62 (MQ=31)

        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                               <  1:2609494/70‑1 (MQ=25)

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

         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:6318659/1‑71 (MQ=21)

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGGGCCttt                                            <  1:3716574/71‑1 (MQ=11)

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCttt                                            <  1:6848372/71‑1 (MQ=11)

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCttt                                            <  1:2666860/71‑1 (MQ=11)

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:3558627/1‑70 (MQ=35)

          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:4818023/1‑70 (MQ=35)

           cGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCttt                                            <  1:6856590/70‑1 (MQ=11)

           cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                           <  1:5540168/71‑1 (MQ=34)

           cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                           <  1:2758587/71‑1 (MQ=34)

             tgtCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTAcc                                         <  1:1863689/71‑1 (MQ=11)

               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcccc                                       <  1:2834785/71‑1 (MQ=34)

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

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    >  1:2634188/1‑70 (MQ=25)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    >  1:6501513/1‑70 (MQ=25)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    >  1:6682766/1‑70 (MQ=25)

                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    >  1:3623070/1‑70 (MQ=25)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      <  1:1234643/50‑1 (MQ=255)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      <  1:4605383/50‑1 (MQ=255)

                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      <  1:2184870/50‑1 (MQ=255)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                <  1:5079652/70‑1 (MQ=255)

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

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

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

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                    >  1:1002540/1‑50 (MQ=255)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                    >  1:6685291/1‑50 (MQ=255)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                    >  1:6297612/1‑50 (MQ=255)

                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                    >  1:1010523/1‑50 (MQ=255)

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

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

                          gtgtGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:5351272/71‑1 (MQ=255)

                           tgtgGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCATTTACCtcacca                                   <  1:4904356/63‑1 (MQ=11)

                           tgtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTTGTGAGcc                                               <  1:3728714/51‑1 (MQ=11)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     <  1:3078123/43‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     <  1:1204011/43‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     <  1:1959860/43‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     <  1:6818050/43‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     <  1:3898897/43‑1 (MQ=255)

                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     <  1:3916158/43‑1 (MQ=255)

                               tcTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      >  1:5732985/2‑40 (MQ=37)

                               gCTTGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      >  1:3705274/1‑40 (MQ=255)

                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      >  1:5830579/1‑40 (MQ=37)

                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      >  1:4559303/1‑40 (MQ=37)

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

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

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGa                                                  >  1:1832984/1‑42 (MQ=255)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    <  1:3239841/56‑1 (MQ=255)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    <  1:4512604/56‑1 (MQ=255)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    <  1:2417011/56‑1 (MQ=255)

                                 tGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    <  1:5815870/56‑1 (MQ=255)

                                  ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                 <  1:6696546/42‑1 (MQ=255)

                                  ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                 <  1:5165486/42‑1 (MQ=255)

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

                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                >  1:6562542/1‑71 (MQ=255)

                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                >  1:4406056/1‑71 (MQ=255)

                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                 >  1:4600794/1‑70 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:4182043/1‑41 (MQ=38)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:6387257/1‑71 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:324328/1‑71 (MQ=255)

                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:5985915/1‑71 (MQ=255)

                                       tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               <  1:6563228/71‑1 (MQ=255)

                                       tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               <  1:3887904/71‑1 (MQ=255)

                                        ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTGCCCCACCAACAAGCTAATCCCATCTggg               >  1:436002/1‑70 (MQ=14)

                                        ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:3798902/1‑70 (MQ=255)

                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc              <  1:1995088/70‑1 (MQ=255)

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

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:905746/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:925043/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:4780432/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:4623624/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:4296124/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:4159876/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:3859846/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:2222523/1‑38 (MQ=255)

                                            tcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:5363663/1‑38 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                      >  1:1572595/1‑42 (MQ=38)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc              >  1:3075915/1‑66 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc              >  1:1005194/1‑66 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc              >  1:3045631/1‑66 (MQ=255)

                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc              >  1:3837632/1‑66 (MQ=255)

                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc         >  1:6828435/1‑71 (MQ=255)

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

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

                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGAACATcc        <  1:899174/71‑1 (MQ=255)

                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc         >  1:1037817/1‑70 (MQ=255)

                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa      >  1:1523031/1‑71 (MQ=255)

                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  <  1:4766656/71‑1 (MQ=255)

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  Exported/1763502‑1763624

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