BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	38.8%	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	38.8%	108.9 / 139.7	121	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (39/35); new base T (30/17); total (69/52)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 126 total aligned reads displayed.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGC  >  minE/2359815‑2359946

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

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

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

accGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATCCCTCACAGCACACCTTca                                                               >  1:4624706/2‑71 (MQ=9)

 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               <  1:3428501/70‑1 (MQ=31)

 ccGCTTTATCCTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                              <  1:2277939/71‑1 (MQ=17)

    cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                             >  1:720232/1‑69 (MQ=255)

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

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

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

      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                          >  1:4423592/1‑70 (MQ=255)

      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                         >  1:3865647/1‑71 (MQ=255)

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                          <  1:2731700/69‑1 (MQ=255)

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTAACAGTACACCTTCACAGGCt                                                         >  1:2850615/1‑70 (MQ=11)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCGAGCATGCCTCACAGCACACCTTCACAGGc                                                          >  1:2954818/1‑68 (MQ=25)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                                        >  1:4585254/1‑69 (MQ=35)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTAc                                                      <  1:2655697/71‑1 (MQ=255)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                                          >  1:3705812/1‑66 (MQ=35)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                      <  1:1932400/70‑1 (MQ=35)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                      <  1:2321423/70‑1 (MQ=35)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                             >  1:4546977/1‑63 (MQ=255)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                             >  1:4403618/1‑63 (MQ=255)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                             >  1:3319936/1‑63 (MQ=255)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                                   <  1:807577/71‑1 (MQ=255)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                  >  1:1675785/1‑70 (MQ=35)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                  <  1:2099603/70‑1 (MQ=35)

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

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                               <  1:343959/69‑1 (MQ=35)

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

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                              <  1:192641/70‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                            >  1:2524948/1‑71 (MQ=34)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                            <  1:2577540/71‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                            <  1:490733/71‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                            <  1:3929110/71‑1 (MQ=255)

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

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

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

                        gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                  >  1:4256625/1‑60 (MQ=11)

                        gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                                        <  1:2680048/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                    <  1:2288325/71‑1 (MQ=34)

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

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                    >  1:4184093/1‑53 (MQ=35)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                               >  1:2419177/1‑42 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                               >  1:789320/1‑42 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                               >  1:978434/1‑42 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                            <  1:4721474/61‑1 (MQ=11)

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

                             cGCACTTCTGATACATCCAGCATACCTCACAGTACACCTTca                                                               >  1:2662720/1‑42 (MQ=11)

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

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

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

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                                >  1:3561457/1‑71 (MQ=34)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                                >  1:2900274/1‑71 (MQ=34)

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

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                                >  1:3206874/1‑71 (MQ=34)

                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                                <  1:2128450/70‑1 (MQ=35)

                                aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                               <  1:55815/71‑1 (MQ=255)

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

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

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                  <  1:1511020/66‑1 (MQ=35)

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

                                    cTGATACCTCCAGCATGCCTCCCAGCACACCTTCACAGGCTTACa                                                     >  1:3724186/1‑45 (MQ=255)

                                    cTGATACCTCCAGCATGCCTCACAGCACAGCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                            <  1:2825577/70‑1 (MQ=255)

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

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

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

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                                <  1:680589/63‑1 (MQ=35)

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

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

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

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                                   >  1:2237681/1‑39 (MQ=14)

                                             ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt                  >  1:1214617/1‑71 (MQ=255)

                                             ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt                  >  1:3904697/1‑71 (MQ=255)

                                             ccAGCAAGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt                  >  1:1775/1‑71 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                  >  1:64506/1‑38 (MQ=40)

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

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

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

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCACCTACCCAACAACACATAGTGTCGCTg                 >  1:4469297/1‑71 (MQ=255)

                                                 cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg              >  1:396123/1‑71 (MQ=255)

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

                                                       tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt        >  1:2279324/1‑71 (MQ=255)

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

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

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

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACACAACAACACATAGTGTCGCTGCCGCAGCTTCg      >  1:4255550/1‑71 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTca      <  1:4123923/67‑2 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg      <  1:4515706/67‑1 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg      <  1:4080304/67‑1 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGc  <  1:192350/71‑1 (MQ=255)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGC  >  minE/2359815‑2359946

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