BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	214,885	C→T	T18I (ACC→ATC)	arfB →	peptidyl‑tRNA hydrolase, ribosome rescue factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	214,885	0	C	T	100.0%	187.7 / NA	47	T18I (ACC→ATC)	arfB	peptidyl‑tRNA hydrolase, ribosome rescue factor
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (24/23); total (24/23)

TGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCACCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGTTAATAAGAC  >  NC_000913/214838‑214930

tgtgATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCg                                             >  1:3644116/1‑50 (MQ=255)

 gcgATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGc                                            >  1:3255260/3‑50 (MQ=255)

  tgATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGcc                                           <  1:3000420/50‑1 (MQ=255)

  tgATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGc                                            >  1:1559392/1‑49 (MQ=255)

  tgATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGc                                            <  1:3322427/49‑1 (MQ=255)

  tgATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGc                                            >  1:3113792/1‑49 (MQ=255)

     tttCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCAtt                                        >  1:852192/1‑50 (MQ=255)

     tttCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCAtt                                        >  1:526521/1‑50 (MQ=255)

       tCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTc                                       <  1:3051366/49‑1 (MQ=255)

       tCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTc                                       <  1:2616843/49‑1 (MQ=255)

       tCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCg                                      >  1:1634153/1‑50 (MQ=255)

        cccGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCg                                      <  1:2392946/49‑1 (MQ=255)

         ccGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGt                                     <  1:3309116/49‑1 (MQ=255)

         ccGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTg                                    >  1:235344/1‑50 (MQ=255)

          cGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTgc                                   >  1:1499578/1‑50 (MQ=255)

          cGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTgc                                   >  1:3285886/1‑50 (MQ=255)

            aCATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTgcg                                  >  1:3798086/1‑49 (MQ=255)

            aCATGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTgcg                                  >  1:1990328/1‑49 (MQ=255)

             cATGTTGCTATTCCCGATGGTGAGTTTGAGATCATCGCCATTCGTgcg                                  >  1:2859480/1‑48 (MQ=38)

              aTGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAg                               >  1:3541736/1‑50 (MQ=255)

              aTGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAg                               >  1:1358660/1‑50 (MQ=255)

               tGTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAgg                              <  1:1527012/50‑1 (MQ=255)

                gTTGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAggg                             >  1:2250363/1‑50 (MQ=255)

                 ttGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGgc                            >  1:2915212/1‑50 (MQ=255)

                 ttGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGgc                            <  1:3964293/50‑1 (MQ=255)

                 ttGCTATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGgc                            <  1:3201858/50‑1 (MQ=255)

                     tATTCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCgcgg                        >  1:371755/1‑50 (MQ=255)

                       ttCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCgcgggc                      <  1:2751747/50‑1 (MQ=255)

                       ttCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCgcgggc                      <  1:1849380/50‑1 (MQ=255)

                        tCCCGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCgcgggcg                     >  1:3707208/1‑50 (MQ=255)

                         cccGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCgcgggcgg                    >  1:488522/1‑50 (MQ=255)

                         cccGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCgcgggcgg                    >  1:2419116/1‑50 (MQ=255)

                          ccGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCgcgggcgg                    <  1:234690/49‑1 (MQ=255)

                          ccGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGcgggcggg                   <  1:3041784/50‑1 (MQ=255)

                           cGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGcgggcggg                   >  1:1495284/1‑49 (MQ=255)

                           cGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGCgggcgggc                  <  1:3813716/50‑1 (MQ=255)

                           cGATGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGCgggcgggc                  <  1:1818236/50‑1 (MQ=255)

                              tGGTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGCGGGCGGgcagc               <  1:3617326/50‑1 (MQ=255)

                                gTGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCAt             <  1:3143366/50‑1 (MQ=255)

                                 tGAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATg            <  1:813960/50‑1 (MQ=255)

                                  gAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGt           >  1:3852120/1‑50 (MQ=255)

                                  gAGCTTGAGATCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGt           >  1:2528176/1‑50 (MQ=255)

                                     cTTGAGATCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGTtaa        <  1:996796/50‑1 (MQ=255)

                                         agaTCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGTTAATAAg    <  1:3723068/50‑1 (MQ=255)

                                         agaTCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGTTAATAAg    <  1:3074942/50‑1 (MQ=255)

                                         agaTCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGTTAATAAg    <  1:1890812/50‑1 (MQ=255)

                                          gaTCATCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGTTAATAAGa   <  1:3234686/50‑1 (MQ=255)

                                           atcatcGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGTTAATAAGAc  <  1:3753325/50‑1 (MQ=255)

TGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCACCGCCATTCGTGCGCAGGGCGCGGGCGGGCAGCATGTTAATAAGAC  >  NC_000913/214838‑214930

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`