BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	38.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	38.5%	96.0 / 107.9	109	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (30/37); new base T (24/18); total (54/55)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 110 total aligned reads displayed.

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359818‑2359945

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

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                           >  1:3216557/1‑71 (MQ=255)

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                           >  1:2000973/1‑71 (MQ=255)

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                           >  1:4602559/1‑71 (MQ=255)

 cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              <  1:1836757/67‑1 (MQ=31)

 cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              <  1:5153493/67‑1 (MQ=31)

 cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                                          >  1:331885/1‑71 (MQ=34)

  tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                                         >  1:5623932/1‑71 (MQ=34)

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

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

    tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                        >  1:4457524/1‑70 (MQ=255)

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

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

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

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

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

         ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   >  1:5405943/1‑70 (MQ=35)

         ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   >  1:2651374/1‑70 (MQ=35)

         ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   >  1:2527219/1‑70 (MQ=35)

            cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               <  1:468211/71‑1 (MQ=34)

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

                tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               <  1:1179536/67‑1 (MQ=35)

                tGTCAGCATTCGCACTTCTGATACCTCCAGCATACGTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                            <  1:3012951/70‑1 (MQ=11)

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

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

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

                  tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      <  1:3607846/58‑1 (MQ=35)

                  tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                         <  1:5434838/71‑1 (MQ=255)

                   cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               >  1:1059741/1‑64 (MQ=35)

                   cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               >  1:786363/1‑64 (MQ=35)

                   cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               >  1:4955679/1‑64 (MQ=35)

                    aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                        <  1:1678427/70‑1 (MQ=35)

                    aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                        <  1:3078047/70‑1 (MQ=35)

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

                    aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTAACAGGCTTACAGAACGCTCCCCTa                                       <  1:1551669/71‑1 (MQ=11)

                     gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   <  1:3139555/58‑1 (MQ=11)

                      cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                        >  1:1774057/1‑68 (MQ=255)

                        ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:2011287/1‑68 (MQ=35)

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

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

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

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

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

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

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGACCGCTCCCCTACCcaa                                  >  1:2506401/1‑70 (MQ=11)

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

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

                           gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                                <  1:2474110/71‑1 (MQ=255)

                           gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:1425669/70‑1 (MQ=255)

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

                            cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                   <  1:3668574/67‑1 (MQ=255)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:2378503/49‑1 (MQ=38)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:4719608/49‑1 (MQ=38)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:1491225/49‑1 (MQ=38)

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

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                             >  1:720663/1‑71 (MQ=35)

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

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

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

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

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                     <  1:3222575/46‑1 (MQ=38)

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                         <  1:4416140/58‑1 (MQ=35)

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            <  1:1190031/55‑1 (MQ=35)

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                             >  1:5214529/1‑70 (MQ=37)

                               ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                              <  1:1609437/69‑1 (MQ=255)

                                 cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            <  1:1627755/53‑1 (MQ=35)

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

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                          <  1:1393065/68‑1 (MQ=40)

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                          <  1:1131736/68‑1 (MQ=40)

                                      aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:827228/1‑61 (MQ=255)

                                      aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:2301972/1‑61 (MQ=255)

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

                                        cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg                   <  1:4109123/71‑1 (MQ=255)

                                         tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                              <  1:3853896/59‑1 (MQ=35)

                                         tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                              <  1:1947825/59‑1 (MQ=35)

                                         tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc                  >  1:170487/1‑71 (MQ=255)

                                         tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:4084914/1‑35 (MQ=37)

                                         tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:5323336/1‑35 (MQ=37)

                                         tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:2651291/1‑35 (MQ=37)

                                          ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:4063043/1‑61 (MQ=40)

                                          ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:5630974/1‑61 (MQ=40)

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

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

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

                                                 accTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    <  1:5269815/44‑1 (MQ=255)

                                                  ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc         <  1:2498512/71‑1 (MQ=255)

                                                  ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc         <  1:5274295/71‑1 (MQ=255)

                                                  ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc         <  1:1258950/71‑1 (MQ=255)

                                                   cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:4439936/1‑45 (MQ=255)

                                                   cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:1419491/1‑45 (MQ=255)

                                                    tCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:4770333/45‑1 (MQ=255)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTcccc                                         >  1:3275099/1‑35 (MQ=40)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTcccc                                         >  1:3378094/1‑35 (MQ=40)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTcccc                                         >  1:2286564/1‑35 (MQ=40)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg                      <  1:1282251/54‑1 (MQ=255)

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

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

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

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

                                                          cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTg  >  1:746457/1‑70 (MQ=255)

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359818‑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: `‑`