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	99.4%	721.4 / ‑3.7	178	coding (488/642 nt)	ribC	riboflavin synthase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (87/90); minor base . (1/0); total (88/90)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGGT  >  NC_000913/1742708‑1742801

gTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGa                                              <  2:310128/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTa                                          <  2:4784871/50‑1 (MQ=255)

    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTa                                          <  1:3042793/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

            ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                  >  2:1419296/1‑50 (MQ=255)

             tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                  <  2:5520612/49‑1 (MQ=255)

             tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                  <  1:4669861/49‑1 (MQ=255)

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

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

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

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

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

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

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

                 ttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                              >  1:5010670/1‑49 (MQ=255)

                 ttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                              >  2:2662872/1‑49 (MQ=255)

                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                            <  2:4104862/49‑1 (MQ=255)

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

                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                          <  1:2141721/50‑1 (MQ=255)

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

                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                           <  1:2772181/49‑1 (MQ=255)

                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                           <  2:2118418/49‑1 (MQ=255)

                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                           <  1:4216063/49‑1 (MQ=255)

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

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

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

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

                      aaGAGTCGTGCGGTCCAGTGTTTCCCGAATTAAATGAACGCAAAAAcgcg                        >  2:3333524/1‑50 (MQ=255)

                       agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                       >  1:3942797/1‑50 (MQ=255)

                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                      <  2:1366088/50‑1 (MQ=255)

                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                      <  2:3864318/50‑1 (MQ=255)

                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                      <  1:4429978/50‑1 (MQ=255)

                         agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                     <  1:3227504/50‑1 (MQ=255)

                         agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                     <  2:749534/50‑1 (MQ=255)

                         agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                     <  2:1125940/50‑1 (MQ=255)

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

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

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

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

                             gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                 >  2:4712220/1‑50 (MQ=255)

                             gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                 >  1:1668720/1‑50 (MQ=255)

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

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

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

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

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

                               gCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGa               >  1:3839736/1‑50 (MQ=255)

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

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

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

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

                                   tCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc           >  2:879657/1‑50 (MQ=255)

                                   tCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc           >  1:5005883/1‑50 (MQ=255)

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

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

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

                                     cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGc         <  1:1148140/50‑1 (MQ=255)

                                     cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGc         >  1:2562548/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                          tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTGCGCCGAc     <  2:1678755/49‑1 (MQ=39)

                                            tCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGGt  >  2:1251440/1‑50 (MQ=255)

                                            tCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGGt  >  1:4466508/1‑50 (MQ=255)

GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGGT  >  NC_000913/1742708‑1742801

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