BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,180,200	T→G	V8G (GTG→GGG)	yjaH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,180,200	0	T	G	98.0%	201.5 / ‑3.9	49	V8G (GTG→GGG)	yjaH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (27/21); total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAG  >  REL606/4180129‑4180272

tGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGtt                                                                      >  1:1059140/1‑76 (MQ=255)

tGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGtt                                                                      <  2:1185140/76‑1 (MQ=255)

 gTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGttt                                                                     >  1:552281/1‑76 (MQ=255)

  tttCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTc                                                                    >  1:1322110/1‑76 (MQ=255)

  tttCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTc                                                                    >  1:1034613/1‑76 (MQ=255)

   ttCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCa                                                                   <  2:1227573/76‑1 (MQ=255)

   ttCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCa                                                                   >  1:147322/1‑76 (MQ=255)

     cTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCAcc                                                                 <  1:47236/76‑1 (MQ=255)

         cAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCtt                                                             >  1:288201/1‑76 (MQ=255)

           aGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCtttt                                                           <  2:1514384/76‑1 (MQ=255)

           aGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCtttt                                                           <  1:1488223/76‑1 (MQ=255)

            ggCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTg                                                          <  2:975706/76‑1 (MQ=255)

            ggCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTg                                                          <  2:1260197/76‑1 (MQ=255)

             gCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGt                                                         <  1:1223072/76‑1 (MQ=255)

             gCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGt                                                         >  1:638933/1‑76 (MQ=255)

               aCTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCt                                                       >  2:1316791/1‑76 (MQ=255)

               aCTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCt                                                       >  1:940647/1‑76 (MQ=255)

                cTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCtt                                                      >  1:583027/1‑76 (MQ=255)

                 tGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCttt                                                     >  1:511550/1‑76 (MQ=255)

                 tGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCttt                                                     >  2:967651/1‑76 (MQ=255)

                    aaGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTg                                                  >  1:235857/1‑76 (MQ=255)

                    aaGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTg                                                  <  2:1541887/76‑1 (MQ=255)

                    aaGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTg                                                  >  2:1247154/1‑76 (MQ=255)

                      gACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGc                                                >  2:1190517/1‑76 (MQ=255)

                            gTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGAt                                          >  1:1237412/1‑76 (MQ=255)

                              taagtaagATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCa                                        <  1:905175/76‑1 (MQ=255)

                               aagtaagATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAt                                       >  2:538233/1‑76 (MQ=255)

                                 gtaagATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTg                                     <  1:106258/76‑1 (MQ=255)

                                  taagATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGa                                    >  2:1601660/1‑76 (MQ=255)

                                   aagATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGAt                                   <  2:983460/76‑1 (MQ=255)

                                         gCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGc                             >  2:556234/1‑76 (MQ=255)

                                         gCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGc                             >  1:698919/1‑76 (MQ=255)

                                          cGTGGCAGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCa                            <  1:1294849/76‑1 (MQ=255)

                                             ggCAGTGAACAGTTTTAACGAAGGGGGGGGTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAggg                         >  1:899065/1‑75 (MQ=255)

                                                aGTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcg                      >  1:1099013/1‑76 (MQ=255)

                                                 gTGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcgc                     <  2:1369199/76‑1 (MQ=255)

                                                  tGAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCt                    >  1:844103/1‑76 (MQ=255)

                                                   gAACAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCGGGTGCGCTg                   <  2:1497119/76‑1 (MQ=255)

                                                      cAGTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGttt                <  2:147322/76‑1 (MQ=255)

                                                       aGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtc               >  2:1550894/1‑76 (MQ=255)

                                                        gTTTTAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtct              <  2:1562077/76‑1 (MQ=255)

                                                          tttAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCa            <  2:1001444/76‑1 (MQ=255)

                                                           ttAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCAc           >  1:1222200/1‑76 (MQ=255)

                                                            tAACGAAGGGGGGGTTTCACCCCTTGTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACt          <  1:1486013/76‑1 (MQ=255)

                                                             aaCGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACTg         >  2:220090/1‑76 (MQ=255)

                                                               cGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGcc       >  1:243657/1‑76 (MQ=255)

                                                                gAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCt      <  2:794609/76‑1 (MQ=255)

                                                                 aaGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTg     <  1:1200611/76‑1 (MQ=255)

                                                                 aaGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTg     >  2:587752/1‑76 (MQ=255)

                                                                    ggggggTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg  <  2:1435134/76‑1 (MQ=255)

TGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAG  >  REL606/4180129‑4180272

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`