BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,504,637	G→A	intergenic (+381/‑344)	insB‑25 → / → ECB_04162	b4563(b4576); IS1 protein insB/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,504,637	0	G	A	100.0%	406.7 / NA	109	intergenic (+381/‑344)	insB‑25/ECB_04162	b4563(b4576); IS1 protein insB/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (64/45); total (64/45)

Only 100 of 110 total aligned reads displayed.

AGCAAAATCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCGGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATGGCCGT  >  REL606/4504564‑4504708

aGCAAAATCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGt                                                                       >  2:1114851/1‑76 (MQ=255)

   aaaaTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGtctc                                                                    >  1:399624/1‑76 (MQ=255)

   aaaaTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGtctc                                                                    >  2:1100522/1‑76 (MQ=255)

    aaaTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATGGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCa                                                                   <  2:1310410/76‑1 (MQ=255)

    aaaTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCa                                                                   <  2:1034047/76‑1 (MQ=255)

    aaaTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCa                                                                   <  1:280800/76‑1 (MQ=255)

    aaaTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCa                                                                   >  2:590409/1‑76 (MQ=255)

    aaaTCAGCAAACGTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCa                                                                   <  1:731246/76‑1 (MQ=255)

     aaTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAg                                                                  <  1:1518413/76‑1 (MQ=255)

      aTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGt                                                                 <  1:1338087/76‑1 (MQ=255)

      aTCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTGGGTGTGGCAGTCTCAGt                                                                 >  1:166327/1‑76 (MQ=255)

       tCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTc                                                                <  2:17271/76‑1 (MQ=255)

         aGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCg                                                              >  2:88000/1‑76 (MQ=255)

            aaaCTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCANTCTCAGTCCGGaa                                                           >  2:1400120/1‑76 (MQ=255)

            aaaCTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGaa                                                           >  2:842456/1‑76 (MQ=255)

            aaaCTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGaa                                                           >  2:23124/1‑76 (MQ=255)

             aaCTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAAc                                                          <  1:1275424/76‑1 (MQ=255)

             aaCTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAAc                                                          <  2:531048/76‑1 (MQ=255)

              acttacttAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACa                                                         <  1:1495393/76‑1 (MQ=255)

              acttacttAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACa                                                         <  2:1329178/76‑1 (MQ=255)

              acttacttAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACa                                                         <  1:219102/76‑1 (MQ=255)

               cttacttaACTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTTCGGAACaa                                                        <  1:58935/76‑1 (MQ=21)

                ttacttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAAt                                                       >  1:1130545/1‑76 (MQ=255)

                ttacttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAAt                                                       >  1:601191/1‑76 (MQ=255)

                 tacttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAAtt                                                      >  2:502506/1‑76 (MQ=255)

                 tacttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAAtt                                                      >  1:1235367/1‑76 (MQ=255)

                 tacttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAAtt                                                      >  1:1200578/1‑76 (MQ=255)

                 tacttaATTATATCGAGTGCGATTACAATCGCTGGCCATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAAtt                                                      >  2:720718/1‑76 (MQ=255)

                  acttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAAttt                                                     <  2:1236699/76‑1 (MQ=255)

                   cttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCTGAACAATTTg                                                    <  1:1417989/76‑1 (MQ=255)

                    ttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGa                                                   >  2:970272/1‑76 (MQ=255)

                    ttaATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGa                                                   >  1:1232464/1‑76 (MQ=255)

                    ttaACTCTATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGa                                                   <  1:1505355/76‑1 (MQ=18)

                     taATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGaa                                                  >  2:694624/1‑76 (MQ=255)

                      aaTTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGaaa                                                 >  2:388386/1‑76 (MQ=255)

                       aTTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGaaaa                                                <  1:970273/76‑1 (MQ=255)

                        ttATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAAt                                               >  2:828328/1‑76 (MQ=255)

                        ttATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGACAGTCTCAGTCCGGAACAATTTGAAAAt                                               <  2:657548/76‑1 (MQ=255)

                         tataTCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATc                                              <  2:1481048/76‑1 (MQ=255)

                         tataTCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATc                                              >  2:949223/1‑76 (MQ=255)

                         tataTCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATc                                              <  1:794856/76‑1 (MQ=255)

                         tataTCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATc                                              >  2:789985/1‑76 (MQ=255)

                         tataTCGAGTGCGATTACAATCGCTGGAGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTCAGAATc                                              >  1:1519093/1‑76 (MQ=18)

                          atatCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGCAAATCa                                             >  2:1494830/1‑76 (MQ=255)

                          atatCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCa                                             <  2:66856/76‑1 (MQ=255)

                           tatCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAg                                            >  2:595508/1‑76 (MQ=255)

                           tatCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAg                                            >  2:100209/1‑76 (MQ=255)

                            atCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGa                                           <  1:543709/76‑1 (MQ=255)

                             tCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGaa                                          <  2:674944/76‑1 (MQ=255)

                              cGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGNACAATTTGAAAATCAGAAt                                         <  1:1423630/76‑1 (MQ=255)

                              cGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAAt                                         >  2:325700/1‑76 (MQ=255)

                              cGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAAt                                         >  2:517400/1‑76 (MQ=255)

                              cGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAAt                                         >  1:1291611/1‑76 (MQ=255)

                               gAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAAtc                                        >  2:1332787/1‑76 (MQ=255)

                                 gTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAAtctc                                      <  2:537971/76‑1 (MQ=255)

                                  tGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCa                                     <  1:1130476/76‑1 (MQ=255)

                                   gCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCAc                                    >  2:419929/1‑76 (MQ=255)

                                    cGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACt                                   <  1:1095009/76‑1 (MQ=255)

                                    cGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACt                                   <  1:593442/76‑1 (MQ=255)

                                     gATTGCAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGATCAATTTGAAAATCTGAATTTCANtt                                  >  2:1448752/1‑76 (MQ=18)

                                     gATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACtt                                  >  2:464742/1‑76 (MQ=255)

                                     gATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACtt                                  >  2:295346/1‑76 (MQ=255)

                                      aTTACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTa                                 >  1:536702/1‑76 (MQ=255)

                                      aTTACAATAGCTGGCGNTCGTACAGTTGGTCTNGCAGTCTGAGTCCGGAACAATTTGAAAATCAGNATCTCANTNa                                 >  1:112400/1‑76 (MQ=18)

                                        tACAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTaaa                               <  1:396168/76‑1 (MQ=255)

                                         aCAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTaaaa                              >  1:836877/1‑76 (MQ=255)

                                         aCAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTaaaa                              >  2:1170180/1‑76 (MQ=255)

                                         aCAATCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTAAAAATCAGAATCTCACTTaaaa                              >  1:836879/1‑76 (MQ=255)

                                           aaTCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACg                            <  1:311764/76‑1 (MQ=255)

                                             tCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGc                          >  1:564322/1‑76 (MQ=255)

                                             tCGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGc                          <  2:176195/76‑1 (MQ=255)

                                             tCGCTGGAGATGGTACAGTTGGTGTGGCAGGCTCAGTCCGGCACAATTTGAAAATCAGAATCTCACTTAAAACGGc                          >  2:934737/1‑76 (MQ=21)

                                              cGCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCa                         <  2:581268/76‑1 (MQ=255)

                                               gCTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCAt                        <  2:596993/76‑1 (MQ=255)

                                                cTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATc                       >  2:454801/1‑76 (MQ=255)

                                                cTGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATc                       >  1:236760/1‑76 (MQ=255)

                                                 tGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATcc                      >  2:267797/1‑76 (MQ=255)

                                                 tGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATcc                      >  2:267793/1‑76 (MQ=255)

                                                 tGGCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATcc                      >  2:267805/1‑76 (MQ=255)

                                                 tGGCGATGGNACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATcc                      >  1:40295/1‑76 (MQ=255)

                                                  ggCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCa                     >  2:851849/1‑76 (MQ=255)

                                                   gCGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCAt                    <  2:526771/76‑1 (MQ=255)

                                                    cGATGGTACAGTTGGTGTGGCAGTCTCAGTCCNGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCNCtt                   >  2:1421959/1‑76 (MQ=25)

                                                    cGATGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCAtt                   >  1:998884/1‑76 (MQ=255)

                                                       tGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTa                <  2:20660/76‑1 (MQ=255)

                                                       tGGTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTa                <  1:1102101/76‑1 (MQ=255)

                                                        ggTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTAc               <  1:806570/76‑1 (MQ=255)

                                                         gTACAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACa              >  2:627340/1‑76 (MQ=255)

                                                           aCAGTTGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACaaa            <  1:1332788/76‑1 (MQ=255)

                                                           aCAGCTGGTGTGGCAGTATCAGTCCGGAACAATTTGAAAATCAGAATCTCGCTTAAAACGGCATGCATTATACnaa            >  2:1439342/1‑76 (MQ=18)

                                                            cAGTTGCTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACaaaa           <  1:1395691/76‑1 (MQ=255)

                                                              gTTGGTGTGGCAGTCTCAATCAGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGa         <  1:1437660/76‑1 (MQ=37)

                                                               ttGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGAt        <  1:222709/76‑1 (MQ=255)

                                                                tGGTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATg       >  1:638209/1‑76 (MQ=255)

                                                                 ggTGTGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATgg      >  2:546446/1‑76 (MQ=255)

                                                                  gtgtGGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATGGc     >  1:154586/1‑76 (MQ=255)

                                                                   tgtgGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATGGcc    >  2:1330363/1‑76 (MQ=255)

                                                                    gtgGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATGGCCg   <  2:161571/76‑1 (MQ=255)

                                                                     tgGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATGGCCGt  >  2:790017/1‑76 (MQ=255)

                                                                     tgGCAGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATGGCCGt  >  2:759983/1‑76 (MQ=255)

AGCAAAATCAGCAAACTTACTTAATTATATCGAGTGCGATTACAATCGCTGGCGATGGTACAGTTGGTGTGGCGGTCTCAGTCCGGAACAATTTGAAAATCAGAATCTCACTTAAAACGGCATCCATTTTACAAAAGATGGCCGT  >  REL606/4504564‑4504708

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