BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	42.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	42.7%	64.9 / 92.4	75	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (32/11); new base A (20/12); total (52/23)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTG  >  minE/2359818‑2359931

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:178057/1‑47 (MQ=38)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                           >  1:139268/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:1271579/68‑1 (MQ=255)

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

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

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

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

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

                              cTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                >  1:268285/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:68344/1‑35 (MQ=33)

                                 cTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                >  1:1906932/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:128415/1‑56 (MQ=35)

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

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

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

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

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

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

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

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTG  >  minE/2359818‑2359931

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