BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	24.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	24.1%	237.3 / 69.2	108	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (45/37); new base T (15/11); total (60/48)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 109 total aligned reads displayed.

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359816‑2359945

ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              <  1:1219821/70‑1 (MQ=255)

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

        tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      >  1:1920024/1‑70 (MQ=255)

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

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

           ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                                  <  1:2912684/71‑1 (MQ=34)

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

             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    >  1:1947817/1‑67 (MQ=35)

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

             aCTTATGTCAGAATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    >  1:1806242/1‑67 (MQ=25)

               ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                        <  1:1864834/61‑1 (MQ=255)

                  tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                           <  1:1015999/71‑1 (MQ=34)

                  tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            <  1:1264742/70‑1 (MQ=35)

                    tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                          <  1:945879/70‑1 (MQ=35)

                    tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTAACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                         <  1:2477947/71‑1 (MQ=21)

                    tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                         >  1:95498/1‑71 (MQ=255)

                      aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    >  1:599781/1‑58 (MQ=35)

                      aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                                    >  1:1049703/1‑58 (MQ=35)

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

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

                      aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                                       >  1:2996662/1‑71 (MQ=255)

                       gCATTCGCACTTCTGATACCTCC‑‑CATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    <  1:3304396/71‑1 (MQ=11)

                        cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:2696427/1‑46 (MQ=32)

                        cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:2930220/1‑46 (MQ=32)

                        cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:1399405/1‑46 (MQ=32)

                        cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      >  1:2549612/1‑70 (MQ=255)

                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:3256678/1‑70 (MQ=35)

                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:205052/1‑70 (MQ=35)

                          ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    >  1:2420337/1‑70 (MQ=35)

                          ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                    >  1:2221000/1‑70 (MQ=35)

                           tCGCACTTCTGATACCTCCAGCATGCCTCCCAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:1482346/1‑71 (MQ=255)

                            cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                           >  1:3148184/1‑61 (MQ=11)

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

                             gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                      <  1:2095182/65‑1 (MQ=35)

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      <  1:331304/49‑1 (MQ=11)

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      <  1:2027009/49‑1 (MQ=11)

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      <  1:342486/49‑1 (MQ=11)

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      <  1:2103543/49‑1 (MQ=11)

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      <  1:1308939/49‑1 (MQ=11)

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                      <  1:107829/49‑1 (MQ=11)

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

                               aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:382897/1‑70 (MQ=35)

                               aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:1490019/1‑70 (MQ=35)

                                cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACacaacaaca                             >  1:1840053/1‑71 (MQ=25)

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

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

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

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

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

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

                                   cTGATACCTCCGGCATGCCTCACAGCACACCTTCACAGGCt                                                        <  1:2216411/41‑1 (MQ=255)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        <  1:2068936/41‑1 (MQ=38)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        <  1:350818/41‑1 (MQ=38)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        <  1:2568887/41‑1 (MQ=38)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               >  1:1292132/1‑50 (MQ=38)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               >  1:577036/1‑50 (MQ=38)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:667999/64‑1 (MQ=35)

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

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

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

                                   cTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                          >  1:2408508/1‑39 (MQ=14)

                                   cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:1197054/1‑40 (MQ=14)

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

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

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

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

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

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

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

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

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

                                           tCCAGCATGCCTAACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc                  >  1:1405040/1‑71 (MQ=255)

                                           tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                          >  1:3177353/1‑47 (MQ=11)

                                           tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                          >  1:3303068/1‑47 (MQ=11)

                                           tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                          >  1:1183944/1‑47 (MQ=11)

                                           tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:709563/1‑58 (MQ=11)

                                           tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               >  1:159023/1‑58 (MQ=11)

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

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

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

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

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

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

                                                cATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                            <  1:64735/56‑3 (MQ=11)

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

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

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

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

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

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

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

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

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

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

                                                         caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCt        <  1:1696432/67‑1 (MQ=255)

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

                                                         caGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGAGTCGCTGCCGCAGCTTCgg    >  1:1409340/1‑71 (MQ=255)

                                                           gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTTTCGCTGCCg             >  1:536702/1‑60 (MQ=255)

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

                                                            cacaCCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGt   >  1:143002/1‑69 (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: `‑`