BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,958,877	2 bp→CG	coding (586‑587/2022 nt)	rep →	rep helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,958,877	0	G	C	100.0%	272.9 / NA	102	A196R (GCC→CGC)	rep	rep helicase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (52/50); total (52/50)
*	NC_000913	3,958,878	0	C	G	100.0%	270.2 / NA	102	A196R (GCC→CGC)	rep	rep helicase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (52/50); total (52/50)

Only 100 of 102 total aligned reads displayed.

CACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGAC  >  NC_000913/3958809‑3958945

||

cccacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATga                                                               <  2:1817156/74‑1 (MQ=255)

 acacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                              <  1:3568610/76‑1 (MQ=255)

  cacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                             <  1:3341749/76‑1 (MQ=255)

   acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                            >  1:408518/1‑76 (MQ=255)

    ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                           <  1:1580346/76‑1 (MQ=255)

    ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                           <  1:3513156/76‑1 (MQ=255)

    ccTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                           <  2:1417872/76‑1 (MQ=255)

     cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGaagaag                                                          >  2:1281562/1‑76 (MQ=255)

     cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGaagaag                                                          <  2:2786375/76‑1 (MQ=255)

      tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGt                                                         <  2:1131345/76‑1 (MQ=255)

       gAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                        <  2:1261201/76‑1 (MQ=255)

        aaaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTcc                                                       >  1:3465705/1‑76 (MQ=255)

        aaaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTcc                                                       <  1:2618605/76‑1 (MQ=255)

        aaaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTcc                                                       <  1:2513566/76‑1 (MQ=255)

        aaaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTcc                                                       <  2:2118225/76‑1 (MQ=255)

         aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                      >  2:3412900/1‑76 (MQ=255)

          tgCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGc                                                     >  2:4726783/2‑76 (MQ=255)

          aGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGCAGTCCGc                                                     >  1:3220406/1‑76 (MQ=255)

           gCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                    >  1:1947135/1‑76 (MQ=255)

           gCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                    <  1:3919963/76‑1 (MQ=255)

           accTGTAGCGTTCACGTCTTGGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                    <  2:2132761/75‑1 (MQ=255)

              tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                 <  1:4236387/76‑1 (MQ=255)

              tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                 <  2:1818091/76‑1 (MQ=255)

              tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                 <  2:4003485/76‑1 (MQ=255)

              tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                 <  1:3202813/76‑1 (MQ=255)

              tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                 >  1:2692359/1‑76 (MQ=255)

               gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcg                                                <  2:3397628/76‑1 (MQ=255)

               gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcg                                                >  1:3609488/1‑76 (MQ=255)

                tAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                               >  2:3706798/1‑76 (MQ=255)

                 atcGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCt                                              <  1:4252013/74‑1 (MQ=255)

                 aaagTTCTAGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCCATGAAGAAGTGCGCAAGCGCt                                              >  2:604771/4‑76 (MQ=255)

                  aCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                             <  2:2838292/76‑1 (MQ=255)

                    gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGc                                           >  1:3056120/1‑76 (MQ=255)

                    gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGc                                           >  1:3849407/1‑76 (MQ=255)

                     ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                          <  2:4435428/76‑1 (MQ=255)

                     ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                          >  2:3055544/1‑76 (MQ=255)

                     ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                          <  1:4414020/76‑1 (MQ=255)

                     ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCa                                          >  1:1822991/1‑76 (MQ=255)

                     ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCAGGCa                                          >  2:682838/1‑76 (MQ=255)

                      tctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGCAGAAGTCCGCAAGCGCTGGCAg                                         >  2:4862833/1‑76 (MQ=255)

                      tctcGACTTCGATGATCTGATTTTATTGCCGACGCTGCTGCTGCAACGCAATGAAGAACTCCTCAACCGCTGGccg                                         >  2:518000/1‑74 (MQ=255)

                       ctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                        >  1:3872184/1‑76 (MQ=255)

                       ctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                        >  1:1503970/1‑76 (MQ=255)

                       ctcGACTTCGATGATCGGATTTTATTGCCGGCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                        <  2:568469/76‑1 (MQ=255)

                         cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCCACGCAATGAAGAAGTCCCCACGCGCTTGCACAAc                                      >  2:1591931/1‑76 (MQ=255)

                         cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAAc                                      >  1:2456552/1‑76 (MQ=255)

                           aCTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaa                                    >  2:2783970/1‑76 (MQ=255)

                            cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaa                                   >  1:2555515/1‑76 (MQ=255)

                              tCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAt                                 >  1:3381686/1‑76 (MQ=255)

                               cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                <  2:1969997/76‑1 (MQ=255)

                                gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                               <  2:576317/76‑1 (MQ=255)

                                gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                               <  2:4285094/76‑1 (MQ=255)

                                gatgatCTGACTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                               >  2:2276244/1‑76 (MQ=255)

                                 atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCg                              <  2:3299434/76‑1 (MQ=255)

                                 atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCg                              >  2:1525139/1‑76 (MQ=255)

                                  tgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACCAAATTCGc                             >  1:555912/1‑76 (MQ=255)

                                  tgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                             <  1:2621865/76‑1 (MQ=255)

                                  tgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                             >  1:3023448/1‑76 (MQ=255)

                                   gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCt                            >  1:4137146/1‑76 (MQ=255)

                                    atCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTa                           <  1:3570251/76‑1 (MQ=255)

                                     tCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGGCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                          >  2:3242346/1‑76 (MQ=255)

                                      cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                         <  1:1412804/76‑1 (MQ=255)

                                       tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATct                        >  2:1857561/1‑76 (MQ=255)

                                        gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctg                       <  1:1407604/76‑1 (MQ=255)

                                        gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctg                       <  1:824896/76‑1 (MQ=255)

                                          ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                     <  1:74854/76‑1 (MQ=255)

                                          ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                     >  2:2552128/1‑76 (MQ=255)

                                           tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgctg                    <  2:3262931/76‑1 (MQ=255)

                                            ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtgg                   <  2:597261/76‑1 (MQ=255)

                                            ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtgg                   >  1:3526874/1‑76 (MQ=255)

                                             tATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggt                  <  1:4593993/76‑1 (MQ=255)

                                             tATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggt                  <  2:3486407/76‑1 (MQ=255)

                                              aTTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGAGCTGGCAGAATAAAATTCGCTATCTGCtggtg                 >  1:513712/1‑76 (MQ=255)

                                               ttGCCGACGTTGCTGCTGCAGCGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                <  1:4417631/76‑1 (MQ=255)

                                               ttGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                >  2:2167168/1‑76 (MQ=255)

                                               ttGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                >  1:3665902/1‑76 (MQ=255)

                                                tGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa               >  2:1376879/1‑76 (MQ=255)

                                                tGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa               >  2:1598165/1‑76 (MQ=255)

                                                 gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGAt              <  2:1031003/76‑1 (MQ=255)

                                                 gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGAt              >  1:2141965/1‑76 (MQ=255)

                                                  ccGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg             <  2:2670064/76‑1 (MQ=255)

                                                  ccGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg             >  2:4720394/1‑76 (MQ=255)

                                                  ccGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg             >  1:3391831/1‑76 (MQ=255)

                                                    gtcGTTGCTGCTGCAACGTGCTGATGCAGTCCTCAAGGGCTGGCAGTACAGAATTCTCTGTCTGCTGGTGGATGAg           <  2:1466906/74‑1 (MQ=255)

                                                    gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg           >  2:1745928/1‑76 (MQ=255)

                                                    gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg           >  2:1251263/1‑76 (MQ=255)

                                                     aCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt          <  1:3946991/76‑1 (MQ=255)

                                                     aCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGt          >  1:1662946/1‑76 (MQ=255)

                                                      cGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTa         >  1:725741/1‑76 (MQ=255)

                                                       gTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGGGTAt        >  1:251403/1‑76 (MQ=255)

                                                       gTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt        >  1:2228516/1‑76 (MQ=255)

                                                       gTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt        >  2:2509494/1‑76 (MQ=255)

                                                       gTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt        >  2:2928481/1‑76 (MQ=255)

                                                         tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa      <  2:3763886/76‑1 (MQ=255)

                                                         tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa      <  2:4171884/76‑1 (MQ=255)

                                                         tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa      <  1:2931926/76‑1 (MQ=255)

                                                          gctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg     >  2:3479070/1‑76 (MQ=255)

                                                            tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa   >  2:3156094/1‑76 (MQ=255)

                                                             gctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac  <  1:2040959/76‑1 (MQ=255)

                                                             gctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac  <  2:2059461/76‑1 (MQ=255)

||

CACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGAC  >  NC_000913/3958809‑3958945

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`