BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,863	A→G	6.4%	Q256Q (CAA→CAG)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,863	0	A	G	6.4%	314.9 / 8.2	109	Q256Q (CAA→CAG)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (54/48); new base G (3/4); total (57/52)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 126 total aligned reads displayed.

CTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCGGAAG  >  Exported/2570795‑2570933

cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                                                      <  1:3691402/71‑1 (MQ=255)

cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                       <  1:3172040/70‑1 (MQ=255)

   aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa                                                                     <  1:5462102/69‑1 (MQ=255)

   aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc                                                                    <  1:4349864/70‑1 (MQ=255)

    tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACgg                                                                  >  1:4744885/1‑71 (MQ=255)

      cGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCa                                                                <  1:2631997/71‑1 (MQ=255)

      cGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCa                                                                <  1:2159368/71‑1 (MQ=255)

        cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCa                                                                <  1:1483664/69‑1 (MQ=255)

            ggTACATACGCAGCTGCCTTGCCGCTGCGTACGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGt                                                          >  1:1046851/1‑71 (MQ=255)

             gTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTc                                                         <  1:1125899/71‑1 (MQ=255)

               aCATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGc                                                       >  1:4283257/1‑71 (MQ=255)

                  tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                       <  1:1697807/52‑1 (MQ=255)

                   aCGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCaa                                                   >  1:5118750/1‑71 (MQ=255)

                    cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc                                                                    <  1:4641982/53‑1 (MQ=255)

                    cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                                  >  1:4691553/1‑71 (MQ=255)

                    cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                                  >  1:2123252/1‑71 (MQ=255)

                     gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                                       <  1:4283010/49‑1 (MQ=255)

                      cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                                >  1:1563921/1‑71 (MQ=255)

                      cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                                <  1:2737867/71‑1 (MQ=255)

                      cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                                <  1:2527354/71‑1 (MQ=255)

                       aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTc                                               <  1:5590133/71‑1 (MQ=255)

                       aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTc                                               <  1:447287/71‑1 (MQ=255)

                         cTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTctg                                             <  1:5453895/71‑3 (MQ=255)

                            ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                          <  1:4999601/71‑1 (MQ=255)

                            ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                          <  1:4446137/71‑1 (MQ=255)

                             gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAAc                                                                    >  1:3872822/3‑44 (MQ=12)

                             gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCag                                          <  1:4087924/68‑2 (MQ=255)

                             gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCag                                          <  1:729192/68‑2 (MQ=255)

                               tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCTAAAACGGCAAATTGTCGCCCAAttt                                                >  1:2746410/1‑62 (MQ=255)

                               tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                                >  1:3171260/1‑62 (MQ=255)

                               tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                                >  1:1634008/1‑62 (MQ=255)

                               tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                        >  1:4206615/1‑70 (MQ=255)

                               tGCCGCTGCGTCCGGTGATGGTCTTCAT‑CCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                                >  1:1560527/1‑61 (MQ=255)

                                gCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGt                                                          >  1:2903263/1‑51 (MQ=255)

                                gCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACt                                      >  1:1614448/1‑71 (MQ=255)

                                 ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTc                                     <  1:1531971/71‑1 (MQ=255)

                                 ccGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTc                                     >  1:352963/1‑71 (MQ=255)

                                  cGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTcc                                    <  1:1330329/71‑2 (MQ=255)

                                   gCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAtt                                                            <  1:1895200/46‑1 (MQ=255)

                                   gCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGc                                   <  1:774436/71‑1 (MQ=255)

                                    cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt                                  >  1:3564061/1‑71 (MQ=255)

                                     tacGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCag                                          <  1:1854955/60‑2 (MQ=255)

                                       cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCg                                    <  1:3251621/66‑1 (MQ=255)

                                       cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCg                                    <  1:2480910/66‑1 (MQ=255)

                                         tCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                             <  1:3943/71‑1 (MQ=255)

                                          ccGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                             <  1:5338924/70‑1 (MQ=255)

                                           cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                           >  1:406426/1‑71 (MQ=255)

                                           cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                           <  1:2256062/71‑1 (MQ=255)

                                             gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt                         <  1:2653371/71‑1 (MQ=255)

                                             gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt                         <  1:5055378/71‑1 (MQ=255)

                                              tGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt                         >  1:3074010/1‑70 (MQ=255)

                                              tGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTg                        <  1:5162461/71‑1 (MQ=255)

                                                aTGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                          >  1:155038/1‑51 (MQ=255)

                                                aTGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGa                       <  1:2546995/70‑1 (MQ=255)

                                                 tGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGaaa                     >  1:3445468/1‑71 (MQ=255)

                                                  ggTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAc                    >  1:852979/1‑71 (MQ=255)

                                                  ggTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAc                    >  1:1499124/1‑71 (MQ=255)

                                                   gTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAcc                   <  1:4823359/71‑1 (MQ=255)

                                                    tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt                  >  1:5240482/1‑71 (MQ=255)

                                                     cTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                        >  1:211175/1‑48 (MQ=255)

                                                     cTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTg                 >  1:5380065/1‑71 (MQ=255)

                                                     cTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTg                 >  1:1113757/1‑71 (MQ=255)

                                                      ttCATCCCCGAGCAAAAACGGCAAATTGTCGCCCa                                                    <  1:660700/35‑1 (MQ=255)

                                                      ttCATCCCCGAGCAAAAACGGCAAATTGTCGCCCa                                                    <  1:759311/35‑1 (MQ=255)

                                                      ttCATCCCCGAGCAAAAACGGCAAATTGTCGCCCa                                                    <  1:4878030/35‑1 (MQ=255)

                                                       tCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc                          >  1:4280645/1‑60 (MQ=255)

                                                       tCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc                          >  1:3492808/1‑60 (MQ=255)

                                                        cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                        >  1:260624/1‑45 (MQ=255)

                                                        cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                             <  1:3973375/56‑1 (MQ=255)

                                                        cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                             <  1:359928/56‑1 (MQ=255)

                                                        cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGc               >  1:4302455/1‑70 (MQ=255)

                                                        cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGc               >  1:4865766/1‑70 (MQ=255)

                                                        cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGACTGAAACCTgg                >  1:2235055/1‑69 (MQ=255)

                                                         aTCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCg              <  1:2788367/70‑1 (MQ=255)

                                                          tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt                           >  1:3639720/1‑56 (MQ=255)

                                                          tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat            <  1:2224427/71‑1 (MQ=255)

                                                          tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat            <  1:3762874/71‑1 (MQ=255)

                                                          tCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat            <  1:1789805/71‑1 (MQ=255)

                                                           ccccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat            >  1:5520658/1‑70 (MQ=255)

                                                            cccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAAc                                       >  1:3986705/1‑42 (MQ=255)

                                                            cccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgatg           <  1:3959626/70‑1 (MQ=255)

                                                             ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                        >  1:3021805/1‑40 (MQ=255)

                                                             ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                        >  1:3451386/1‑40 (MQ=255)

                                                             ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                        >  1:5564383/1‑40 (MQ=255)

                                                             ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                                        >  1:662394/1‑40 (MQ=255)

                                                             ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCgat            <  1:2358549/68‑1 (MQ=255)

                                                             ccGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCCCTGTTACGTCTGAAACCTGGCgat            <  1:5389464/68‑1 (MQ=255)

                                                              cGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                                          <  1:4974498/37‑1 (MQ=255)

                                                              cGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGaaa                     >  1:3267301/1‑58 (MQ=255)

                                                                aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTg                                 >  1:2620415/1‑44 (MQ=255)

                                                                aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTg                                 >  1:5443221/1‑44 (MQ=255)

                                                                aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCg      >  1:3725461/1‑71 (MQ=255)

                                                                aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCg      >  1:4182288/1‑71 (MQ=255)

                                                                aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCg      >  1:5579984/1‑71 (MQ=255)

                                                                aGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGATGTTACGTCTGAAACCTGGCGATGATGCg      >  1:635513/1‑71 (MQ=255)

                                                                   aaaaaCGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGatgatg        <  1:1068369/66‑1 (MQ=255)

                                                                    aaaaCGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACg                            <  1:3264233/45‑1 (MQ=255)

                                                                    aaaaCGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGaa                      <  1:5561314/51‑1 (MQ=255)

                                                                    aaaaCGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGc               >  1:3295029/1‑58 (MQ=255)

                                                                    aaaaCGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCGGAAg  >  1:5067464/1‑71 (MQ=255)

CTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGATGATGCGGAAG  >  Exported/2570795‑2570933

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