BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,774	A→G	20.0%	N227D (AAT→GAT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,774	0	A	G	20.0%	216.4 / 45.3	107	N227D (AAT→GAT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (41/44); new base G (14/7); total (56/51)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 113 total aligned reads displayed.

TGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGG  >  Exported/2570706‑2570840

tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGATCAGACTGTCATTCGCGCGCCGAGCaa                                                                   <  1:4146578/70‑1 (MQ=255)

tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                                   <  1:2134918/70‑1 (MQ=255)

  aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATg                                                                 <  1:5068480/70‑1 (MQ=255)

   gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                               <  1:4654431/71‑1 (MQ=255)

   gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                               <  1:2380771/71‑1 (MQ=255)

   gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                               <  1:4609709/71‑1 (MQ=255)

   gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                               <  1:1826526/71‑1 (MQ=255)

   gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                               <  1:4739520/71‑1 (MQ=255)

   gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                               <  1:249362/71‑1 (MQ=255)

     ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                                  >  1:750550/1‑66 (MQ=255)

     ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTa                                                             >  1:2094604/1‑71 (MQ=255)

       aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                           >  1:2652771/1‑71 (MQ=255)

       aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                           >  1:3452316/1‑71 (MQ=255)

       aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                           >  1:5260611/1‑71 (MQ=255)

       aaGAGCGCAACTTACTGAAGCAAAATAAAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                           >  1:4828358/1‑71 (MQ=255)

        agagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGt                                                          <  1:5429040/71‑1 (MQ=255)

          agCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTa                                                        >  1:4617417/1‑71 (MQ=255)

           gcgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTAc                                                       >  1:217300/1‑71 (MQ=255)

              cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                                      >  1:2620716/1‑69 (MQ=255)

               aaCTTACTGAAGCAAAATATAACCTTGAGTAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                              <  1:1854577/60‑1 (MQ=255)

               aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                              <  1:2430738/60‑1 (MQ=255)

               aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                              <  1:5539252/60‑1 (MQ=255)

                aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCa                                                                    <  1:2932552/53‑1 (MQ=255)

                aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgt                                                  <  1:1579372/71‑2 (MQ=255)

                aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                  >  1:766332/1‑71 (MQ=255)

                aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                  >  1:5504789/1‑71 (MQ=255)

                 cTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                                   >  1:2418141/1‑53 (MQ=255)

                    aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                                   >  1:2963613/1‑50 (MQ=255)

                    aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCATGTa                                                <  1:3341062/69‑1 (MQ=255)

                    aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                              >  1:129133/1‑71 (MQ=255)

                     cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                                                <  1:1242373/68‑1 (MQ=255)

                      tGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATg                                                                 >  1:3481438/1‑50 (MQ=255)

                       gAAGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATggtt                                                              <  1:311883/52‑3 (MQ=255)

                        aaGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                                          >  1:5325897/1‑71 (MQ=255)

                        aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                                          >  1:5628654/1‑70 (MQ=255)

                        aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                                          >  1:2660324/1‑70 (MQ=255)

                           cAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                               <  1:5421041/47‑1 (MQ=255)

                           cAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCg                                        <  1:930966/70‑1 (MQ=255)

                            aaaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGAc                                                       >  1:4144447/1‑54 (MQ=255)

                            aaaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgt                                       <  1:212453/70‑2 (MQ=255)

                             aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGtt                                                         <  1:3532462/51‑1 (MQ=255)

                             aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGAt                                                                  <  1:1039785/42‑1 (MQ=12)

                             aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGAt                                                                  <  1:2680923/42‑1 (MQ=12)

                             aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                                     >  1:1253642/1‑71 (MQ=255)

                             aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                                     <  1:3208415/71‑1 (MQ=255)

                              aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCa                                    <  1:1698509/71‑1 (MQ=255)

                              aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTa                                                >  1:3456040/1‑59 (MQ=255)

                              aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                                    >  1:2880584/1‑70 (MQ=255)

                              aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                                    >  1:2733471/1‑70 (MQ=255)

                              aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                                    >  1:2610651/1‑70 (MQ=255)

                                tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGt                                                          <  1:2890290/47‑1 (MQ=255)

                                tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAgg                                  >  1:3649775/1‑71 (MQ=255)

                                  taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                >  1:4119278/1‑71 (MQ=255)

                                  taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                >  1:1397815/1‑71 (MQ=255)

                                  taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                <  1:3749721/71‑1 (MQ=255)

                                    aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                              <  1:2085988/71‑1 (MQ=255)

                                     ccTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgtc                                      >  1:3467639/1‑60 (MQ=255)

                                      cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                  <  1:449845/49‑1 (MQ=255)

                                      cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                  <  1:2222215/49‑1 (MQ=255)

                                      cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                                  <  1:2274309/49‑1 (MQ=255)

                                      cTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTa                            >  1:1271195/1‑71 (MQ=255)

                                       ttGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGa                                            >  1:562770/1‑54 (MQ=255)

                                        tGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACg                          >  1:2174325/1‑71 (MQ=255)

                                         gAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                            >  1:2387129/1‑52 (MQ=255)

                                          aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc                         <  1:2295133/70‑1 (MQ=255)

                                          aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc                         <  1:1028155/70‑1 (MQ=255)

                                             aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt                     <  1:3225794/71‑1 (MQ=255)

                                             aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt                     <  1:3218372/71‑1 (MQ=255)

                                              gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGc                                       >  1:359309/1‑52 (MQ=255)

                                              gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                              <  1:1560887/61‑1 (MQ=255)

                                              gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                    <  1:2670348/71‑1 (MQ=255)

                                              gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                    >  1:3098634/1‑71 (MQ=255)

                                                cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                      >  1:5506930/1‑67 (MQ=255)

                                                cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                      >  1:917187/1‑67 (MQ=255)

                                                cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                      >  1:3899776/1‑67 (MQ=255)

                                                   tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCtt                <  1:124460/70‑1 (MQ=255)

                                                     aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc             >  1:2726988/1‑71 (MQ=255)

                                                      ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg            <  1:2016074/71‑1 (MQ=255)

                                                      ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg            <  1:4652072/71‑1 (MQ=255)

                                                      ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg            <  1:96965/71‑1 (MQ=255)

                                                      ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg            <  1:1277911/71‑1 (MQ=255)

                                                      ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg            <  1:3065483/71‑1 (MQ=255)

                                                      ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg            <  1:5155473/71‑1 (MQ=255)

                                                      ttCGCGCGCCGAGCAAAGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg            <  1:1593534/71‑1 (MQ=255)

                                                       tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc           >  1:4405786/1‑71 (MQ=255)

                                                        cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg            >  1:4819656/1‑69 (MQ=255)

                                                           gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                               >  1:974992/1‑47 (MQ=255)

                                                           gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                               >  1:3886713/1‑47 (MQ=255)

                                                           gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                               >  1:3157193/1‑47 (MQ=255)

                                                           gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                               >  1:2616414/1‑47 (MQ=255)

                                                           gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                               >  1:2029958/1‑47 (MQ=255)

                                                            cgcCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGt      >  1:4854762/1‑71 (MQ=255)

                                                             gcagAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc     >  1:5095792/4‑71 (MQ=255)

                                                             gcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc     >  1:5027913/1‑71 (MQ=255)

                                                                gAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg  <  1:4523482/71‑1 (MQ=255)

                                                                gAGCAATGGCTACGTTACTCAGGTACTGAACCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTcc    >  1:2462052/1‑69 (MQ=255)

                                                                 aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                    >  1:3905850/1‑52 (MQ=255)

                                                                 aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                    >  1:1658442/1‑52 (MQ=255)

                                                                  gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg  <  1:3771478/69‑1 (MQ=255)

                                                                  gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg  <  1:2791965/69‑1 (MQ=255)

TGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGG  >  Exported/2570706‑2570840

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