BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	484,413	C→T	12.4%	R73R (CGC→CGT)	tolB →	periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	484,413	0	C	T	12.4%	263.1 / 22.7	97	R73R (CGC→CGT)	tolB	periplasmic protein
Reads supporting (aligned to +/- strand): ref base C (51/34); new base T (9/3); total (60/37)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 101 total aligned reads displayed.

CACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAT  >  Exported/484343‑484481

cACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATcgc                                                                      <  1:3783020/71‑1 (MQ=255)

     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCt                                                                   >  1:2156131/1‑69 (MQ=255)

     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:973481/1‑71 (MQ=255)

     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:2346452/1‑71 (MQ=255)

     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:3990787/1‑71 (MQ=255)

     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:1441452/1‑71 (MQ=255)

     gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCg                                                                 >  1:1293722/1‑71 (MQ=255)

      aaGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTc                                                                  <  1:2251711/69‑1 (MQ=255)

      aaGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTc                                                                  <  1:2526762/69‑1 (MQ=255)

      aaGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTc                                                                  <  1:183243/69‑1 (MQ=255)

          tatTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGc                                                            >  1:5032444/1‑71 (MQ=255)

            ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:5024136/71‑1 (MQ=255)

            ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:3934714/71‑1 (MQ=255)

            ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:1911562/71‑1 (MQ=255)

            ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:3968939/71‑1 (MQ=255)

             tGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:415750/70‑1 (MQ=255)

                  gcATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCAcagcag                                                    <  1:2230236/71‑1 (MQ=255)

                   cATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACagcagc                                                   <  1:3470395/71‑1 (MQ=255)

                    aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACagcagc                                                   <  1:4497911/70‑1 (MQ=255)

                    aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  <  1:3116447/71‑1 (MQ=255)

                    aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  <  1:304255/71‑1 (MQ=255)

                        ttGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGt                                              <  1:134081/71‑1 (MQ=255)

                              ctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTgc                                         >  1:3200793/1‑70 (MQ=255)

                              ctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTgc                                         <  1:1512489/70‑1 (MQ=255)

                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:1405030/1‑69 (MQ=255)

                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:2568564/1‑69 (MQ=255)

                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:1248139/1‑69 (MQ=255)

                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:417388/1‑69 (MQ=255)

                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:3322310/1‑69 (MQ=255)

                                    tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:5124175/1‑69 (MQ=255)

                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:3794876/69‑1 (MQ=255)

                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:3697408/69‑1 (MQ=255)

                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:4171059/69‑1 (MQ=255)

                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 >  1:1386852/1‑71 (MQ=255)

                                     gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 >  1:3111809/1‑71 (MQ=255)

                                      cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4077991/57‑1 (MQ=255)

                                      cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4765968/57‑1 (MQ=255)

                                      cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4174140/57‑1 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCAcagca                                                     >  1:2256892/1‑48 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:3124706/1‑57 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:3605645/1‑57 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:1793924/1‑57 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:2316276/1‑57 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:4326772/1‑57 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:188709/1‑57 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:2189797/1‑57 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              <  1:3891415/71‑1 (MQ=255)

                                        tAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              <  1:3475548/71‑1 (MQ=255)

                                         aaCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTg                                                             >  1:2156493/1‑39 (MQ=255)

                                         aaCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              >  1:3781454/1‑70 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:474664/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:5114070/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:767597/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3999399/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:1080455/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:1439391/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:2354194/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:248188/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3026030/1‑43 (MQ=255)

                                          aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3312105/1‑43 (MQ=255)

                                            aGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:4912169/71‑1 (MQ=255)

                                                gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:3290992/67‑1 (MQ=255)

                                                gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:4653797/67‑1 (MQ=255)

                                                gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:5161206/67‑1 (MQ=255)

                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:531551/54‑1 (MQ=255)

                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:4182551/54‑1 (MQ=255)

                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:5189153/54‑1 (MQ=255)

                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACaa                            <  1:4624919/62‑1 (MQ=255)

                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACaa                            <  1:172070/62‑1 (MQ=255)

                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCa                   >  1:217982/1‑71 (MQ=255)

                                                   aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCa                   >  1:572602/1‑71 (MQ=255)

                                                     tttAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAc                           >  1:3086278/1‑61 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  >  1:3363271/1‑35 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:744349/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:627231/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:706034/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1408754/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:76171/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:537950/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1571011/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1582753/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:105250/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1794182/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1921216/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:2102796/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:398656/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:3733995/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:378018/1‑71 (MQ=255)

                                                        aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:3799750/1‑71 (MQ=255)

                                                          tCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCg            >  1:536670/1‑71 (MQ=255)

                                                            cGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGc           <  1:1703316/70‑1 (MQ=255)

                                                              ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 <  1:1191845/46‑1 (MQ=255)

                                                              ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:3202075/71‑1 (MQ=255)

                                                              ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:3689644/71‑1 (MQ=255)

                                                              ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:58573/71‑1 (MQ=255)

                                                                aGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTgg      >  1:3463048/1‑71 (MQ=255)

                                                                aGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTgg      >  1:3944476/1‑71 (MQ=255)

                                                                    cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAt  >  1:2375499/1‑71 (MQ=255)

                                                                    cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAt  >  1:802516/1‑71 (MQ=255)

                                                                    cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGAATGGTCCGCACTGGGCAt  >  1:1718268/1‑71 (MQ=255)

CACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAT  >  Exported/484343‑484481

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