BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,456,456	A→C	5.4%	I298M (ATT→ATG)	yfaQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,456,456	0	A	C	5.4%	327.0 / 13.8	129	I298M (ATT→ATG)	yfaQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (63/59); new base C (5/2); total (68/61)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 139 total aligned reads displayed.

GAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTAC  >  minE/1456387‑1456526

gagaTGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTaa                                                                       >  1:3961112/1‑71 (MQ=255)

gagaTGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTaa                                                                       <  1:3944700/71‑1 (MQ=255)

 agaTGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAt                                                                      <  1:3755709/71‑1 (MQ=255)

  gaTGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAtg                                                                     <  1:1143789/71‑1 (MQ=255)

   aTGCGACTGCAATTTCAGTTTGCCCTTTTCGGTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAtg                                                                     <  1:4032494/70‑1 (MQ=255)

   aTGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAtgt                                                                    <  1:4367491/71‑1 (MQ=255)

   aTGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAtgt                                                                    >  1:3537808/1‑71 (MQ=255)

    tGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAtgtg                                                                   >  1:4536590/1‑71 (MQ=255)

    tGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAtgt                                                                    >  1:4541993/1‑70 (MQ=255)

      cGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGt                                                                 >  1:4193420/1‑71 (MQ=255)

      cGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGt                                                                 >  1:129690/1‑71 (MQ=255)

         cTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTAc                                                              <  1:145055/71‑1 (MQ=255)

            cAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTAcc                                                             <  1:1909374/69‑1 (MQ=255)

            cAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCAt                                                           <  1:931660/71‑1 (MQ=255)

             aaTTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCAtt                                                          <  1:212899/71‑1 (MQ=255)

              aTTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTAc                                                              <  1:2554613/66‑1 (MQ=255)

               tttCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAtgt                                                                    <  1:1742321/59‑1 (MQ=255)

               tttCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTAc                                                              <  1:505525/65‑1 (MQ=255)

               tttCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCa                                                            >  1:1669817/1‑67 (MQ=255)

               tttCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAg                                                        >  1:1116040/1‑71 (MQ=255)

               tttCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAg                                                        >  1:3221115/1‑71 (MQ=255)

                ttCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCa                                                            >  1:1101631/1‑66 (MQ=255)

                 tCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAt                                                                      >  1:2940310/1‑55 (MQ=255)

                   aGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTaaa                                                    <  1:197072/71‑1 (MQ=255)

                   aGTTTGCCCTTTACCGTCAGTAACGTTGTCTCACCGTGGCTGACAAATGTCATGTGGTTACCATTGGCaaa                                                    >  1:1680379/1‑71 (MQ=255)

                   aGTTTGCCCTTTACCGTCAGTAACGTTGTCTCACCGTGGCTGACAAATGTCATGTGGTTACCATTGGCaaa                                                    >  1:2315169/1‑71 (MQ=255)

                     tttGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACg                                                  <  1:13248/71‑1 (MQ=255)

                      ttGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGt                                                 >  1:3505648/1‑71 (MQ=255)

                       tGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAAt                                                                      >  1:3515183/1‑49 (MQ=255)

                       tGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTgg                                                                  >  1:954697/1‑53 (MQ=255)

                       tGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGt                                                 >  1:2697533/1‑70 (MQ=255)

                       tGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGt                                                 <  1:2134998/70‑1 (MQ=255)

                       tGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGt                                                 >  1:2447968/1‑70 (MQ=255)

                       tGCCCTTTACCGTCAGTAACGTTGTCTCACCGTGGCTGACAAATGTCATGTGGTTACCATTGGCAAACGt                                                 >  1:1692362/1‑70 (MQ=255)

                       tGCCCTTTACCGTCAGTAACGTTGTCTCACCGTGGCTGACAAATGTCATGTGGTTACCATTGGCAAACGt                                                 >  1:2971241/1‑70 (MQ=255)

                          ccTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAg                                                        <  1:3021676/60‑1 (MQ=255)

                           cTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCAt                                                           <  1:349392/56‑1 (MQ=255)

                           cTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCacac                                            <  1:1152635/71‑1 (MQ=255)

                           cTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCacac                                            >  1:3888163/1‑71 (MQ=255)

                            ttttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCacac                                            >  1:3697555/1‑70 (MQ=255)

                            ttttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCacac                                            >  1:1601895/1‑70 (MQ=255)

                             tttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCa                                                            >  1:2892318/1‑53 (MQ=255)

                              ttCTCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTa                                                                        >  1:49327/1‑40 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTaa                                                                       <  1:3189253/41‑1 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTaa                                                                       >  1:520937/1‑41 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTaaa                                                    >  1:321178/1‑60 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGt                                                 <  1:3959193/63‑1 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTca                                               >  1:1224652/1‑65 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCacac                                            >  1:18325/1‑68 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGa                                          >  1:4334142/1‑70 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGa                                          <  1:375351/70‑1 (MQ=255)

                              ttCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGa                                          >  1:137050/1‑70 (MQ=255)

                               tCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGa                                          >  1:3259148/1‑69 (MQ=255)

                                cGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAg                                       >  1:4270240/1‑71 (MQ=255)

                                 gCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAgc                                      <  1:3394554/71‑1 (MQ=255)

                                     aGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTaaa                                                    >  1:4490122/1‑53 (MQ=255)

                                     aGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTaaa                                                    >  1:1270617/1‑53 (MQ=255)

                                     aGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGcc                                   <  1:3348072/70‑1 (MQ=255)

                                     aGTAACGTTGTCTCACCGTGGCTGACAAATGTCATGTGGTTACCATTGGCAAACGTCACACGATAGCGccc                                  <  1:1940605/71‑1 (MQ=255)

                                      gTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTaa                                                     <  1:2869499/51‑1 (MQ=255)

                                      gTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGt                                                 <  1:2397824/55‑1 (MQ=255)

                                       tAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAAc                                                   <  1:2532098/52‑1 (MQ=255)

                                        aaGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGt                                                                 >  1:3560791/1‑37 (MQ=255)

                                        aaGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATa                                        >  1:2395052/1‑62 (MQ=255)

                                         aGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGt                                                                 >  1:3926951/1‑36 (MQ=255)

                                         aGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTa                                                         >  1:3256163/1‑44 (MQ=255)

                                         aGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTaaa                                                    >  1:1235974/1‑49 (MQ=255)

                                         aGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCAtt                               <  1:3998669/70‑1 (MQ=255)

                                         aGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCAtt                               <  1:1787811/70‑1 (MQ=255)

                                         aGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCAtt                               <  1:677331/70‑1 (MQ=255)

                                          ggTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTa                                                      <  1:3028805/46‑1 (MQ=255)

                                           gtcgtcTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAAc                                                   >  1:285598/1‑48 (MQ=255)

                                           gtcgtcTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAg                            >  1:1385361/1‑71 (MQ=255)

                                           gtcgtcTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAg                            >  1:4048335/1‑71 (MQ=255)

                                            tcgtcTCACCGTGGCTGACAAACGTAATGTGGTTACCAt                                                           >  1:626329/1‑39 (MQ=255)

                                             cgtcTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTa                          >  1:1801367/1‑71 (MQ=255)

                                              gtcTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTAc                         <  1:569516/71‑1 (MQ=255)

                                               tctcACCGTGTCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCa                                 <  1:750713/62‑1 (MQ=255)

                                               tctcACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTa                              <  1:4150697/65‑1 (MQ=255)

                                                 tcACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAAc                                                   <  1:1480045/42‑1 (MQ=255)

                                                 tcACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTca                                               <  1:698543/46‑1 (MQ=255)

                                                 tcACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGcc                      >  1:964757/1‑71 (MQ=255)

                                                  cACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCAtt                               <  1:1066110/61‑1 (MQ=255)

                                                  cACCGTGCCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGAt                                         >  1:4695925/1‑51 (MQ=255)

                                                   aCCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTca                                               >  1:614383/1‑44 (MQ=255)

                                                   aCCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGAt                                         >  1:483005/1‑50 (MQ=255)

                                                   aCCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGAt                                         >  1:153586/1‑50 (MQ=255)

                                                   aCCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCgg                    >  1:2673482/1‑71 (MQ=255)

                                                    ccGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCgg                    >  1:3444388/1‑70 (MQ=255)

                                                         gCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGAttt               <  1:2443786/70‑1 (MQ=255)

                                                         gCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTa              >  1:963760/1‑71 (MQ=255)

                                                          cTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCgg                    <  1:3922471/64‑1 (MQ=255)

                                                           tGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGAttt               <  1:3843245/68‑1 (MQ=255)

                                                            gACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCAtt                               <  1:862224/51‑1 (MQ=255)

                                                             aCAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACt            <  1:2684853/69‑1 (MQ=255)

                                                               aaaCGTAATGTGGTTACCATTAGTAAACGTCACACGAt                                         >  1:2058863/1‑38 (MQ=255)

                                                                aaCGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGa                  >  1:952922/1‑60 (MQ=255)

                                                                  cGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAg                            <  1:1326156/48‑1 (MQ=255)

                                                                   gTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCcgtcgt    <  1:2058921/71‑1 (MQ=255)

                                                                    tAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTa   <  1:259839/71‑1 (MQ=255)

GAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGTAAACGTCACACGATAGCGCCCATTAAGTACGCCGGGATTTACTGCCGTCGTAC  >  minE/1456387‑1456526

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