BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,359,876	C→T	27.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
*	Exported	2,359,876	0	C	T	27.4%	124.4 / 54.2	73	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (38/15); new base T (10/10); total (48/25)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCG  >  Exported/2359818‑2359942

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

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

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                        <  1:66456/71‑1 (MQ=255)

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

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

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

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

     aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                   >  1:1897381/1‑71 (MQ=34)

        gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                   <  1:1976499/68‑1 (MQ=35)

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

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

          tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                   >  1:502575/1‑66 (MQ=35)

          tACATATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                           >  1:148121/1‑58 (MQ=9)

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

             ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                              >  1:878157/1‑68 (MQ=35)

             ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                              >  1:963883/1‑68 (MQ=35)

               aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                  >  1:995360/1‑62 (MQ=35)

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

               aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                          >  1:1343740/1‑70 (MQ=255)

               aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                         >  1:617124/1‑71 (MQ=255)

                     gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                          <  1:383399/64‑1 (MQ=11)

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

                      cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                      >  1:1598157/1‑67 (MQ=35)

                      cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                  <  1:117699/71‑1 (MQ=34)

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

                        ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                           <  1:1477282/60‑1 (MQ=11)

                        ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                 >  1:300083/1‑70 (MQ=255)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                           >  1:1403569/1‑42 (MQ=32)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                       >  1:93021/1‑62 (MQ=35)

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

                           gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                     <  1:779038/47‑1 (MQ=38)

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

                             aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                         >  1:1672312/1‑57 (MQ=35)

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

                               ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:1823441/1‑65 (MQ=255)

                                tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                          >  1:161980/1‑37 (MQ=37)

                                tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                          >  1:494932/1‑37 (MQ=37)

                                tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                          >  1:444005/1‑37 (MQ=37)

                                tCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                           >  1:807939/1‑36 (MQ=33)

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

                                 cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                        >  1:1729478/1‑38 (MQ=40)

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

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

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

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                       <  1:1512889/68‑1 (MQ=255)

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

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

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg                     >  1:353240/1‑70 (MQ=255)

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

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

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

                                           cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:524157/1‑53 (MQ=35)

                                           cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                          >  1:520875/1‑58 (MQ=38)

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

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

                                           aaGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                               >  1:514411/2‑37 (MQ=255)

                                             gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                            <  1:2050315/38‑1 (MQ=40)

                                             gCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                       <  1:80719/43‑1 (MQ=14)

                                              cATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                              >  1:879528/1‑35 (MQ=40)

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

                                                tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCa        <  1:1081925/71‑1 (MQ=255)

                                                   cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                           <  1:2088661/49‑1 (MQ=255)

                                                   cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                          <  1:1620731/50‑1 (MQ=40)

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

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

                                                      acaGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:1038217/1‑42 (MQ=255)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                        >  1:433677/1‑49 (MQ=255)

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

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  <  1:1783151/71‑1 (MQ=255)

                                                      acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  <  1:789733/71‑1 (MQ=255)

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

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

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

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCG  >  Exported/2359818‑2359942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`