BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	41.1%	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	41.1%	131.2 / 166.7	146	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (42/44); new base T (33/27); total (75/71)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

Only 100 of 147 total aligned reads displayed.

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359816‑2359945

ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:1477468/1‑70 (MQ=31)

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

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

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

  gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                           >  1:2124850/1‑71 (MQ=34)

   cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCGAGCATACCTCACAGTACACCTTCACAgg                                                          >  1:1392461/1‑71 (MQ=11)

   cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                          >  1:692578/1‑71 (MQ=255)

    tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:911028/1‑71 (MQ=255)

    tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCAAGCATGCCTCACAGCACACCTTCACAgg                                                          >  1:3451088/1‑70 (MQ=25)

      tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                        <  1:3176459/70‑1 (MQ=255)

       aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         <  1:2520470/68‑1 (MQ=255)

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

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

         cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:1092255/1‑66 (MQ=255)

         cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:3078911/1‑66 (MQ=255)

          gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                           <  1:3822841/63‑1 (MQ=35)

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

            tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 >  1:507290/1‑71 (MQ=34)

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

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

             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                                  >  1:1792964/1‑69 (MQ=255)

              cTTATTTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                                <  1:608686/70‑1 (MQ=11)

              cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                               >  1:2935102/1‑71 (MQ=255)

                tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                             >  1:3637354/1‑71 (MQ=255)

                 aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                                          >  1:2850442/1‑57 (MQ=35)

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

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

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

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

                 aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                              <  1:1500201/69‑1 (MQ=255)

                 aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                              <  1:3350054/69‑1 (MQ=255)

                      aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                                       <  1:3128590/71‑1 (MQ=255)

                       gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              <  1:2001499/47‑1 (MQ=11)

                       gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              <  1:1265550/47‑1 (MQ=11)

                       gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              <  1:1885316/47‑1 (MQ=11)

                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            <  1:2373821/63‑1 (MQ=35)

                         aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:186656/1‑53 (MQ=11)

                         aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:528521/1‑53 (MQ=11)

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

                          ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        <  1:1702775/50‑1 (MQ=38)

                           tCGCACTTCAGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        <  1:1685869/49‑1 (MQ=255)

                            cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      >  1:2996997/1‑50 (MQ=38)

                            cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:3887420/1‑55 (MQ=11)

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

                             gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                                <  1:2799559/71‑1 (MQ=34)

                             gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:848145/70‑1 (MQ=35)

                              cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                   <  1:1340900/67‑1 (MQ=35)

                              cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:2169384/1‑71 (MQ=255)

                              cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:2077832/1‑71 (MQ=255)

                                cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:1440228/1‑46 (MQ=11)

                                cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                             <  1:2903473/71‑2 (MQ=255)

                                cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                             <  1:1051206/71‑2 (MQ=255)

                                 ttCTGATACCTCCAGCATGCCTCACAGCACACCTTcaca                                                            <  1:808932/39‑1 (MQ=40)

                                 ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                             >  1:1573399/1‑69 (MQ=11)

                                  tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                            >  1:333701/1‑70 (MQ=39)

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

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACTGAACGCTcc                                           <  1:2674783/54‑1 (MQ=255)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                       >  1:1054903/1‑58 (MQ=35)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                       <  1:1193091/58‑1 (MQ=35)

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

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           <  1:865284/70‑1 (MQ=40)

                                    tGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                         <  1:363063/55‑1 (MQ=35)

                                    tGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                             <  1:1923532/67‑1 (MQ=37)

                                     gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 <  1:2844535/46‑1 (MQ=38)

                                     gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:1235587/1‑64 (MQ=35)

                                     gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:1571656/1‑64 (MQ=35)

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:1040399/1‑56 (MQ=35)

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                         <  1:2651780/69‑1 (MQ=40)

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                         <  1:3578356/69‑1 (MQ=40)

                                       tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 >  1:2184248/1‑44 (MQ=38)

                                       tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg                      <  1:2119030/71‑1 (MQ=255)

                                       tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg                      >  1:3067228/1‑71 (MQ=255)

                                       tACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:822565/1‑55 (MQ=11)

                                        aCCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:438091/1‑35 (MQ=14)

                                         ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:1090963/58‑1 (MQ=35)

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

                                          cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:2566768/1‑56 (MQ=35)

                                          cTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 >  1:1916586/1‑57 (MQ=11)

                                          cTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 >  1:14631/1‑57 (MQ=11)

                                            ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               >  1:1344074/1‑41 (MQ=38)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                          >  1:1386564/1‑61 (MQ=40)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg                      >  1:3036903/1‑65 (MQ=255)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                    >  1:2250758/1‑67 (MQ=255)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc                  <  1:427270/69‑1 (MQ=255)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                <  1:79979/71‑1 (MQ=255)

                                              aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                >  1:2209717/1‑70 (MQ=255)

                                              aGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                         <  1:3605521/45‑1 (MQ=11)

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

                                               gCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                        <  1:3495681/45‑1 (MQ=11)

                                                 aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg             >  1:240438/1‑70 (MQ=255)

                                                     cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                                <  1:3283897/47‑1 (MQ=255)

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

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

                                                      tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                        <  1:3827521/38‑1 (MQ=40)

                                                       cacaGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                              >  1:3364657/1‑47 (MQ=255)

                                                       cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg             >  1:2815191/1‑64 (MQ=255)

                                                        acaGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:1522815/43‑1 (MQ=255)

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

                                                         caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCgg    >  1:1216643/1‑71 (MQ=255)

                                                          aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGt   <  1:415577/71‑1 (MQ=255)

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359816‑2359945

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