BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,854	C→A	9.7%	I253I (ATC→ATA)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,854	0	C	A	9.7%	321.7 / 12.7	103	I253I (ATC→ATA)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (41/51); new base A (4/6); total (46/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 107 total aligned reads displayed.

TTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGAT  >  Exported/2570784‑2570923

ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATc                                                                       >  1:508147/1‑71 (MQ=255)

  aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATcc                                                                      >  1:415778/1‑70 (MQ=255)

      aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                                                 >  1:2878984/1‑71 (MQ=255)

       ggTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                                                 >  1:3890918/1‑70 (MQ=255)

        gTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCa                                                               <  1:2235105/71‑1 (MQ=255)

          aCTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCag                                                              >  1:5457882/1‑69 (MQ=255)

           cTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCaga                                                             <  1:1211790/70‑3 (MQ=255)

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

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

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

             gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                     >  1:1758223/1‑60 (MQ=255)

             gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa                                                               >  1:1108039/1‑66 (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)

                        gTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                                            <  1:2733484/58‑1 (MQ=255)

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

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

                             tACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGc                                                                <  1:2466572/49‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                     tGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGa                                                                  >  1:3548325/1‑39 (MQ=255)

                                     tGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGa                                                                  >  1:3208350/1‑39 (MQ=255)

                                       ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCgg                                   >  1:5018281/1‑68 (MQ=255)

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

                                       ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCa                                 >  1:364113/1‑70 (MQ=255)

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

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

                                        gctGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCag                                <  1:4087924/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:2480910/66‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                         tGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGc                                  <  1:2582970/51‑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)

                                                            tGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                        <  1:1138752/42‑1 (MQ=255)

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

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

                                                             ggTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAc          >  1:964811/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:1113757/1‑71 (MQ=255)

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

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

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

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

                                                                 ttCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTCCGGCaa                                <  1:733199/45‑1 (MQ=255)

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

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

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

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

                                                                   cATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc                   <  1:3973375/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:1789805/71‑1 (MQ=255)

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

TTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGGCGAT  >  Exported/2570784‑2570923

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