BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	34.5%	noncoding (1219/2905 nt)	rrlA ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,359,876	0	C	T	34.5%	232.3 / 139.5	139	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (41/50); new base T (32/16); total (73/66)
Fisher's exact test for biased strand distribution p-value = 2.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 139 total aligned reads displayed.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359815‑2359945

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:982880/1‑71 (MQ=31)

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              <  1:2382342/71‑1 (MQ=31)

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              <  1:212637/71‑1 (MQ=31)

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:1733050/1‑71 (MQ=31)

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:1395476/1‑71 (MQ=31)

   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                           <  1:1999949/71‑1 (MQ=34)

   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTACACAg                                                           <  1:3194566/71‑1 (MQ=21)

      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        >  1:1158943/1‑71 (MQ=34)

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

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACTTCACAGTACACCTTCACAGGCt                                                        >  1:3237563/1‑70 (MQ=11)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGttt                                                       <  1:1648548/70‑1 (MQ=11)

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

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

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

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:1474561/70‑1 (MQ=35)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:3828601/70‑1 (MQ=35)

             tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:3583028/1‑71 (MQ=255)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                             >  1:3135225/1‑54 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:3166627/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:2779957/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:3739743/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:3340452/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:551030/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:580225/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:1068439/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                            >  1:925775/1‑71 (MQ=255)

                  aTGTCAGCATTCGAACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                              <  1:1031413/69‑1 (MQ=25)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            <  1:3102318/70‑1 (MQ=35)

                    gTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                          >  1:1723789/1‑71 (MQ=34)

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    >  1:249889/1‑59 (MQ=35)

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    >  1:3513471/1‑59 (MQ=35)

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                         >  1:2916741/1‑70 (MQ=35)

                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCACTa                                       >  1:822378/1‑71 (MQ=11)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      <  1:2800807/71‑1 (MQ=34)

                        gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:1601524/1‑71 (MQ=255)

                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:795348/1‑71 (MQ=34)

                          aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:1478719/1‑70 (MQ=255)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      <  1:2955148/52‑1 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:1764836/1‑51 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:1344490/1‑51 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                                    >  1:1659931/1‑52 (MQ=11)

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

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACAcaa                                  >  1:1305419/1‑70 (MQ=11)

                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 >  1:827773/1‑70 (MQ=35)

                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                              <  1:2952889/71‑1 (MQ=34)

                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                             <  1:1142725/71‑1 (MQ=35)

                                 cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                             >  1:1334627/1‑70 (MQ=255)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                        <  1:3809034/43‑1 (MQ=14)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:3484645/48‑1 (MQ=11)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:1826415/48‑1 (MQ=11)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:551332/48‑1 (MQ=11)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:673706/48‑1 (MQ=11)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      >  1:2673193/1‑43 (MQ=38)

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

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

                                    cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                            <  1:1372797/69‑3 (MQ=11)

                                      gATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:1772659/1‑38 (MQ=14)

                                      gATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:1010879/1‑46 (MQ=11)

                                      gATACCTCCAGCATACCTCACAGTACAACTTCACAGGCTTACAGaa                                                 >  1:3446275/1‑46 (MQ=11)

                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                     >  1:3725823/1‑71 (MQ=255)

                                         aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                          >  1:3076725/1‑50 (MQ=11)

                                         aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                          >  1:1091262/1‑50 (MQ=11)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                    <  1:313531/71‑1 (MQ=255)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                    <  1:1110712/71‑1 (MQ=255)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                    <  1:1189670/71‑1 (MQ=255)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                    <  1:1839285/71‑1 (MQ=255)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                    <  1:1920107/71‑1 (MQ=255)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               <  1:2894417/58‑1 (MQ=35)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               <  1:2375323/58‑1 (MQ=35)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                         <  1:996272/64‑1 (MQ=40)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                         <  1:644766/64‑1 (MQ=40)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:2878327/38‑1 (MQ=14)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:1284265/38‑1 (MQ=14)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                >  1:2152503/1‑71 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                >  1:3704924/1‑71 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                >  1:3285568/1‑71 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                >  1:1693216/1‑71 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc              <  1:2041832/71‑1 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc              <  1:3625279/71‑1 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc              <  1:1479399/71‑1 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc            <  1:758109/71‑1 (MQ=255)

                                                    gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAg          >  1:3601/1‑71 (MQ=255)

                                                     ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAg          >  1:1621646/1‑70 (MQ=255)

                                                     ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc         >  1:1010757/1‑71 (MQ=255)

                                                     ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc         >  1:467730/1‑71 (MQ=255)

                                                      cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCt        <  1:396696/71‑1 (MQ=255)

                                                      cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCt        >  1:3849326/1‑71 (MQ=255)

                                                      cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCt        >  1:1342330/1‑71 (MQ=255)

                                                        cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                       >  1:2447867/1‑54 (MQ=255)

                                                        cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:3386166/1‑71 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt       <  1:2640541/69‑1 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt       <  1:2408686/69‑1 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      <  1:1569023/70‑1 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      <  1:1836567/70‑1 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      <  1:3545079/70‑1 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     >  1:524489/1‑71 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     >  1:1203112/1‑71 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     >  1:3174175/1‑71 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     >  1:88356/1‑71 (MQ=255)

                                                          caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt       <  1:2757883/68‑1 (MQ=255)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359815‑2359945

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