BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,359,876	0	C	T	28.0%	173.7 / 71.6	93	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (37/29); new base T (9/17); total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGC  >  Exported/2359815‑2359946

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                               >  1:1046526/1‑71 (MQ=9)

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

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

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

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                         <  1:2966414/70‑1 (MQ=35)

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

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                         <  1:2408087/70‑1 (MQ=35)

       tATCGTTACTTATGTCAGCACTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                         <  1:2020251/70‑1 (MQ=25)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                                        <  1:326997/69‑1 (MQ=35)

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

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

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

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

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                                 >  1:2504036/1‑71 (MQ=255)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                                 >  1:2676685/1‑71 (MQ=255)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                                 >  1:2766405/1‑71 (MQ=255)

               cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                  <  1:998016/69‑1 (MQ=35)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGCACGCTc                                             >  1:2641624/1‑71 (MQ=21)

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

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                       >  1:1509450/1‑61 (MQ=255)

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

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

                     tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                          <  1:2612728/71‑1 (MQ=34)

                     tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                         <  1:2958273/56‑1 (MQ=11)

                     tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                         <  1:1993999/56‑1 (MQ=11)

                     tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                         <  1:675534/56‑1 (MQ=11)

                      cAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                           <  1:2523171/69‑1 (MQ=255)

                      cAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                           <  1:1109508/69‑1 (MQ=255)

                       aGCATTCGCACTTCTGATACCTCCAGCTTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                        >  1:375764/1‑71 (MQ=255)

                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                        >  1:1782706/1‑71 (MQ=34)

                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                    <  1:992068/59‑1 (MQ=11)

                         cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                    <  1:1229615/57‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTAc                                                      <  1:405087/53‑1 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               >  1:2771094/1‑43 (MQ=32)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               >  1:217311/1‑43 (MQ=32)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               >  1:1291736/1‑43 (MQ=32)

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccaac                                   <  1:2823863/70‑4 (MQ=25)

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

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                   <  1:1262444/70‑1 (MQ=255)

                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                               >  1:1624030/1‑41 (MQ=32)

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

                              gCACTTCAGATAACTCCAGCATGACTCACAGCACACCTTca                                                               >  1:2447905/1‑41 (MQ=11)

                                 cTTCTGATACCTCCAGCATGCCTCACAGGACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                              <  1:1000026/71‑1 (MQ=25)

                                 cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                       <  1:130573/46‑1 (MQ=11)

                                 cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                       <  1:1624995/46‑1 (MQ=11)

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

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

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

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

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                         <  1:2591330/59‑1 (MQ=11)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                         <  1:1122216/59‑1 (MQ=11)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                              >  1:1759987/1‑68 (MQ=37)

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

                                     tGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                  <  1:2782922/47‑1 (MQ=11)

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                         <  1:1687080/54‑1 (MQ=35)

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                               >  1:2231490/1‑64 (MQ=35)

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

                                       aTACCTCCAGCATGCCCCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                    <  1:873169/59‑1 (MQ=255)

                                        tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                        >  1:424156/1‑70 (MQ=255)

                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                  >  1:3017903/1‑59 (MQ=35)

                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                  >  1:285310/1‑59 (MQ=35)

                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                  >  1:2738737/1‑59 (MQ=35)

                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                  >  1:1524632/1‑59 (MQ=35)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                      <  1:2368637/70‑1 (MQ=255)

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                      <  1:2750759/70‑1 (MQ=255)

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

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                       >  1:369239/1‑51 (MQ=35)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                            <  1:599108/62‑1 (MQ=40)

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

                                             ccAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                     >  1:1844681/1‑52 (MQ=11)

                                             ccAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                     >  1:1803726/1‑52 (MQ=11)

                                             ccAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                     >  1:2283097/1‑52 (MQ=11)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt                  <  1:802483/70‑1 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt                  <  1:1517082/70‑1 (MQ=255)

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt                  <  1:1438000/70‑1 (MQ=255)

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

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

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

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

                                                 cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                               >  1:2950088/1‑38 (MQ=40)

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

                                                 cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg              <  1:44345/71‑1 (MQ=255)

                                                   tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg              >  1:2606261/1‑69 (MQ=255)

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

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

                                                      cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg              >  1:1582709/1‑66 (MQ=255)

                                                      cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAg           >  1:986350/1‑69 (MQ=255)

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

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg                       >  1:686832/1‑54 (MQ=255)

                                                         acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg                       >  1:675006/1‑54 (MQ=255)

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

                                                            gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                            >  1:2059822/1‑46 (MQ=255)

                                                            gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTg   >  1:2317072/1‑71 (MQ=255)

                                                             cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCtt                 <  1:668622/56‑2 (MQ=255)

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

                                                             cacaCCTTAACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGc  <  1:922163/71‑1 (MQ=255)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTGC  >  Exported/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: `‑`