BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	42.0%	noncoding (1228/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,867	0	G	A	42.0%	75.5 / 141.1	113	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (33/32); new base A (30/17); total (64/49)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 113 total aligned reads displayed.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  minE/2359815‑2359934

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   <  1:293808/71‑1 (MQ=31)

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   <  1:1898262/71‑1 (MQ=31)

accGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATCCCTCACAGCACACCTTca                                                   >  1:4624706/2‑71 (MQ=9)

 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   <  1:3428501/70‑1 (MQ=31)

 ccGCTTTATCCTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                  <  1:2277939/71‑1 (MQ=17)

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

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

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

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

      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                              >  1:4423592/1‑70 (MQ=255)

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

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCGAGCATGCCTCACAGCACACCTTCACAGGc                                              >  1:2954818/1‑68 (MQ=25)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                            >  1:4585254/1‑69 (MQ=35)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                             <  1:803164/68‑1 (MQ=255)

         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTAc                                          <  1:2655697/71‑1 (MQ=255)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                              >  1:3705812/1‑66 (MQ=35)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                          <  1:1932400/70‑1 (MQ=35)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                          <  1:2321423/70‑1 (MQ=35)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                 >  1:4403618/1‑63 (MQ=255)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                 >  1:3319936/1‑63 (MQ=255)

          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                 >  1:4546977/1‑63 (MQ=255)

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

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

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                      >  1:1675785/1‑70 (MQ=35)

              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAAc                                     >  1:2908960/1‑71 (MQ=34)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                   <  1:708691/69‑1 (MQ=35)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                   <  1:343959/69‑1 (MQ=35)

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

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                  <  1:192641/70‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                >  1:2524948/1‑71 (MQ=34)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                <  1:490733/71‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                <  1:3929110/71‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                <  1:2577540/71‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 <  1:4001728/70‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 <  1:351761/70‑1 (MQ=255)

                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 <  1:656003/70‑1 (MQ=255)

                        gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      >  1:4256625/1‑60 (MQ=11)

                        gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                            <  1:2680048/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:2957546/50‑1 (MQ=38)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:3751875/50‑1 (MQ=38)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:3990786/50‑1 (MQ=38)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:1088287/50‑1 (MQ=38)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:1735435/50‑1 (MQ=38)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:3099804/50‑1 (MQ=38)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:2974978/50‑1 (MQ=38)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:4564270/50‑1 (MQ=38)

                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                        <  1:2288325/71‑1 (MQ=34)

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

                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           >  1:4582491/1‑50 (MQ=38)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   >  1:789320/1‑42 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   >  1:978434/1‑42 (MQ=11)

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   >  1:2419177/1‑42 (MQ=11)

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

                             cGCACTTCTGATACATCCAGCATACCTCACAGTACACCTTca                                                   >  1:2662720/1‑42 (MQ=11)

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

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

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

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

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    >  1:1889380/1‑71 (MQ=34)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    >  1:2715963/1‑71 (MQ=34)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    >  1:2900274/1‑71 (MQ=34)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    >  1:3206874/1‑71 (MQ=34)

                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    >  1:3561457/1‑71 (MQ=34)

                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    <  1:2128450/70‑1 (MQ=35)

                                aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                   <  1:55815/71‑1 (MQ=255)

                                 cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                  <  1:888462/71‑2 (MQ=255)

                                 cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                  <  1:4667499/71‑2 (MQ=255)

                                  ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                      <  1:1511020/66‑1 (MQ=35)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                  >  1:3516219/1‑69 (MQ=11)

                                    cTGATACCTCCAGCATGCCTCCCAGCACACCTTCACAGGCTTACa                                         >  1:3724186/1‑45 (MQ=255)

                                    cTGATACCTCCAGCATGCCTCACAGCACAGCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                <  1:2825577/70‑1 (MQ=255)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                         >  1:1879628/1‑45 (MQ=38)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                         >  1:4397127/1‑45 (MQ=38)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                         >  1:2380907/1‑45 (MQ=38)

                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    <  1:680589/63‑1 (MQ=35)

                                        tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg           <  1:3423752/71‑1 (MQ=255)

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

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

                                           cTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                       >  1:2482750/1‑40 (MQ=14)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                       >  1:2237681/1‑39 (MQ=14)

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

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

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

                                             ccAGCAAGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt      >  1:1775/1‑71 (MQ=255)

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

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

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

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

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCACCTACCCAACAACACATAGTGTCGCTg     >  1:4469297/1‑71 (MQ=255)

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

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg  >  1:4000935/1‑70 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg  >  1:3590141/1‑70 (MQ=255)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  minE/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: `‑`