BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,359,867	G→A	31.3%	noncoding (1228/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,867	0	G	A	31.3%	139.1 / 72.8	83	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (28/29); new base A (9/17); total (37/46)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  Exported/2359815‑2359934

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:2408087/70‑1 (MQ=35)

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             >  1:2470685/1‑70 (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:217311/1‑43 (MQ=32)

                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   >  1:2771094/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:1624995/46‑1 (MQ=11)

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

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

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

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

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                             <  1:420758/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:2738737/1‑59 (MQ=35)

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

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

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

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

                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt          <  1:2368637/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:2283097/1‑52 (MQ=11)

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

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

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

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

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

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

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

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

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

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

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

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

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  Exported/2359815‑2359934

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: `‑`