BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,756,138	Δ1 bp	intergenic (+183/‑21)	ssrA → / → intA	tmRNA/CP4‑57 prophage; integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,756,137	0	C	.	100.0%	299.4 / NA	60	intergenic (+182/‑22)	ssrA/intA	tmRNA/CP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (27/33); total (27/33)

TTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTA  >  NC_000913/2756089‑2756182

ttAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                                               <  1:2425093/49‑1 (MQ=255)

 tAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                              >  1:4021676/1‑48 (MQ=255)

  aGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGcta                                             >  1:665411/1‑47 (MQ=255)

  aGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGcta                                             >  1:5653541/1‑47 (MQ=255)

    gCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAAc                                         <  1:3725599/50‑1 (MQ=255)

     cccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                        >  1:793951/1‑50 (MQ=255)

     cccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                        <  1:5757635/50‑1 (MQ=255)

      ccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                        >  1:3086512/1‑49 (MQ=255)

      ccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGc                                       <  1:5636357/50‑1 (MQ=255)

        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAg                                     >  1:126058/1‑50 (MQ=255)

        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAg                                     >  1:5911151/1‑50 (MQ=255)

        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAg                                     >  1:4552440/1‑50 (MQ=255)

         ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                    <  1:3901792/50‑1 (MQ=255)

         ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                    <  1:5257983/50‑1 (MQ=255)

         ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                    <  1:4648167/50‑1 (MQ=255)

         ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                    <  1:1354630/50‑1 (MQ=255)

         ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                    >  1:1833969/1‑50 (MQ=255)

           gATAGGCCCAACGAAAAGCTCTATTGTTTACTTTGGG‑CTAAACGCAgg                                    >  1:3220959/1‑48 (MQ=255)

           gATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                    >  1:322410/1‑48 (MQ=255)

           gATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                    >  1:6540139/1‑48 (MQ=255)

            aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGgag                                 <  1:5460617/50‑1 (MQ=255)

            aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGga                                  <  1:3435540/49‑1 (MQ=255)

            aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGga                                  <  1:6572779/49‑1 (MQ=255)

             tAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGgaga                                <  1:876926/50‑1 (MQ=255)

               ggCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGAc                               <  1:2207066/49‑1 (MQ=255)

               ggCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACt                              <  1:144562/50‑1 (MQ=255)

               ggCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACt                              <  1:5691046/50‑1 (MQ=255)

                gCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTc                             >  1:2666533/1‑50 (MQ=255)

                gCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTc                             <  1:3781355/50‑1 (MQ=255)

                 cccAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTcc                            >  1:205878/1‑50 (MQ=255)

                  ccAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTccc                           <  1:208306/50‑1 (MQ=255)

                   cAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTcccc                          <  1:1612320/50‑1 (MQ=255)

                   cAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTccc                           >  1:268536/1‑49 (MQ=255)

                    aaCGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCA                         <  1:5762875‑M1/50‑1 (MQ=255)

                    aaCGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCA                         >  1:5342871‑M1/1‑50 (MQ=255)

                    aaCGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCA                         >  1:4759907‑M1/1‑50 (MQ=255)

                     aCGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAc                        >  1:1678943‑M1/1‑49 (MQ=255)

                       gAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcct                      >  1:5329416‑M1/1‑47 (MQ=255)

                       gAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcct                      >  1:285086‑M1/1‑47 (MQ=255)

                       gAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcct                      >  1:3863651‑M1/1‑47 (MQ=255)

                        aaaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctt                     >  1:3441088‑M1/1‑46 (MQ=255)

                        aaaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctt                     <  1:4871773‑M1/50‑5 (MQ=255)

                        aaaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctt                     <  1:1647995‑M1/50‑5 (MQ=255)

                           aGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttga                  <  1:3929431‑M1/50‑8 (MQ=255)

                           aGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttg                   >  1:1081368‑M1/1‑43 (MQ=255)

                             ctctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaa                <  1:4743418‑M1/50‑10 (MQ=255)

                              tctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaa               >  1:3312649‑M1/1‑40 (MQ=255)

                              tctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaa               >  1:3490451‑M1/1‑40 (MQ=255)

                              tctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaa               >  1:3561393‑M1/1‑40 (MQ=255)

                                tATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaa               <  1:4051792‑M1/48‑11 (MQ=255)

                                 aTTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacag            >  1:5130114‑M1/1‑37 (MQ=255)

                                  ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacag            <  1:215941‑M1/49‑14 (MQ=255)

                                  ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacag            <  1:5971334‑M1/49‑14 (MQ=255)

                                  ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaaca             <  1:894970‑M1/48‑13 (MQ=255)

                                  ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaaca             <  1:5965453‑M1/48‑13 (MQ=255)

                                   tGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagga          >  1:3886204‑M1/1‑35 (MQ=255)

                                    gTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat         >  1:2877409‑M1/1‑34 (MQ=255)

                                     tttACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg        >  1:423275‑M1/1‑33 (MQ=255)

                                     tttACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat         <  1:3793458‑M1/49‑17 (MQ=255)

                                     tttACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat         <  1:4462481‑M1/49‑17 (MQ=255)

                                        aCGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgta      <  1:4543365‑M1/49‑20 (MQ=255)

                                         cGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtag     <  1:3201436‑M1/49‑21 (MQ=255)

                                          gTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcg   <  1:6424450‑M1/50‑23 (MQ=255)

                                           ttGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcga  <  1:1926970‑M1/50‑24 (MQ=255)

TTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTA  >  NC_000913/2756089‑2756182

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