BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,571,214	T→C	12.5%	Y373Y (TAT→TAC)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,571,214	0	T	C	12.5%	282.4 / 25.8	112	Y373Y (TAT→TAC)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (53/45); new base C (4/10); total (57/55)
Fisher's exact test for biased strand distribution p-value = 9.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 114 total aligned reads displayed.

ACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAA  >  Exported/2571144‑2571277

aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTAt                                                                 >  1:682438/1‑71 (MQ=255)

   ccGACCATTTCAGCCATTTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                                                              >  1:1245131/1‑71 (MQ=255)

   ccGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                                                               >  1:5092416/1‑70 (MQ=255)

    cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                                                              <  1:1151396/70‑1 (MQ=255)

      atcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCACTACCttt                                                             <  1:1088495/67‑1 (MQ=255)

       tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtc                                                          <  1:5010097/70‑2 (MQ=255)

       tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtc                                                          <  1:3609129/70‑2 (MQ=255)

       ccATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt                                                          >  1:4091093/1‑71 (MQ=255)

        cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtct                                                         >  1:2350337/1‑69 (MQ=255)

         aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                                        <  1:1718755/71‑1 (MQ=255)

         aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                                        <  1:3148849/71‑1 (MQ=255)

         aTTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTaac                                                          >  1:5161912/1‑67 (MQ=255)

         aTTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtctc                                                        >  1:1381078/1‑68 (MQ=255)

          tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                                        <  1:4417888/70‑1 (MQ=255)

            tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                     >  1:472049/1‑71 (MQ=255)

            tCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGAt                                                     >  1:5009621/1‑71 (MQ=255)

             cAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                     <  1:956006/70‑1 (MQ=255)

              aGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCa                                                   <  1:1321618/71‑1 (MQ=255)

              aGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCa                                                   <  1:2616552/71‑1 (MQ=255)

                ccATGTTTCGGTGATGCGGAAATTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:4346763/71‑1 (MQ=255)

                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt                                                          <  1:3158190/62‑1 (MQ=255)

                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt                                                          <  1:377226/62‑1 (MQ=255)

                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:1570369/71‑1 (MQ=255)

                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:5384325/71‑1 (MQ=255)

                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:1733253/71‑1 (MQ=255)

                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 >  1:4362758/1‑71 (MQ=255)

                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                                                  >  1:3999890/1‑70 (MQ=255)

                ccATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtt                                                 <  1:2614570/71‑1 (MQ=255)

                 cATGTTTCGGTGATTCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                                <  1:5186671/71‑1 (MQ=255)

                 cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                                <  1:5332099/71‑1 (MQ=255)

                 cATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtta                                                <  1:2087464/71‑2 (MQ=255)

                 cATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtta                                                <  1:5103990/71‑2 (MQ=255)

                  aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 >  1:603299/1‑69 (MQ=255)

                  aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 >  1:3875052/1‑69 (MQ=255)

                  aTGTTTCGGTGATGCGGAAAGTCCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                                        >  1:1980416/1‑62 (MQ=255)

                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                                       <  1:353479/62‑1 (MQ=255)

                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                                       <  1:4432622/62‑1 (MQ=255)

                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                                       <  1:4604329/62‑1 (MQ=255)

                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                                              <  1:3632032/71‑1 (MQ=255)

                    gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                                               >  1:2894390/1‑69 (MQ=255)

                      ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                                        >  1:924756/1‑58 (MQ=255)

                      ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                                               >  1:665949/1‑67 (MQ=255)

                       tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                                               >  1:4006483/1‑66 (MQ=255)

                        cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                                                               <  1:5112401/49‑1 (MQ=255)

                        cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                                          <  1:3664760/70‑1 (MQ=255)

                        cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                                          <  1:3658348/70‑1 (MQ=255)

                        cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                                          <  1:5288122/70‑1 (MQ=255)

                        cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                           >  1:4349559/1‑69 (MQ=255)

                         ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                         <  1:42438/70‑1 (MQ=255)

                             ctGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                                                               <  1:5495212/43‑1 (MQ=255)

                             aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                     <  1:3269735/54‑1 (MQ=255)

                             aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                     <  1:4818097/54‑1 (MQ=255)

                             aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                     <  1:1691149/54‑1 (MQ=255)

                             aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:2620025/58‑1 (MQ=255)

                              tGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtta                                                <  1:4467704/58‑2 (MQ=255)

                               gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt                                                         <  1:5496116/48‑1 (MQ=255)

                                 ggAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                                               <  1:4302414/56‑1 (MQ=255)

                                   aaaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt                               >  1:616329/1‑70 (MQ=255)

                                    aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa                             >  1:2923462/1‑71 (MQ=255)

                                    aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa                             >  1:249509/1‑71 (MQ=255)

                                    aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa                             >  1:2522013/1‑71 (MQ=255)

                                    aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa                             <  1:1695461/71‑1 (MQ=255)

                                        gctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                                          >  1:4843058/1‑54 (MQ=255)

                                         ctgctAAGAATGACCAGCTGGATGCATTATCTTTa                                                            >  1:5485272/1‑35 (MQ=255)

                                         ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa                        >  1:2258770/1‑71 (MQ=255)

                                          tgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggagg                                       >  1:2292323/1‑55 (MQ=255)

                                          tgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAc                       >  1:4026501/1‑71 (MQ=255)

                                           gctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                                                   <  1:724465/42‑1 (MQ=255)

                                              aaGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACtt                     >  1:611319/1‑69 (MQ=255)

                                               aGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                     >  1:5504686/1‑36 (MQ=255)

                                                  aTGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTcc               >  1:843344/1‑71 (MQ=255)

                                                  aTGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTcc               >  1:332961/1‑71 (MQ=255)

                                                  aTGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTcc               >  1:5250554/1‑71 (MQ=255)

                                                   tGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:1620436/36‑1 (MQ=255)

                                                   tGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:294048/36‑1 (MQ=255)

                                                   tGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:5580565/36‑1 (MQ=255)

                                                   tGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                 <  1:3713627/36‑1 (MQ=255)

                                                   tGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                                <  1:5005324/37‑1 (MQ=255)

                                                   tGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                                <  1:838098/37‑1 (MQ=255)

                                                   tGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                                <  1:4901180/37‑1 (MQ=255)

                                                    gACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGgaga                                >  1:4602360/1‑52 (MQ=255)

                                                    gACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGgaga                                >  1:2485115/1‑52 (MQ=255)

                                                    gACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGgaga                                >  1:3714040/1‑52 (MQ=255)

                                                     aCCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                         >  1:3385383/1‑42 (MQ=255)

                                                     aCCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                         >  1:3707547/1‑42 (MQ=255)

                                                       cAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                         >  1:1400839/1‑40 (MQ=255)

                                                       cAGCTGGATGCATTATCTTTATTTGGATCATTgag                                              >  1:5615586/1‑35 (MQ=255)

                                                         gCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGc                  >  1:2083339/1‑61 (MQ=255)

                                                         gCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGc                  >  1:615632/1‑61 (MQ=255)

                                                             gATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGAc     >  1:5389856/1‑70 (MQ=255)

                                                               tGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTaa  >  1:809252/1‑71 (MQ=255)

                                                                gCATTATCTTTATTTGGATCATTGagggggaggggg                                    >  1:875681/1‑36 (MQ=255)

                                                                gCATTATCTTTAGTTGGATCATTGagggggaggggg                                    >  1:4927003/1‑36 (MQ=255)

                                                                     aTCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATgg                          >  1:640543/1‑41 (MQ=255)

                                                                     aTCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACtt                     >  1:5452028/1‑46 (MQ=255)

                                                                     aTCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACtt                     >  1:3864963/1‑46 (MQ=255)

                                                                     aTCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTg                   >  1:1434035/1‑48 (MQ=255)

                                                                     aTCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTg                   >  1:3726639/1‑48 (MQ=255)

                                                                     aTCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCt                 >  1:1481853/1‑50 (MQ=255)

                                                                      tCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa                        >  1:5307168/1‑42 (MQ=255)

ACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAA  >  Exported/2571144‑2571277

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