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	100.0%	174.1 / NA	44	P163R (CCG→CGG)	ribC	riboflavin synthase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (17/27); total (17/27)

GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742708‑1742800

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                            tCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACgg  <  1:2108965/49‑1 (MQ=255)

GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742708‑1742800

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