BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	42.9%	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.9%	126.9 / 144.0	112	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (26/38); new base A (32/16); total (58/54)
Fisher's exact test for biased strand distribution p-value = 7.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 113 total aligned reads displayed.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAG  >  minE/2359815‑2359937

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                      >  1:982880/1‑71 (MQ=31)

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

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                      >  1:1733050/1‑71 (MQ=31)

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

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

cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                      >  1:1395476/1‑71 (MQ=31)

 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                      >  1:2238087/1‑70 (MQ=31)

   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTACACAg                                                   <  1:3194566/71‑1 (MQ=21)

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

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

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

       tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACTTCACAGTACACCTTCACAGGCt                                                >  1:3237563/1‑70 (MQ=11)

        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGttt                                               <  1:1648548/70‑1 (MQ=11)

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

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

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

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

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                           <  1:1474561/70‑1 (MQ=35)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                           <  1:3828601/70‑1 (MQ=35)

            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                           <  1:1103909/70‑1 (MQ=35)

             tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                         >  1:3583028/1‑71 (MQ=255)

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

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                     >  1:3135225/1‑54 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                      >  1:580225/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                      >  1:3166627/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                      >  1:2779957/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                      >  1:3500572/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                      >  1:1068439/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                      >  1:3739743/1‑53 (MQ=11)

                  aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                      >  1:551030/1‑53 (MQ=11)

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

                  aTGTCAGCATTCGAACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                      <  1:1031413/69‑1 (MQ=25)

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

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

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

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                            >  1:3513471/1‑59 (MQ=35)

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                            >  1:249889/1‑59 (MQ=35)

                      cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                 >  1:2916741/1‑70 (MQ=35)

                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCACTa                               >  1:822378/1‑71 (MQ=11)

                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                              <  1:2800807/71‑1 (MQ=34)

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

                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                             >  1:2833997/1‑71 (MQ=34)

                          aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                             >  1:1478719/1‑70 (MQ=255)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                              <  1:3282223/52‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                              <  1:2955148/52‑1 (MQ=11)

                           ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                           <  1:3585397/55‑1 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                              >  1:1344490/1‑51 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                              >  1:1317250/1‑51 (MQ=11)

                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                              >  1:1764836/1‑51 (MQ=11)

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

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

                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACAcaa                          >  1:1305419/1‑70 (MQ=11)

                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                         >  1:797339/1‑70 (MQ=35)

                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                         >  1:827773/1‑70 (MQ=35)

                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                      <  1:2952889/71‑1 (MQ=34)

                                 cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                     >  1:1334627/1‑70 (MQ=255)

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

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                           <  1:551332/48‑1 (MQ=11)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                           <  1:3484645/48‑1 (MQ=11)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                           <  1:1826415/48‑1 (MQ=11)

                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                           <  1:673706/48‑1 (MQ=11)

                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                              >  1:2673193/1‑43 (MQ=38)

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

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

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

                                      gATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                 >  1:1772659/1‑38 (MQ=14)

                                      gATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                              >  1:1752956/1‑41 (MQ=14)

                                      gATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                         >  1:1010879/1‑46 (MQ=11)

                                      gATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                         >  1:657911/1‑46 (MQ=11)

                                      gATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                         >  1:1961719/1‑46 (MQ=11)

                                      gATACCTCCAGCATACCTCACAGTACAACTTCACAGGCTTACAGaa                                         >  1:3446275/1‑46 (MQ=11)

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

                                         aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                  >  1:1091262/1‑50 (MQ=11)

                                         aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                  >  1:3076725/1‑50 (MQ=11)

                                          ccTCCATCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc            <  1:1896869/71‑1 (MQ=255)

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

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

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

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

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

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                       <  1:2894417/58‑1 (MQ=35)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                       <  1:2375323/58‑1 (MQ=35)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                 <  1:1032138/64‑1 (MQ=40)

                                            tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                 <  1:644766/64‑1 (MQ=40)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                           <  1:1284265/38‑1 (MQ=14)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                           <  1:2878327/38‑1 (MQ=14)

                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                           <  1:1371655/38‑1 (MQ=14)

                                             ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACTCAACAACACATAGTGTCGCt         >  1:942745/1‑71 (MQ=255)

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

                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt         >  1:3147321/1‑70 (MQ=255)

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

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

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

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

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc      <  1:551823/71‑1 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc      <  1:2041832/71‑1 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc      <  1:3625279/71‑1 (MQ=255)

                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc      <  1:1479399/71‑1 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc    <  1:758109/71‑1 (MQ=255)

                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc    <  1:31929/71‑1 (MQ=255)

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAG  >  minE/2359815‑2359937

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