BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,359,867	G→A	28.4%	noncoding (1228/2905 nt)	rrlA ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,359,867	0	G	A	28.4%	150.6 / 72.2	95	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (39/29); new base A (11/16); total (50/45)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCC  >  Exported/2359815‑2359933

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                  >  1:2859990/1‑71 (MQ=9)

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                  >  1:2718009/1‑71 (MQ=9)

 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                 <  1:908601/71‑1 (MQ=255)

  cGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                  >  1:1200510/1‑69 (MQ=255)

      ttgtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACATTCACAGGc                                             <  1:2692419/67‑1 (MQ=11)

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                            >  1:1231051/1‑70 (MQ=35)

        gtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCtt                                           >  1:3217660/2‑70 (MQ=11)

        gtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCtt                                           >  1:4022429/2‑70 (MQ=11)

        gtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCtt                                           >  1:2057340/2‑70 (MQ=11)

        gtCGTTACTCATGTCAGCATTCGCACTTCTGATACCTCCAGCAAGCTTCTCAACTCACCTTCACAGGCt                                            >  1:1578302/2‑69 (MQ=11)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                          >  1:2999006/1‑70 (MQ=35)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                          <  1:2065365/70‑1 (MQ=255)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                          <  1:1635439/70‑1 (MQ=255)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                          <  1:1470147/69‑1 (MQ=255)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                        <  1:1590780/70‑1 (MQ=35)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                        <  1:4146891/70‑1 (MQ=35)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                       <  1:131959/71‑1 (MQ=34)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                       <  1:4355287/71‑1 (MQ=255)

           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                       <  1:612577/71‑1 (MQ=255)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                      <  1:3441075/71‑1 (MQ=34)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                       <  1:1411050/70‑1 (MQ=255)

             tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                  >  1:3052806/1‑58 (MQ=31)

             tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                  >  1:2958641/1‑58 (MQ=31)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAAc                                    >  1:1037015/1‑71 (MQ=34)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     <  1:3019401/70‑1 (MQ=255)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                >  1:4001399/1‑71 (MQ=34)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                <  1:3961388/70‑1 (MQ=255)

                     tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                              >  1:1655544/1‑70 (MQ=35)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                         >  1:1654707/1‑57 (MQ=35)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            >  1:1609990/1‑70 (MQ=35)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                           >  1:1111928/1‑71 (MQ=34)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                           >  1:960316/1‑71 (MQ=34)

                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:1042726/70‑1 (MQ=255)

                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                           >  1:1680830/1‑71 (MQ=255)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                  <  1:2185006/47‑1 (MQ=32)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                  <  1:364539/47‑1 (MQ=32)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                  <  1:949872/47‑1 (MQ=32)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                  <  1:2038990/47‑1 (MQ=32)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                <  1:1020027/65‑1 (MQ=35)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                <  1:1106051/65‑1 (MQ=35)

                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                              >  1:886297/1‑65 (MQ=35)

                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                              >  1:3275944/1‑65 (MQ=35)

                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                              >  1:2956437/1‑65 (MQ=35)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                        >  1:1925689/1‑70 (MQ=35)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                        <  1:984305/54‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                        <  1:883135/54‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                        <  1:705658/54‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:1960939/66‑1 (MQ=255)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:3457662/66‑1 (MQ=255)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                      <  1:2806292/55‑1 (MQ=35)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                      <  1:1996104/55‑1 (MQ=35)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                      <  1:3916790/55‑1 (MQ=35)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:2815496/1‑71 (MQ=34)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:1133134/1‑71 (MQ=34)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:3927596/1‑71 (MQ=34)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                      <  1:1535631/54‑1 (MQ=35)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:4351365/1‑70 (MQ=35)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:4289060/1‑70 (MQ=255)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                          >  1:1700398/1‑65 (MQ=255)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                          >  1:406779/1‑65 (MQ=255)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                          >  1:3984578/1‑65 (MQ=255)

                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                          >  1:2186018/1‑65 (MQ=255)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                   <  1:2427504/71‑1 (MQ=34)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                   <  1:2292556/71‑1 (MQ=34)

                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                           >  1:2854243/1‑61 (MQ=35)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                  <  1:3665446/69‑1 (MQ=35)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                  <  1:578017/69‑1 (MQ=35)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                  <  1:156259/69‑1 (MQ=35)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                  <  1:2249574/69‑1 (MQ=35)

                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca               >  1:406991/1‑71 (MQ=39)

                                   tCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                 >  1:1493806/1‑37 (MQ=11)

                                   tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                >  1:3198516/1‑68 (MQ=11)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                  >  1:3510110/1‑35 (MQ=33)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                        >  1:700750/1‑45 (MQ=38)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                        >  1:2754388/1‑45 (MQ=38)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                        >  1:348022/1‑45 (MQ=38)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                 <  1:3994897/68‑1 (MQ=37)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                <  1:1753502/69‑1 (MQ=39)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca               >  1:2736779/1‑70 (MQ=40)

                                    cTGATACCTCCACCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                 <  1:2057722/68‑1 (MQ=255)

                                      gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg            >  1:3777192/1‑71 (MQ=40)

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt           >  1:210117/1‑71 (MQ=255)

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg            <  1:2685644/70‑1 (MQ=255)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                          <  1:2498028/51‑1 (MQ=11)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACTCTccc                              >  1:3006355/1‑45 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca               >  1:3673547/1‑60 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca               >  1:1618810/1‑60 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca               >  1:699718/1‑60 (MQ=255)

                                               aGCATGCCTCACAGCACACCTTCACATGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc   >  1:2651377/1‑71 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                   <  1:552288/38‑1 (MQ=40)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                   <  1:2607487/38‑1 (MQ=40)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                   <  1:1648801/54‑1 (MQ=35)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                   <  1:4282818/54‑1 (MQ=35)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc  >  1:1661703/1‑71 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                 >  1:4442921/1‑54 (MQ=38)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCC  >  Exported/2359815‑2359933

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