BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	455,040	A→C	31.3%	V15V (GTT→GTG)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	455,040	0	A	C	31.3%	95.0 / 105.7	112	V15V (GTT→GTG)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (37/40); new base C (12/23); total (49/63)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 115 total aligned reads displayed.

TTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAA  >  Exported/454972‑455107

ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                                                                   >  1:2081188/1‑71 (MQ=255)

tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                                                                   <  1:3353582/69‑1 (MQ=255)

tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                                                                   <  1:3263095/69‑1 (MQ=255)

tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                                                                   <  1:3862299/69‑1 (MQ=255)

      agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                                              <  1:3771237/70‑1 (MQ=255)

      agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                                              <  1:3291639/70‑1 (MQ=255)

       gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg                                                            >  1:3082007/1‑71 (MQ=255)

       gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg                                                            <  1:265496/71‑1 (MQ=255)

        aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                             <  1:2225637/69‑1 (MQ=255)

        aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                             >  1:3039509/1‑69 (MQ=255)

         tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                             <  1:5071689/68‑1 (MQ=255)

         tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt                                                           <  1:4606702/70‑1 (MQ=255)

         tAGAGAATCAGCAGTGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                          <  1:1399617/71‑1 (MQ=255)

         tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCTCAACGACATAGTa                                                          <  1:1553546/71‑1 (MQ=255)

         tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                          <  1:5403146/71‑1 (MQ=255)

         tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTACCAGATTTGCAGCGCAACGACATAGTa                                                          <  1:4030605/71‑1 (MQ=255)

          agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                                          >  1:136720/1‑70 (MQ=255)

          agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                          <  1:4914839/70‑1 (MQ=255)

           gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                                                               >  1:3129678/1‑62 (MQ=255)

           gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                                                               >  1:2970162/1‑62 (MQ=255)

           gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt                                                           >  1:4002479/1‑68 (MQ=255)

           gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTata                                                        <  1:700401/71‑3 (MQ=255)

           gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTat                                                         <  1:700044/70‑2 (MQ=255)

           gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTat                                                         <  1:5468723/70‑2 (MQ=255)

           gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        >  1:328063/1‑71 (MQ=255)

           gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        <  1:3796277/71‑1 (MQ=255)

           gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        <  1:762769/71‑1 (MQ=255)

           gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        <  1:4334378/71‑1 (MQ=255)

           gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        <  1:2283969/71‑1 (MQ=255)

           gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        <  1:3352644/71‑1 (MQ=255)

           gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        <  1:2087903/71‑1 (MQ=255)

           gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        <  1:4845549/71‑1 (MQ=255)

            agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatat                                                       <  1:386759/71‑1 (MQ=255)

              aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                                             >  1:4607506/1‑63 (MQ=255)

              aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                                             >  1:1265823/1‑63 (MQ=255)

               aTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                                     <  1:1708427/70‑1 (MQ=255)

                 cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                                                  >  1:4318729/1‑71 (MQ=255)

                 cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                                                  >  1:2903474/1‑71 (MQ=255)

                  agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                                 >  1:626790/1‑71 (MQ=255)

                  agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                                 >  1:1768082/1‑71 (MQ=255)

                   gcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        >  1:494033/1‑63 (MQ=255)

                    taaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATcgcgcgcg                                               >  1:4027102/4‑71 (MQ=255)

                    cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                                               >  1:925877/1‑71 (MQ=255)

                     agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                                               <  1:283032/54‑1 (MQ=255)

                     agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                              >  1:1929373/1‑71 (MQ=255)

                     aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                                              >  1:1049272/3‑71 (MQ=255)

                     aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                                              >  1:7948/3‑71 (MQ=255)

                         cACGCATGCCGTAAAAGCTAAAATATTTCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                          <  1:1217570/71‑1 (MQ=255)

                         cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                          <  1:4509645/71‑1 (MQ=255)

                         cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                          <  1:4503315/71‑1 (MQ=255)

                         cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                          <  1:2786885/71‑1 (MQ=255)

                         cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                           <  1:1933962/70‑1 (MQ=255)

                         cACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                          <  1:3778408/71‑1 (MQ=255)

                           cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                                        <  1:294826/71‑1 (MQ=255)

                            gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                       >  1:2119683/1‑70 (MQ=255)

                             cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                                                                   >  1:3208739/1‑42 (MQ=255)

                             cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGTTGTGATg                                      >  1:3128724/1‑71 (MQ=255)

                             cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAg                                                            <  1:3572438/49‑1 (MQ=255)

                             cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                      <  1:121246/71‑3 (MQ=255)

                               tGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                          <  1:4810159/49‑1 (MQ=255)

                               tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatat                                                       <  1:4975880/52‑1 (MQ=255)

                               tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                      <  1:1998365/69‑3 (MQ=255)

                                 ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                                                     >  1:4173295/1‑36 (MQ=255)

                                 ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                                                     >  1:4221607/1‑36 (MQ=255)

                                 ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                                                     >  1:4484290/1‑36 (MQ=255)

                                    tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                                             <  1:4168569/41‑1 (MQ=255)

                                    tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                                           <  1:941405/43‑1 (MQ=255)

                                    tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                                           <  1:1091625/43‑1 (MQ=255)

                                    tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                        <  1:4982778/46‑1 (MQ=255)

                                    tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                          <  1:5459242/60‑1 (MQ=255)

                                      aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAg                                                            <  1:236756/40‑1 (MQ=255)

                                      aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                             <  1:3652356/70‑3 (MQ=255)

                                      aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                             <  1:1798406/70‑3 (MQ=255)

                                      aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                             <  1:4056179/70‑3 (MQ=255)

                                      aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt                                                              >  1:776803/1‑38 (MQ=255)

                                      aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                                     >  1:1839905/1‑47 (MQ=255)

                                         gCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGc                                             <  1:2532398/52‑1 (MQ=255)

                                           tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                          <  1:515501/53‑1 (MQ=255)

                                            aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                              >  1:4483838/1‑64 (MQ=255)

                                             aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat                      <  1:4807160/71‑1 (MQ=255)

                                                tatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt                                <  1:1532654/58‑1 (MQ=255)

                                                     ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                >  1:4394748/1‑69 (MQ=255)

                                                     ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac              <  1:5599699/71‑1 (MQ=255)

                                                     ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca               >  1:3020669/1‑70 (MQ=255)

                                                     ccAGATTTGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca               >  1:1175136/1‑70 (MQ=255)

                                                      cAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca             <  1:1510223/71‑1 (MQ=255)

                                                       aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                        >  1:2925206/1‑43 (MQ=255)

                                                       aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                        >  1:3796198/1‑43 (MQ=255)

                                                        gATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAt           >  1:5553338/1‑71 (MQ=255)

                                                         aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt          <  1:5153035/71‑1 (MQ=255)

                                                         aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt          >  1:1924153/1‑71 (MQ=255)

                                                         aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt          <  1:3458865/71‑1 (MQ=255)

                                                           ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa        >  1:5031569/1‑71 (MQ=255)

                                                           ttGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa        >  1:737611/1‑71 (MQ=255)

                                                             gCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa        <  1:5020296/69‑1 (MQ=255)

                                                              cAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACgg     >  1:3548632/1‑71 (MQ=255)

                                                               aGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc    <  1:1001850/71‑1 (MQ=255)

                                                                gcgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc    <  1:4602369/70‑1 (MQ=255)

                                                                 cgcCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                             >  1:2435057/1‑41 (MQ=12)

                                                                 cgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCaa  >  1:4018329/1‑71 (MQ=255)

TTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAA  >  Exported/454972‑455107

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