BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,907,399	C→T	E13K (GAG→AAG)	cspC ←	stress protein, member of the CspA family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,907,399	0	C	T	100.0%	476.4 / NA	120	E13K (GAG→AAG)	cspC	stress protein, member of the CspA family
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (64/56); total (64/56)

Only 100 of 121 total aligned reads displayed.

ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907352‑1907445

aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt                                              <  1:2254794/50‑1 (MQ=255)

 tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             >  1:4661695/1‑50 (MQ=255)

  cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             >  1:3148176/1‑49 (MQ=255)

  cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             >  1:4581858/1‑49 (MQ=255)

  cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             >  1:4123581/1‑49 (MQ=255)

     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg                                            <  1:611705/47‑1 (MQ=255)

     tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                         >  1:362935/1‑50 (MQ=255)

      gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                         >  1:3847043/1‑49 (MQ=255)

       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                        <  1:3055698/49‑1 (MQ=255)

       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                        <  1:2808726/49‑1 (MQ=255)

       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                       >  1:465753/1‑50 (MQ=255)

       ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa                                       <  1:3002409/50‑1 (MQ=255)

       ctgcCATCAGCCGGAGTAAAGAAGCCAAAACCTTTAGACTTGTTGAAcc                                        <  1:1660855/49‑1 (MQ=37)

        tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  1:192348/1‑50 (MQ=255)

        tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  1:81196/1‑50 (MQ=255)

        tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  1:1443683/1‑50 (MQ=255)

        tgcCATCAGCCGGAGGAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc                                      >  1:4006169/1‑50 (MQ=38)

         gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     >  1:1670479/1‑50 (MQ=255)

          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     <  1:868348/49‑1 (MQ=255)

          ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACAACtt                                    <  1:64063/50‑1 (MQ=37)

           cATCAGCCGGAGTAATGAAGCCAAAATCTTTAGACTTGTTGAACCACtt                                    >  1:3796991/1‑49 (MQ=38)

           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    >  1:4699962/1‑49 (MQ=255)

           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    >  1:1613071/1‑49 (MQ=255)

           cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                   <  1:3439037/50‑1 (MQ=255)

            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  <  1:3695006/50‑1 (MQ=255)

            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  >  1:4127143/1‑50 (MQ=255)

            aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  <  1:4334807/50‑1 (MQ=255)

             tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAc                                 >  1:1678035/1‑50 (MQ=255)

             tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAc                                 >  1:3420586/1‑50 (MQ=255)

              cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                >  1:4085910/1‑50 (MQ=255)

               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc                                <  1:2489919/49‑1 (MQ=255)

               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:3625975/1‑50 (MQ=255)

               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:104365/1‑50 (MQ=255)

               aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:2080975/1‑50 (MQ=255)

                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:868972/1‑49 (MQ=255)

                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:4357584/1‑49 (MQ=255)

                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                              <  1:3124377/50‑1 (MQ=255)

                gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                              <  1:3536051/50‑1 (MQ=255)

                 ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                             <  1:3669222/50‑1 (MQ=255)

                  cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            <  1:1201600/50‑1 (MQ=255)

                  cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            >  1:4505019/1‑50 (MQ=255)

                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           >  1:6008478/1‑50 (MQ=255)

                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           >  1:5724592/1‑50 (MQ=255)

                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           <  1:2684317/50‑1 (MQ=255)

                   ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc                           <  1:3619419/50‑1 (MQ=255)

                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          >  1:1257927/1‑50 (MQ=255)

                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:1080080/50‑1 (MQ=255)

                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          >  1:126505/1‑50 (MQ=255)

                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:4098517/50‑1 (MQ=255)

                    gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:2626192/50‑1 (MQ=255)

                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         <  1:836028/50‑1 (MQ=255)

                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         <  1:855700/50‑1 (MQ=255)

                     aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:494347/1‑50 (MQ=255)

                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        <  1:2989506/50‑1 (MQ=255)

                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        <  1:704618/50‑1 (MQ=255)

                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:2056734/1‑49 (MQ=255)

                      gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:726229/1‑49 (MQ=255)

                       tAATGAAGCCAAAACCTTTAGACTTGTTGAGCCACTTAACCTGACCtttg                       <  1:683491/50‑2 (MQ=37)

                       tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        >  1:2098226/1‑49 (MQ=255)

                        aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:1049544/1‑50 (MQ=255)

                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      <  1:1029868/49‑1 (MQ=255)

                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      <  1:2452056/49‑1 (MQ=255)

                         aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt                     <  1:4077331/50‑1 (MQ=255)

                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                    >  1:584878/1‑50 (MQ=255)

                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                    >  1:548131/1‑50 (MQ=255)

                          tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc                    >  1:4349740/1‑50 (MQ=255)

                           gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   >  1:5873645/1‑50 (MQ=255)

                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  >  1:626948/1‑50 (MQ=255)

                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  <  1:2786752/50‑1 (MQ=255)

                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  <  1:3418878/50‑1 (MQ=255)

                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   <  1:596188/49‑1 (MQ=255)

                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   <  1:3808638/49‑1 (MQ=255)

                            aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   <  1:3643593/49‑1 (MQ=255)

                             aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 >  1:1915445/1‑50 (MQ=255)

                              gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                >  1:1136965/1‑50 (MQ=255)

                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:5276322/1‑50 (MQ=255)

                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:378132/1‑50 (MQ=255)

                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               <  1:4805869/50‑1 (MQ=255)

                               ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               <  1:5580373/50‑1 (MQ=255)

                                cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc              <  1:5891597/50‑1 (MQ=255)

                                cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:3659465/1‑49 (MQ=255)

                                cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:3689257/1‑49 (MQ=255)

                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:95441/1‑50 (MQ=255)

                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             <  1:2218433/50‑1 (MQ=255)

                                 aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa             >  1:4895608/1‑50 (MQ=255)

                                  aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt            <  1:2828712/50‑1 (MQ=255)

                                   aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           <  1:4462099/50‑1 (MQ=255)

                                   aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           >  1:3443021/1‑50 (MQ=255)

                                    aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          <  1:4948406/50‑1 (MQ=255)

                                    aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          >  1:5226172/1‑50 (MQ=255)

                                    aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          <  1:5706677/50‑1 (MQ=255)

                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          >  1:4518825/1‑49 (MQ=255)

                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         <  1:5670383/50‑1 (MQ=255)

                                     ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         <  1:1489589/50‑1 (MQ=255)

                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg         >  1:1321270/1‑49 (MQ=255)

                                      cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        <  1:5811376/50‑1 (MQ=255)

                                         tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa      <  1:3450815/49‑1 (MQ=255)

                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa    <  1:5335401/49‑1 (MQ=255)

                                           gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt   <  1:624050/50‑1 (MQ=255)

                                            aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  >  1:571088/1‑50 (MQ=255)

ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907352‑1907445

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