BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	38.5%	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	38.5%	96.9 / 91.5	83	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (37/14); new base T (20/12); total (57/26)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGT  >  minE/2359818‑2359944

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                          >  1:1882702/1‑71 (MQ=34)

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGACTCACAGCACACCTTCACAg                                                          >  1:365894/1‑71 (MQ=21)

   ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                       >  1:1860690/1‑71 (MQ=255)

    tATCTTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                       >  1:1767217/1‑70 (MQ=25)

    tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                       <  1:681136/70‑1 (MQ=255)

     aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                                      <  1:145267/70‑1 (MQ=35)

     aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                     >  1:1125325/1‑71 (MQ=34)

      tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                     >  1:1385299/1‑70 (MQ=255)

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

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

          tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                >  1:1194540/1‑71 (MQ=34)

          tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                >  1:724806/1‑71 (MQ=34)

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

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

            cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                              <  1:494623/71‑1 (MQ=255)

              tATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                     >  1:1567782/1‑62 (MQ=35)

                  tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                        >  1:1224468/1‑71 (MQ=255)

                   cAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                >  1:1930728/1‑62 (MQ=255)

                   cAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                        >  1:112269/1‑70 (MQ=255)

                    aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcaca                                                           >  1:888889/1‑50 (MQ=38)

                    aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                             >  1:237522/1‑48 (MQ=32)

                    aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                       <  1:883027/70‑1 (MQ=35)

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

                        ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                          <  1:1464258/47‑1 (MQ=38)

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

                        ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                    >  1:854320/1‑69 (MQ=35)

                        ttCGCACTTCTGATACCTACAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                   >  1:1609699/1‑70 (MQ=255)

                         tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                  >  1:2168208/1‑70 (MQ=255)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                                         >  1:1092052/1‑46 (MQ=38)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                                        >  1:139268/1‑47 (MQ=38)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                                        >  1:178057/1‑47 (MQ=38)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                 >  1:1730691/1‑70 (MQ=35)

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

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

                           gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                     >  1:441747/1‑65 (MQ=255)

                             aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                <  1:567237/68‑1 (MQ=255)

                             aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                <  1:1891187/68‑1 (MQ=255)

                             aCTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                <  1:1271579/68‑1 (MQ=255)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCt                                            <  1:852372/55‑1 (MQ=35)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                            <  1:1784140/71‑1 (MQ=35)

                              cTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                             >  1:268285/1‑38 (MQ=11)

                              cTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                             >  1:1619419/1‑38 (MQ=11)

                              cTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                             >  1:253619/1‑38 (MQ=11)

                               ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                           <  1:66421/55‑1 (MQ=11)

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

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

                                tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTA‑CCAACAAcaca                          >  1:414342/1‑70 (MQ=255)

                                tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                           >  1:508960/1‑68 (MQ=11)

                                 cTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                             >  1:1906932/1‑35 (MQ=33)

                                 cTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                             >  1:68344/1‑35 (MQ=33)

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

                                 cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                          <  1:118776/70‑1 (MQ=40)

                                 cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                        >  1:1471162/1‑40 (MQ=14)

                                 cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                >  1:2014475/1‑48 (MQ=11)

                                 cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                >  1:2012143/1‑48 (MQ=11)

                                 cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                           <  1:945315/69‑3 (MQ=11)

                                 cTGATACCTCCAGCATACCTCACAGTACACCATCACAGGc                                                        >  1:959932/1‑40 (MQ=14)

                                 cTGATACCTCCAGCATACCTCACACTACACCTTCACAGGCTTACACaa                                                >  1:1596364/1‑48 (MQ=11)

                                  tGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                        <  1:1172456/55‑1 (MQ=11)

                                   gATACCTCCAGCATACCTCACAGTACACCTTcaca                                                           >  1:1576925/1‑35 (MQ=14)

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                      <  1:1578604/71‑1 (MQ=255)

                                     tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                      >  1:2086621/1‑54 (MQ=35)

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

                                        cTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                <  1:1041566/57‑1 (MQ=11)

                                        cTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                <  1:1302420/57‑1 (MQ=11)

                                          ccAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                >  1:110835/1‑55 (MQ=11)

                                           cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                              >  1:27162/1‑56 (MQ=35)

                                           cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                              >  1:1732113/1‑56 (MQ=35)

                                           cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                              >  1:128415/1‑56 (MQ=35)

                                           cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                              >  1:844509/1‑56 (MQ=35)

                                            aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg               <  1:366324/70‑1 (MQ=255)

                                             gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                          >  1:1273978/1‑42 (MQ=38)

                                               aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                    >  1:1767608/1‑62 (MQ=255)

                                               aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                    >  1:750227/1‑62 (MQ=255)

                                                tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                           <  1:2169006/38‑1 (MQ=40)

                                                    tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc                 <  1:1198644/60‑1 (MQ=255)

                                                    tCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc                 <  1:554817/60‑1 (MQ=255)

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

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

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

                                                        aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGt  >  1:2132578/1‑71 (MQ=255)

                                                        aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGt  <  1:675362/71‑1 (MQ=255)

                                                         gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                          >  1:1218803/1‑46 (MQ=255)

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGT  >  minE/2359818‑2359944

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