BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,803	C→T	11.4%	R236R (CGC→CGT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,803	0	C	T	11.4%	278.6 / 19.9	105	R236R (CGC→CGT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (44/49); new base T (9/3); total (53/52)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 112 total aligned reads displayed.

AAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCA  >  Exported/2570735‑2570871

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

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

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

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

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

  atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                                                                  >  1:1906756/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

                aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt                                                    <  1:3225794/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:3303062/71‑1 (MQ=255)

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

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

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

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

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

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

                         ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                           <  1:2016074/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:3886713/1‑47 (MQ=255)

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

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

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

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

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

                               cgcCGAGCCATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTatg                                                         >  1:5221741/1‑49 (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)

                                        aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg                            <  1:4950043/71‑1 (MQ=255)

                                        aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg                            >  1:544441/1‑71 (MQ=255)

                                        aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                    <  1:5461260/63‑1 (MQ=255)

                                        aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                    <  1:2249279/63‑1 (MQ=255)

                                         tGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATgg                           >  1:622311/1‑71 (MQ=255)

                                         tGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGa                              <  1:198856/68‑1 (MQ=255)

                                          ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                  >  1:4733665/1‑63 (MQ=255)

                                          ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                  >  1:951436/1‑63 (MQ=255)

                                          ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATgg                           >  1:2486290/1‑70 (MQ=255)

                                          ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt                          >  1:4296935/1‑71 (MQ=255)

                                           gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                         <  1:2328269/71‑1 (MQ=255)

                                           gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                         >  1:3392954/1‑71 (MQ=255)

                                               cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                       <  1:1319461/69‑1 (MQ=255)

                                               cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCa                     >  1:2203160/1‑71 (MQ=255)

                                                gTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                                            <  1:4329799/47‑1 (MQ=255)

                                                 ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                                          <  1:3557122/48‑1 (MQ=255)

                                                 ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                    <  1:2937265/54‑1 (MQ=255)

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

                                                   aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTg                                        >  1:2601556/1‑48 (MQ=255)

                                                   aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTg                                        >  1:2121406/1‑48 (MQ=255)

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

                                                    cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                   <  1:3699059/36‑1 (MQ=255)

                                                    cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                   <  1:879172/36‑1 (MQ=255)

                                                       aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCg                                  >  1:126606/1‑50 (MQ=255)

                                                       aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                                          >  1:4093893/1‑42 (MQ=255)

                                                       aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                         >  1:4223055/1‑59 (MQ=255)

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

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

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

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

                                                         gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTg                               <  1:632709/51‑1 (MQ=255)

                                                         gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                       <  1:4955373/59‑1 (MQ=255)

                                                         gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                       <  1:2622058/59‑1 (MQ=255)

                                                         gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                       <  1:261161/59‑1 (MQ=255)

                                                         gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                       <  1:554302/59‑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)

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

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

                                                              gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa           >  1:1108039/1‑66 (MQ=255)

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

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

                                                                tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCAt                    >  1:2717923/1‑55 (MQ=255)

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

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

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

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

AAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCA  >  Exported/2570735‑2570871

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