BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	43.2%	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	43.2%	104.4 / 157.0	132	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (38/37); new base A (32/25); total (70/62)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGC  >  minE/2359816‑2359935

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

ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcac                                                   <  1:3603681/71‑1 (MQ=255)

 cGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                    >  1:2365016/1‑69 (MQ=31)

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

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

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

    tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                               >  1:911028/1‑71 (MQ=255)

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

       aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                               <  1:2520470/68‑1 (MQ=255)

       aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            >  1:665368/1‑71 (MQ=255)

        tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:1560266/70‑1 (MQ=255)

        tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:2816695/70‑1 (MQ=255)

        tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:954384/70‑1 (MQ=255)

         cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                               >  1:1092255/1‑66 (MQ=255)

          gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                 <  1:3822841/63‑1 (MQ=35)

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

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

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

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

             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                       <  1:2933887/70‑1 (MQ=255)

             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                       <  1:3903482/70‑1 (MQ=255)

             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                        >  1:1792964/1‑69 (MQ=255)

              cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                     <  1:2766417/71‑1 (MQ=34)

              cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                     >  1:2935102/1‑71 (MQ=255)

               ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                    >  1:1142065/1‑71 (MQ=255)

                tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                   >  1:3637354/1‑71 (MQ=255)

                 aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                                >  1:2850442/1‑57 (MQ=35)

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

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

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

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

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

                   gTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTcaca                                                  <  1:2261698/53‑1 (MQ=11)

                    tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                     <  1:853698/65‑1 (MQ=11)

                    tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                               <  1:1801904/71‑1 (MQ=11)

                      aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                             <  1:3128590/71‑1 (MQ=255)

                       gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                    <  1:1265550/47‑1 (MQ=11)

                       gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                    <  1:2001499/47‑1 (MQ=11)

                       gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                    <  1:1885316/47‑1 (MQ=11)

                        cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                           >  1:106496/1‑71 (MQ=255)

                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                       <  1:797804/58‑1 (MQ=35)

                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                  <  1:2373821/63‑1 (MQ=35)

                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                          <  1:3834932/71‑1 (MQ=34)

                         aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            >  1:2143760/1‑53 (MQ=11)

                         aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            >  1:528521/1‑53 (MQ=11)

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

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

                           tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                       >  1:1922644/1‑56 (MQ=11)

                           tCGCACTTCAGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              <  1:1685869/49‑1 (MQ=255)

                             gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                      <  1:1465005/71‑1 (MQ=34)

                             gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                      <  1:3064926/71‑1 (MQ=34)

                             gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                      <  1:2799559/71‑1 (MQ=34)

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

                              cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                     >  1:2169384/1‑71 (MQ=255)

                                cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              <  1:211359/44‑1 (MQ=38)

                                cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            >  1:1440228/1‑46 (MQ=11)

                                cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTAc                                           >  1:1322504/1‑47 (MQ=11)

                                cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                   <  1:2903473/71‑2 (MQ=255)

                                cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                   <  1:1051206/71‑2 (MQ=255)

                                 ttCTGATACCTCCAGCATGCCTCACAGCACACCTTcaca                                                  <  1:808932/39‑1 (MQ=40)

                                  tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                  >  1:333701/1‑70 (MQ=39)

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

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACTGAACGCTcc                                 <  1:2674783/54‑1 (MQ=255)

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

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                             <  1:1193091/58‑1 (MQ=35)

                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                             >  1:1054903/1‑58 (MQ=35)

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

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

                                     gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                       <  1:553395/46‑1 (MQ=38)

                                     gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGc                                    >  1:3553366/1‑49 (MQ=38)

                                     gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                     >  1:1571656/1‑64 (MQ=35)

                                     gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                     >  1:1235587/1‑64 (MQ=35)

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                            >  1:1040399/1‑56 (MQ=35)

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa               <  1:2651780/69‑1 (MQ=40)

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa               <  1:3578356/69‑1 (MQ=40)

                                      aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt             >  1:870980/1‑71 (MQ=255)

                                       tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                       >  1:2184248/1‑44 (MQ=38)

                                       tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg            >  1:3067228/1‑71 (MQ=255)

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

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

                                        aCCTCCAGCATACCTCACAGTACACCTTCACAGGc                                               >  1:438091/1‑35 (MQ=14)

                                        aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                    >  1:2653223/1‑46 (MQ=11)

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

                                          cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg         >  1:2527382/1‑71 (MQ=255)

                                          cTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                       >  1:14631/1‑57 (MQ=11)

                                          cTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                       >  1:1916586/1‑57 (MQ=11)

                                            ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                     >  1:1344074/1‑41 (MQ=38)

                                            ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                    >  1:1582434/1‑58 (MQ=35)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                >  1:1386564/1‑61 (MQ=40)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg            >  1:3036903/1‑65 (MQ=255)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc          >  1:2250758/1‑67 (MQ=255)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc        <  1:427270/69‑1 (MQ=255)

                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg      <  1:79979/71‑1 (MQ=255)

                                              aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg      >  1:2209717/1‑70 (MQ=255)

                                              aGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                               <  1:3605521/45‑1 (MQ=11)

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

                                               gCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                              <  1:3495681/45‑1 (MQ=11)

                                                 aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg   >  1:240438/1‑70 (MQ=255)

                                                  tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                              <  1:2827999/42‑1 (MQ=38)

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

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGC  >  minE/2359816‑2359935

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