BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,289,664	G→A	5.5%	L431L (TTG→TTA)	amn →	AMP nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,289,664	0	G	A	5.5%	343.4 / 6.5	109	L431L (TTG→TTA)	amn	AMP nucleosidase
Reads supporting (aligned to +/- strand): ref base G (52/51); new base A (2/4); total (54/55)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 117 total aligned reads displayed.

CCATTGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCAAA  >  minE/1289594‑1289734

ccATTGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTg                                                                        >  1:3799643/1‑71 (MQ=255)

ccATTGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTg                                                                        <  1:2663794/71‑1 (MQ=255)

 cATTGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTg                                                                        >  1:3092386/1‑70 (MQ=255)

 cATTGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTg                                                                        >  1:2188400/1‑70 (MQ=255)

  aTTGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTg                                                                        <  1:800973/69‑1 (MQ=255)

  aTTGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGc                                                                       >  1:368945/1‑70 (MQ=255)

   ttGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCAt                                                                     <  1:2561787/71‑1 (MQ=255)

    tGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATg                                                                    >  1:4228816/1‑71 (MQ=255)

     gCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATgg                                                                   >  1:2002578/1‑71 (MQ=255)

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

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

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

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

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

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

              ggATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTaa                                                          <  1:3739579/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)

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

                   tCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGc                                                                       >  1:700326/1‑53 (MQ=255)

                   tCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCa                                                                      >  1:2339518/1‑54 (MQ=255)

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

                   tCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTc                                                     >  1:2343041/1‑71 (MQ=255)

                   tCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTc                                                     >  1:2946910/1‑71 (MQ=255)

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

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

                        tCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATg                                                                    >  1:2103177/1‑51 (MQ=255)

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

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

                           gcgTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACtt                                                      <  1:1685526/62‑1 (MQ=255)

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

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

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

                                ccATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTa                                                           <  1:891755/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)

                                    tcGGGACACTACTGTGTGTTTCAGATAAACCGTTg                                                                        <  1:1659581/34‑1 (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:2452555/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

                                          cacTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTAt                              <  1:2471279/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:1705543/65‑1 (MQ=255)

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

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

                                             tactGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATg                             >  1:4626606/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)

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

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

                                              actGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAg                          >  1:4025378/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:2212693/1‑71 (MQ=255)

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

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

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

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

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

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

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

                                                               aaCCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACAt         <  1:491830/71‑1 (MQ=255)

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

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

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

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

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

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

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

CCATTGCCGCGCAAGGATATCGTTTCCGCGTGCCATACGGGACACTACTGTGTGTTTCAGATAAACCGTTGCATGGCGAGATTAAACTTCCCGGCCAGGCTAACCGTTTTTATGAAGGCGCTATTTCCGAACATCTGCAAA  >  minE/1289594‑1289734

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