BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	21.7%	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	21.7%	319.1 / 79.8	143	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (65/46); new base A (8/23); total (73/70)
Fisher's exact test for biased strand distribution p-value = 1.98e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 144 total aligned reads displayed.

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCA  >  minE/2359816‑2359936

 cGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                    <  1:1385908/70‑1 (MQ=255)

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

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

       aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                               <  1:600988/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                  tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                     <  1:2188137/68‑1 (MQ=255)

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

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

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

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

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

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

                    tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                <  1:1788556/71‑1 (MQ=34)

                    tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                <  1:3125323/71‑1 (MQ=34)

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

                        cATTCGCACTTCTGATACCTCCAGCATTCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:228072/57‑1 (MQ=255)

                        cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:848868/57‑1 (MQ=35)

                        cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                              >  1:4227763/1‑69 (MQ=35)

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

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

                        cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                 >  1:370716/1‑66 (MQ=11)

                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:1320190/1‑58 (MQ=35)

                          ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACCCCTTCACAGCCTTACAGaa                                        >  1:783567/1‑57 (MQ=11)

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

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

                          ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:300525/1‑57 (MQ=35)

                          ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:127849/1‑57 (MQ=35)

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

                          ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTGACAGGCt                                               <  1:4252523/50‑1 (MQ=11)

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

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

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

                            agCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                             >  1:2712470/2‑50 (MQ=38)

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

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

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

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                        <  1:1499544/70‑1 (MQ=255)

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                        <  1:1079635/70‑1 (MQ=255)

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                        >  1:2300759/1‑70 (MQ=255)

                                cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                >  1:1228740/1‑59 (MQ=35)

                                cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                 <  1:2846076/58‑1 (MQ=35)

                                cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                 <  1:2351209/58‑1 (MQ=35)

                                 ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        <  1:1529179/50‑1 (MQ=38)

                                 ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        <  1:582072/50‑1 (MQ=38)

                                  tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAAc                                       >  1:2005547/1‑50 (MQ=38)

                                  tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAAc                                       >  1:3503669/1‑50 (MQ=38)

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

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:2610784/1‑48 (MQ=38)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:1408856/1‑48 (MQ=38)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:1352238/1‑48 (MQ=38)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:2926736/1‑48 (MQ=38)

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

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

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

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

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

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

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

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

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTGCCCTACCCAACAACACATAg               <  1:963239/70‑1 (MQ=255)

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg               <  1:3145338/70‑1 (MQ=255)

                                      aTACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                <  1:3323050/53‑1 (MQ=11)

                                      aTACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                <  1:1851645/53‑1 (MQ=11)

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

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc           <  1:3559764/67‑1 (MQ=255)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc           <  1:311317/67‑1 (MQ=255)

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

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

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

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

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

                                              aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                         >  1:4037189/1‑52 (MQ=35)

                                              aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                         >  1:2370923/1‑52 (MQ=35)

                                              aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                         >  1:2588380/1‑52 (MQ=35)

                                               gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc     >  1:2938055/1‑71 (MQ=255)

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

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

                                                  tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc     >  1:131299/1‑68 (MQ=255)

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

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

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

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

                                                  tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGACg    >  1:1415893/1‑69 (MQ=255)

                                                   gCCTCACAGCACACCTTCACAGGCTTACAGCACGCTCCCCTAc                             >  1:2543338/1‑43 (MQ=255)

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCA  >  minE/2359816‑2359936

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