BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	20.0%	noncoding (1219/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,876	0	C	T	20.0%	371.0 / 79.6	155	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (74/50); new base T (8/23); total (82/73)
Fisher's exact test for biased strand distribution p-value = 1.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 155 total aligned reads displayed.

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCG  >  minE/2359816‑2359942

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

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

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

      tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTAACAGGCt                                                     <  1:2054615/70‑1 (MQ=11)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                          ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                     <  1:1824158/50‑1 (MQ=11)

                          ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                     <  1:797005/50‑1 (MQ=11)

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

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

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

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

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

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

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

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

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

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

                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCACCCCTACCcaac                              >  1:61513/1‑70 (MQ=11)

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

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

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

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

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

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

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

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

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

                                  tCTGATACCTACAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     >  1:4237762/1‑58 (MQ=255)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                              >  1:2285484/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)

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

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

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                          >  1:2146066/1‑68 (MQ=37)

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

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

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

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

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

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

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

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

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

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

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                      >  1:150587/1‑46 (MQ=38)

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

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

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

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

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

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

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

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

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

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

                                               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:2057843/1‑69 (MQ=255)

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

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

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

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

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

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

                                                    ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAg       <  1:2200890/70‑1 (MQ=255)

                                                    ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc      >  1:3063202/1‑71 (MQ=255)

                                                    ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGc      >  1:2243331/1‑71 (MQ=255)

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

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

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

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