BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	42.4%	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	42.4%	87.2 / 126.5	99	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (36/21); new base T (20/22); total (56/43)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 101 total aligned reads displayed.

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359818‑2359945

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

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

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

  tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                          >  1:306494/1‑70 (MQ=255)

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

     aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                        <  1:2392173/69‑1 (MQ=255)

      tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      <  1:2522956/70‑1 (MQ=35)

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

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

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

        gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                                    <  1:1234611/70‑1 (MQ=255)

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

         ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                   >  1:1860133/1‑70 (MQ=255)

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

          tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:1908192/1‑66 (MQ=255)

          tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:1770748/1‑66 (MQ=255)

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

           aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:1304051/1‑70 (MQ=255)

               aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:927730/1‑53 (MQ=31)

               aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTTACAGTACACCTTCACAGGCTTACAGAACGc                                              <  1:309652/69‑1 (MQ=11)

               aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                           >  1:1086804/1‑56 (MQ=11)

               aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                        >  1:1642377/1‑59 (MQ=11)

               aTGTCAGCATTCCCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                            >  1:305887/1‑71 (MQ=11)

                tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                           >  1:2438209/1‑71 (MQ=255)

                  tCAGCATTCGCACTTCTGATACCTCCAGCATGCGTCACAGCACACCTTCACAGGCTTACAGa                                                  <  1:2934585/62‑1 (MQ=25)

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

                   cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 >  1:1102336/1‑62 (MQ=35)

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

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

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

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

                     gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      <  1:2199996/71‑1 (MQ=255)

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

                      cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcaca                                                            >  1:586703/1‑48 (MQ=38)

                      cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                     >  1:2725492/1‑55 (MQ=35)

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

                      cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:497636/1‑46 (MQ=11)

                      cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:308094/1‑71 (MQ=255)

                      cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:2931378/1‑71 (MQ=255)

                       aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    <  1:1653236/71‑1 (MQ=34)

                        ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                           <  1:1139034/47‑1 (MQ=11)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      <  1:1548328/50‑1 (MQ=38)

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

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:528784/1‑50 (MQ=11)

                           gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                                       <  1:1225384/48‑1 (MQ=11)

                             aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      <  1:1162471/47‑1 (MQ=38)

                              cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 <  1:33267/51‑1 (MQ=11)

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

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 <  1:1713944/50‑1 (MQ=38)

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

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

                               ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                                    <  1:426962/47‑1 (MQ=11)

                               ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                            <  1:1432454/55‑1 (MQ=11)

                               ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                            <  1:108771/55‑1 (MQ=11)

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

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

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

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

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

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

                                  tGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                               <  1:1927841/49‑1 (MQ=11)

                                   gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        <  1:2790287/39‑1 (MQ=40)

                                   gATACCTCCAGCATGCCTCACAGCAAACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg                        >  1:313576/1‑71 (MQ=255)

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

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

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                       <  1:398862/71‑1 (MQ=255)

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                       <  1:1844809/71‑1 (MQ=255)

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

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                       <  1:1723937/71‑1 (MQ=255)

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

                                      aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:630096/1‑43 (MQ=14)

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

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

                                         tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:2862434/1‑40 (MQ=14)

                                            aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    >  1:835440/1‑50 (MQ=38)

                                            aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc               <  1:2711943/71‑1 (MQ=255)

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

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

                                                  ccTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                          >  1:1087820/1‑38 (MQ=255)

                                                   cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    <  1:767326/43‑1 (MQ=255)

                                                   cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     <  1:2771791/42‑1 (MQ=255)

                                                   cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 >  1:2579011/1‑46 (MQ=38)

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

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

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

                                                     cacaGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:2951166/1‑43 (MQ=255)

                                                     cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt       >  1:2922341/1‑70 (MQ=255)

                                                     cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt       >  1:2398022/1‑70 (MQ=255)

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

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

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:487766/1‑38 (MQ=40)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:2600808/1‑38 (MQ=40)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:854628/1‑42 (MQ=38)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg                        >  1:470640/1‑52 (MQ=255)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAg          >  1:1173297/1‑66 (MQ=255)

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

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

                                                        aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      <  1:2416691/68‑1 (MQ=255)

                                                          cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                            >  1:627228/1‑44 (MQ=255)

                                                          cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTg  >  1:205174/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: `‑`