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.3%	531.7 / ‑4.1	135	P163R (CCG→CGG)	ribC	riboflavin synthase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (68/66); minor base T (0/1); total (68/67)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 135 total aligned reads displayed.

GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742708‑1742800

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

 tGCTCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                            <  2:5289722/50‑1 (MQ=255)

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

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

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

    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtga                                         <  2:4207718/50‑3 (MQ=255)

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

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

    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGCTCCAGTGTTTCCCGAAtt                                          >  1:1681899/1‑49 (MQ=255)

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

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

      cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       >  1:5338352/1‑50 (MQ=255)

      cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       >  1:4308229/1‑50 (MQ=255)

      cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       >  2:4189058/1‑50 (MQ=255)

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

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

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

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

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

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

         gTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     >  2:3862360/1‑49 (MQ=255)

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

           tttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGTAc                                  >  2:1427395/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

               tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                               <  2:5631514/49‑1 (MQ=255)

               tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                               <  1:5397482/49‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          >  2:3017473/1‑50 (MQ=255)

                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          >  2:2104653/1‑50 (MQ=255)

                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          >  2:1693987/1‑50 (MQ=255)

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

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

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

                     cAAGAGTCGTGCGTTCCGGTGTTTCCCGAATTAAGTGGACGCAAAAAcgc                        >  2:2221510/1‑50 (MQ=255)

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

                      aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                       <  1:2805173/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

                          gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   <  2:3006037/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                 gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            <  1:5015670/50‑1 (MQ=255)

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

                                  ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           <  2:4728490/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                       gtgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCg      <  2:2211496/50‑1 (MQ=255)

                                       gtgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCg      <  1:1240052/50‑1 (MQ=255)

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

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

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

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

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

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

GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742708‑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: `‑`