BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,289,670	C→T	5.9%	G433G (GGC→GGT)	amn →	AMP nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,289,670	0	C	T	5.9%	316.8 / 5.9	101	G433G (GGC→GGT)	amn	AMP nucleosidase
Reads supporting (aligned to +/- strand): ref base C (45/50); new base T (2/4); total (47/54)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCAAATTGGC  >  minE/1289601‑1289739

cgcgCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGc                                                                       >  1:2303815/1‑70 (MQ=255)

 gcgcAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCga                                                                     <  1:3256473/71‑1 (MQ=255)

     aaGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGAt                                                                  >  1:1579158/1‑70 (MQ=255)

     aaGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGAt                                                                  >  1:2854592/1‑70 (MQ=255)

      aGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                                <  1:2302037/71‑1 (MQ=255)

      aGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                                <  1:2788209/71‑1 (MQ=255)

       ggATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTaa                                                               <  1:1950048/71‑1 (MQ=255)

       ggATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTaa                                                               <  1:984159/71‑1 (MQ=255)

       ggATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTaa                                                               <  1:3739579/71‑1 (MQ=255)

       ggATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTaa                                                               <  1:2759739/71‑1 (MQ=255)

        gATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAACCCGTTGCATGGCGAGATTa                                                                >  1:1050480/1‑69 (MQ=255)

        gATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTaa                                                               >  1:277217/1‑70 (MQ=255)

         atatCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGc                                                                       >  1:2422181/1‑61 (MQ=255)

          tatCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGc                                                                       <  1:2927317/60‑1 (MQ=255)

           atCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                                <  1:4534386/66‑1 (MQ=255)

            tCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACtt                                                           <  1:1352046/70‑1 (MQ=255)

             cGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTcc                                                         >  1:1551728/1‑71 (MQ=255)

             cGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTcc                                                         >  1:2381274/1‑71 (MQ=255)

               tttCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGc                                                                       <  1:4486087/55‑1 (MQ=255)

                   cgcgTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCa                                                   <  1:44926/71‑1 (MQ=255)

                   cgcgTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCa                                                   <  1:23257/71‑1 (MQ=255)

                      gTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGc                                                >  1:1236971/1‑71 (MQ=255)

                         ccATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                                <  1:4078725/52‑1 (MQ=255)

                         ccATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                                <  1:891755/52‑1 (MQ=255)

                         ccATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                                <  1:55404/52‑1 (MQ=255)

                         ccATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTa                                              <  1:646422/70‑1 (MQ=255)

                          cATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGc                                                     <  1:1945255/62‑1 (MQ=255)

                           aTACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAAcc                                           >  1:824775/1‑71 (MQ=255)

                            tACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCg                                          >  1:2279153/1‑71 (MQ=255)

                             aCGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGt                                         >  1:1527425/1‑71 (MQ=255)

                              cGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACtt                                                           <  1:3593614/52‑1 (MQ=255)

                              cGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACtt                                                           <  1:4289592/52‑1 (MQ=255)

                              cGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtt                                        >  1:3242253/1‑71 (MQ=255)

                              cGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtt                                        >  1:1211864/1‑71 (MQ=255)

                              cGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtt                                        >  1:2452555/1‑71 (MQ=255)

                               gggACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCga                                                                     <  1:2691622/41‑1 (MQ=255)

                                ggACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtttt                                      <  1:4655286/71‑1 (MQ=255)

                                ggACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtttt                                      <  1:872064/71‑1 (MQ=255)

                                ggACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtttt                                      <  1:3709662/71‑1 (MQ=255)

                                ggACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGttt                                       <  1:875805/70‑1 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAAc                                            >  1:2365905/1‑62 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGttt                                       >  1:3302726/1‑67 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTAt                                   <  1:1299235/71‑1 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTAt                                   <  1:3516110/71‑1 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTAt                                   <  1:2471279/71‑1 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTAt                                   >  1:3989210/1‑71 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTAt                                   <  1:1926867/71‑1 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTACATGGTGAGATTAAGCTTCCCGGCCAGGCTAACCGTTTTTAt                                   <  1:1267647/71‑1 (MQ=255)

                                   cacTACTGTGTGTTTCAGATAAACCGTTACATGGTGAGATTAAGCTTCCCGGCCAAGCTAACCGTTTTTAt                                   <  1:2373868/71‑1 (MQ=255)

                                    actactGTGTGTTTCTGATAAACCGTTGCATGGCGAGATTAAACt                                                            <  1:4187848/45‑1 (MQ=255)

                                    actactGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTAt                                   >  1:148663/1‑70 (MQ=255)

                                      tactGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtttt                                      <  1:4613717/65‑1 (MQ=255)

                                      tactGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtttt                                      <  1:1705543/65‑1 (MQ=255)

                                      tactGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATg                                  >  1:4626606/1‑69 (MQ=255)

                                      tactGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATg                                  >  1:4551732/1‑69 (MQ=255)

                                      tactGTGTGTTTCAGATAAACCGTTACATGGTGAGATTaagc                                                             >  1:2283517/1‑40 (MQ=37)

                                      tactGTCTGTTTCAGATAAACCGTTACATGGTGAGATTaagc                                                             >  1:2143737/1‑40 (MQ=37)

                                       actGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                                >  1:1772315/1‑38 (MQ=255)

                                       actGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                                >  1:1174857/1‑38 (MQ=255)

                                       actGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGttt                                       >  1:692135/1‑63 (MQ=255)

                                       actGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTAt                                   >  1:2066421/1‑67 (MQ=255)

                                       actGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGaa                                >  1:2174042/1‑70 (MQ=255)

                                       actGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAg                               >  1:3405980/1‑71 (MQ=255)

                                       actGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAg                               >  1:4494379/1‑71 (MQ=255)

                                        ctGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTaa                                                               <  1:1087762/38‑1 (MQ=255)

                                        ctGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTaa                                                               <  1:206467/38‑1 (MQ=255)

                                        ctGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGaa                                <  1:4370242/69‑1 (MQ=255)

                                        ctGTGTGTTTCAGATAAACCGTTACATGGTGAGATTAAGCTTCCCGGCCAGGCTAACCGTTTTTATGAAgg                              <  1:632252/71‑1 (MQ=255)

                                        ctGTGTGTTTCAGATAAACCGTTACATGGTGAGATTAAGCTTCCCGGCCAGGCTAACCGTTTTTATGAAgg                              <  1:278390/71‑1 (MQ=255)

                                         tgtgtgTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTccc                                                        >  1:4070679/1‑44 (MQ=255)

                                         tgtgtgTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTccc                                                        >  1:338132/1‑44 (MQ=255)

                                         tgtgtgTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCg                                          <  1:4693719/58‑1 (MQ=255)

                                           tgtgtTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGcc                                                    >  1:2963138/1‑46 (MQ=255)

                                           tgtgtTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTa                                    >  1:424968/1‑62 (MQ=255)

                                           tgtgtTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGgcg                            <  1:3379990/70‑1 (MQ=255)

                                             tgtTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCTCGGCCAGGCTAACCGTTTTTATGAAGGCGCTa                         >  1:2370320/1‑71 (MQ=255)

                                             tgtTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGttttt                                     <  1:2610626/59‑1 (MQ=255)

                                                 tCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTa                                    >  1:981845/1‑56 (MQ=255)

                                                  cAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAgg                                                 >  1:710171/1‑42 (MQ=255)

                                                  cAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTcc                    >  1:580335/1‑71 (MQ=255)

                                                  cAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTcc                    >  1:2212693/1‑71 (MQ=255)

                                                  cAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTc                     <  1:1359810/70‑1 (MQ=255)

                                                   aGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGtttt                                      >  1:3846399/1‑52 (MQ=255)

                                                   aGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCg                   <  1:716564/71‑1 (MQ=255)

                                                   aGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCg                   <  1:549664/71‑1 (MQ=255)

                                                    gATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCg                   <  1:4028529/70‑1 (MQ=255)

                                                    gATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGa                  <  1:983078/71‑1 (MQ=255)

                                                     aTAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCG‑TTTTATGAAGGCGATATTTCCGaa                 >  1:3180282/1‑70 (MQ=255)

                                                      tAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAAc                <  1:2734341/71‑1 (MQ=255)

                                                          ccGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATc             >  1:3113847/1‑70 (MQ=255)

                                                             ttGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACa               <  1:889487/65‑1 (MQ=255)

                                                             ttGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACa               <  1:295485/65‑1 (MQ=255)

                                                             ttGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCa         <  1:471006/71‑1 (MQ=255)

                                                              tGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACAt              >  1:3744021/1‑65 (MQ=255)

                                                              tGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACAt              >  1:2216951/1‑65 (MQ=255)

                                                              tGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACAt              >  1:2042879/1‑65 (MQ=255)

                                                               gCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCaaa       >  1:1449476/1‑71 (MQ=255)

                                                                   ggCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCa         <  1:1063785/65‑1 (MQ=255)

                                                                    gCGAGATTAAACTTCTCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCAAATTgg   <  1:1521953/70‑1 (MQ=255)

                                                                     cGAGATTAAACTTCCCGGCCAGGCTACCCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCAAATTGGc  >  1:1207490/1‑70 (MQ=255)

CGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCAAATTGGC  >  minE/1289601‑1289739

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