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%	348.2 / NA	69	intergenic (+182/‑22)	ssrA/intA	tmRNA/CP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (31/38); total (31/38)

TTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTT  >  NC_000913/2756089‑2756181

ttAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                                              >  1:3540720/1‑49 (MQ=255)

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

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

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

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

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

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

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

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

    gCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAAc                                        >  1:4053334/1‑50 (MQ=255)

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

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

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

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

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

       cATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                       >  1:3871486/1‑48 (MQ=255)

       cATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                       >  1:6988827/1‑48 (MQ=255)

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

        aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAg                                    <  1:3477286/50‑1 (MQ=255)

          tGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                   >  1:3913761/1‑49 (MQ=255)

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

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

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

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

             tAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGgaga                               >  1:5146339/1‑50 (MQ=255)

             tAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGgaga                               >  1:851751/1‑50 (MQ=255)

             tAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGgaga                               >  1:2077515/1‑50 (MQ=255)

              aGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGAc                              >  1:3207250/1‑50 (MQ=255)

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

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

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

                  ccAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTccc                          >  1:7027314/1‑50 (MQ=255)

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

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

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

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

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

                   cAACGAAAAGCTCTATTGTTGACGTTGGG‑CTAAACGCAGGGAGACTcccc                         <  1:1108423/50‑1 (MQ=255)

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

                     aCGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAc                       <  1:354608‑M1/50‑2 (MQ=255)

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

                      cGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcc                      <  1:4775738‑M1/50‑3 (MQ=255)

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

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

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

                         aaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAccttt                   <  1:2633908‑M1/50‑6 (MQ=255)

                         aaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAccttt                   >  1:5492408‑M1/1‑45 (MQ=255)

                         aaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAccttt                   <  1:261280‑M1/50‑6 (MQ=255)

                         aaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAccttt                   <  1:619703‑M1/50‑6 (MQ=255)

                         aaaGCTCTATTGATTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAccttt                   <  1:7307705‑M1/50‑6 (MQ=255)

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

                            gCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaa                <  1:6745377‑M1/50‑9 (MQ=255)

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

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

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

                                  ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagg          <  1:6008071‑M1/50‑15 (MQ=255)

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

                                  ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacag           >  1:3369407‑M1/1‑36 (MQ=255)

                                  ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacag           >  1:6940793‑M1/1‑36 (MQ=255)

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

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

                                  ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaac             <  1:4455397‑M1/47‑12 (MQ=255)

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

                                    gTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagga         <  1:2135677‑M1/49‑16 (MQ=255)

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

                                       tACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg       >  1:3715142‑M1/1‑31 (MQ=255)

                                       tACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg       >  1:6436045‑M1/1‑31 (MQ=255)

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

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

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

                                         cGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc   <  1:1165538‑M1/50‑22 (MQ=255)

                                         cGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc   <  1:6311017‑M1/50‑22 (MQ=255)

                                         cGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc   >  1:6137964‑M1/1‑29 (MQ=255)

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

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

                                          gTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc   >  1:4814675‑M1/1‑28 (MQ=255)

TTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTT  >  NC_000913/2756089‑2756181

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