BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,742,755	G→C	P163R (CCG→CGG)	ribC ←	riboflavin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,742,755	0	G	C	95.5%	434.0 / 10.1	111	P163R (CCG→CGG)	ribC	riboflavin synthase
Reads supporting (aligned to +/- strand): ref base G (0/5); new base C (52/54); total (52/59)
Fisher's exact test for biased strand distribution p-value = 5.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 111 total aligned reads displayed.

TGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742709‑1742800

tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                            >  2:1735011/1‑50 (MQ=255)

tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                            <  2:2606999/50‑1 (MQ=255)

tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                            >  2:347204/1‑50 (MQ=255)

 gcgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:4786561/1‑50 (MQ=255)

 gcgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:4485348/1‑50 (MQ=255)

 gcgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:213902/1‑50 (MQ=255)

  cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          <  1:4051631/50‑1 (MQ=255)

  cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          <  2:3466590/50‑1 (MQ=255)

  cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  2:4664772/1‑49 (MQ=255)

  cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  2:2399922/1‑49 (MQ=255)

  cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           <  1:5340584/49‑1 (MQ=255)

   gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          >  2:1263477/1‑49 (MQ=255)

   gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          >  1:1968372/1‑49 (MQ=255)

    ccAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaa                                        >  2:793152/1‑50 (MQ=255)

    ccAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaa                                        <  1:3913242/50‑1 (MQ=255)

    ccAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaa                                        <  1:2906706/50‑1 (MQ=255)

     cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       <  1:150846/50‑1 (MQ=255)

     cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       <  2:1244221/50‑1 (MQ=255)

      aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      >  2:255322/1‑50 (MQ=255)

      aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      >  2:2659056/1‑50 (MQ=255)

      aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      >  1:3675197/1‑50 (MQ=255)

      aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      <  2:678891/50‑1 (MQ=255)

       aGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     >  1:1257797/1‑50 (MQ=255)

       aGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     <  2:2903442/50‑1 (MQ=255)

        gTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATgg                                    <  2:2309384/50‑1 (MQ=255)

         ttttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                   <  2:4919103/50‑1 (MQ=255)

         ttttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                   >  1:4861396/1‑50 (MQ=255)

          tttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  >  1:1412403/1‑50 (MQ=255)

           ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGAc                                  <  1:4541940/49‑1 (MQ=255)

           ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  <  2:3757793/49‑1 (MQ=255)

            tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                 >  1:1424450/1‑49 (MQ=255)

             ttcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGc                                >  1:2287741/1‑49 (MQ=255)

             ttcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGc                                <  2:3913143/49‑1 (MQ=255)

              tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                               >  2:2946879/1‑49 (MQ=255)

               cttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaa                              <  1:2099047/49‑1 (MQ=255)

               cttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaa                              >  1:4967779/1‑49 (MQ=255)

                 tcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           >  2:3576490/1‑50 (MQ=255)

                  cccAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCaaaaa                           <  1:708175/49‑1 (MQ=255)

                  cccAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCaaaa                            <  2:3067179/48‑1 (MQ=255)

                  cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaa                            <  2:686821/48‑1 (MQ=255)

                  cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaa                            <  2:2769799/48‑1 (MQ=255)

                  cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaa                            <  1:927812/48‑1 (MQ=255)

                  cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                             <  2:4688610/47‑1 (MQ=255)

                   ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         >  1:1787635/1‑50 (MQ=255)

                    cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        >  2:1387690/1‑50 (MQ=255)

                    cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:2870871/50‑1 (MQ=255)

                    cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  2:461875/50‑1 (MQ=255)

                    cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:866331/50‑1 (MQ=255)

                    cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:3348056/50‑1 (MQ=255)

                     aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                       >  2:3147999/1‑50 (MQ=255)

                     aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                       >  2:558166/1‑50 (MQ=255)

                      agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      <  1:640004/50‑1 (MQ=255)

                      agagTCGTGCGCTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      >  2:2113039/1‑50 (MQ=255)

                       gagTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGtt                     <  2:1705056/50‑1 (MQ=255)

                       gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                     >  1:4636263/1‑50 (MQ=255)

                        agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                    >  2:5379545/1‑50 (MQ=255)

                         gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   >  2:4391417/1‑50 (MQ=255)

                          tCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGc                  <  2:2446157/50‑1 (MQ=255)

                          tCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGc                  >  2:4072002/1‑50 (MQ=255)

                          tCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGc                  <  1:5298120/50‑1 (MQ=255)

                           cGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 >  2:3496763/1‑50 (MQ=255)

                           cGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 >  2:270499/1‑50 (MQ=255)

                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 <  1:3794384/49‑1 (MQ=255)

                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 <  1:4131803/49‑1 (MQ=255)

                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                <  1:3126433/50‑1 (MQ=255)

                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                <  1:3460015/50‑1 (MQ=255)

                             tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:1116162/1‑49 (MQ=255)

                             tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTg               <  2:5244848/50‑1 (MQ=255)

                             tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTg               <  1:4442003/50‑1 (MQ=255)

                               cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGa              <  2:1468196/49‑1 (MQ=255)

                               cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             >  2:2283901/1‑50 (MQ=255)

                               cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             <  1:1392801/50‑1 (MQ=255)

                                gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  2:4425861/1‑50 (MQ=255)

                                gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:2571917/1‑50 (MQ=255)

                                 ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           >  2:2377857/1‑50 (MQ=255)

                                 ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           >  2:2127565/1‑50 (MQ=255)

                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  1:4217062/49‑1 (MQ=255)

                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  2:5227638/49‑1 (MQ=255)

                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  1:63116/49‑1 (MQ=255)

                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  2:3737991/1‑50 (MQ=255)

                                    cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         <  1:4269796/49‑1 (MQ=255)

                                    cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         <  2:1319786/49‑1 (MQ=255)

                                     aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       <  2:4741053/50‑1 (MQ=255)

                                     aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       <  1:1089089/50‑1 (MQ=255)

                                     aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       >  1:482537/1‑50 (MQ=255)

                                     aGTGTTTCCCGAATTAAATGAACGCAAAATCGCGTTGGCGTGACTTCGcc       <  2:1494525/50‑1 (MQ=255)

                                      gtgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       >  2:2200152/1‑49 (MQ=255)

                                       tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     >  2:2320299/1‑50 (MQ=255)

                                       tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     <  2:2070371/50‑1 (MQ=255)

                                       tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     >  1:1812705/1‑50 (MQ=255)

                                       tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     <  2:755346/50‑1 (MQ=255)

                                       tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     <  1:1067119/50‑1 (MQ=255)

                                        gtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    >  1:5075845/1‑50 (MQ=255)

                                        gtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    >  2:902925/1‑50 (MQ=255)

                                         tttCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    <  2:3590076/49‑1 (MQ=255)

                                         tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    <  1:2439417/49‑1 (MQ=255)

                                         tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACg   <  1:3822951/50‑1 (MQ=255)

                                         tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACg   >  1:2183587/1‑50 (MQ=255)

                                         tttCCCGAATGAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    <  1:612130/49‑1 (MQ=255)

                                          ttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACgg  >  2:2228431/1‑50 (MQ=255)

TGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742709‑1742800

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`